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ABSTRACT 


This  is  an  experimental  investigation  of  the  filaments  and  drops  generated  at  the 
free  surface  of  liquid  wall  jets  formed  over  smooth-  and  sand-roughened  surfaces.  The 
jet  characteristics  and  the  geometric  properties  of  the  filaments  and  drops  were  measured 
from  images  captured  using  high-speed  digital  cameras.  A  statistical  investigation  of  the 
various  properties  revealed  the  characteristic  behavior  of  the  filaments  and  drops  as  a 
function  of  the  relative  wall  roughness,  wall  curvature  and  jet  inertia.  For  this 
investigation,  the  wall  jet  Reynolds  number  ranged  from  2.6x10"^  to  4.5x10'*,  the  Froude 
number  from  19  to  33  and  the  Weber  number  from  1600  to  4700. 

The  emphasis  herein  was  on  the  physics  of  the  process  rather  than  the 
development  of  empirical  relationships.  As  such,  the  results  indicate  that  spray 
generation  from  a  wall  jet  is  a  boundary-layer-driven  phenomenon,  requiring  that  the  jet 
be  in  a  highly  supercritical  state  (Fr  »  1).  Wall  roughness  reduces  the  minimum 
necessary  level  of  supercriticality,  but  it  is  not  a  prerequisite  condition  for  the  formation 
of  drops.  While  increasing  the  jet  inertia  enhances  the  drop  formation  process,  concave 
wall  curvature  tends  to  reduce  the  quantity  and  the  energy  of  the  drop  forming  events. 
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INTRODUCTION 


A.  BACKGROUND 

How  is  it  that  something  as  ubiquitous  as  a  drop  of  water  could  be  so 
intellectually  enticing?  The  seemingly  elementary  process  of  a  drop  releasing  from  a 
fixed  object  and  falling  in  a  quiescent  medium  presents  a  myriad  of  formidable  fluid 
dynamic  challenges  even  for  the  most  capable  researcher.  As  a  simple  example,  consider 
the  following  three  questions: 

•  What  influence  does  the  fluid  circulating  within  the  drop  have  on  the  events 
that  will  transpire  and  vice  versa? 

•  At  what  point  in  the  ever-shrinking  neck  of  the  drop  do  the  Navier-Stokes 
equations  become  singular? 

•  When  does  the  dynamic  variation  of  surface  tension  become  a  concern? 

At  present,  a  complete  analytical  description  remains  elusive,  while  computer  models  that 
can  accurately  reproduce  this  simple  laminar  hydrodynamic  marvel  require  a  good  deal  of 
numerical  ingenuity. 

Now  take  this  same  fluid  particle,  join  it  with  a  continuum  of  others,  and  direct  it 
into  a  similar  environment.  As  it  flows  out  to  its  ultimate  destiny  (whatever  that  may  be), 
imbue  it  with  random  velocity  fluctuations,  expose  it  to  shear  forces,  and  perhaps  even 
accelerate  it,  before  once  again  leaving  it  to  its  own  demise.  Such  events  transpire  an 
unaccountable  number  of  times  per  day  as  individuals  drive  their  cars  to  work,  fight  forest 
fires,  or  just  turn  on  faucets.  Obviously,  there  are  many  challenges  associated  with 
understanding  how  a  drop  forms  in  such  harsh  and  stimulating  environments.  In  most 
instances,  the  physical  parameters  of  the  problem  have  such  complex  and  intertwined 
relationships  that  the  engineer  must  rely  on  correlations  that  can  guide  him  within  a 
limited  range.  If  the  individual  ventures  outside  the  range,  then  new  correlations  must  be 
established.  Nonetheless,  many  in-roads  have  been  made  which  enable  us  to  see  the 
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problems  more  clearly,  even  if  they  cannot  be  completely  understood.  For  example,  the 
description  of  the  breakup  of  a  liquid  column  by  Lord  Rayleigh  nearly  100  years  ago  was 
fundamental  to  the  principles  later  developed  for  the  design  of  ink-jet  printers. 

The  investigation  of  drops  (or  spray)  discussed  in  this  writing  is  directed  towards 
understanding  how  drops  are  produced  at  the  free  surface  of  a  sheet  of  water  flowing  over 
a  wall;  that  is,  a  high-speed  liquid  wall  jet.  The  motivation  for  the  investigation  is  the 
desire  to  understand  the  physics  behind  the  spray  created  by  the  bow  of  a  ship  as  it  moves 
through  the  water.  The  bow  forms  bow  sheets  (or  spray  roots)  along  the  hull  that  produce 
drops  at  their  free  surfaces,  and  also  detach  from  the  hull  forming  thin  liquid  sheets, 
which  then  disintegrate  into  spray.  The  consequence  of  this  process  is  that  much  of  this 
water  spray  is  deposited  on  the  superstructure  of  the  ship.  At  a  minimum,  this  frequently 
results  in  treacherous  topside  footing  and  reduced  equipment  operability  or  damage.  Li 
more  severe  cases,  when  accompanied  by  subfreezing  temperatures,  it  can  lead  to  a 
reduction  in  ship  stability  as  the  frozen  spray  builds  up  on  the  superstructure,  raising  the 
ship’s  center  of  gravity  and,  reducing  its  righting  moment.  Li  the  case  of  warships,  there 
is  the  added  concern  that  the  wetted  surfaces  make  the  vessel  more  susceptible  to 
electromagnetic  detection,  which  may  limit  the  ship’s  ability  to  fulfill  its  designated 
mission.  Outside  of  this  specific  interest,  there  are  many  other  practical  applications  for 
understanding  how  a  continuum  becomes  a  multi-phase  flow  via  a  drop  formation 
process:  agriculture,  spillway  hydraulics,  chemical  warfare  and  food  processing,  to  name 
a  few.  As  the  first  exploration  of  this  type  into  this  process,  the  investigation  is  directed 
towards  acquiring  a  fundamental  understanding  of  the  spray  formation  process,  rather 
than  obtaining  correlations  to  describe  the  gross  spray  features. 

B.  LITERATURE  REVIEW 

Despite  the  clear  application  and  intellectual  challenges  of  the  high-speed  liquid 
wall  jet,  wall  jets  do  not  appear  to  have  been  studied  before.  Only  a  number  of  low 
Reynolds  and  Froude  number  investigations  can  be  found  in  the  literature.  Hanratty  and 
Engen  (1957)  examined  the  effects  of  a  turbulent  gas  stream  blowing  over  fairly  thin  (4.1 
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to  5.3  mm)  slow  moving  (11.9  <  Re^  <  508)  smooth-walled  water  channel  flows.  Their 
findings  show  that  the  free  surface  of  a  shallow  water  channel  is  more  difficult  to  disturb 
than  a  deeper  one.  They  also  found  that  the  free  surface  characteristics  of  the  water  are 
more  strongly  influenced  by  the  Reynolds  number  of  the  water  than  that  of  the  air,  but  a 
later  study  of  flow  down  an  inclined  channel  suggests  otherwise  (Roy  and  Jain  1989). 
For  the  deepest  channel  combined  with  the  highest  gas  velocity,  the  maximum  surface 
distortions  noted  by  Hanratty  and  Engen  were  approximately  0.75  mm,  which  are  much 
smaller  than  those  seen  here. 

Finley  et  al  (1966)  also  investigated  the  mean  flow  properties  of  water  flow  in  a 
smooth-walled  shallow  channel  arrangement.  Their  goal  was  to  determine  an  empirical 
equation  for  the  mean  streamwise  velocity  in  the  “fully  developed  region  of  the  flow,” 
which  they  define  as  the  point  where  the  free  surface  is  disturbed  by  the  turbulent 
boundary  layer.  The  Reynolds  numbers  of  the  flows  they  investigated  were  similar  to 
those  encountered  here  (Re  =  3x10'*),  but  the  Froude  numbers  (2  <  Fr  <  9)  were 
substantially  lower.  Using  Preston  Tubes,  Finley  et  al  measured  the  flow  to  within  0.76 
mm  of  the  free  surface  of  a  8.5  mm  deep  channel,  and  found  that  the  mean  flow  in  the 
buffer  and  inner  region  could  be  described  by  a  power  law  relationship. 
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where, 

8.02  <C<  8.14  6.5  <n<  6.61. 

In  the  outer  region  they  claim  that  the  velocity  profiles  followed  the  log-law, 
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They  pointed  out  that  this  form  of  the  pol3mioniial  appears  to  be  independent  of  the 
Reynolds  number  over  the  range  of  values  they  have  examined.  However,  as  was  the 
case  with  Hanratty  and  Engen  (1957),  the  work  of  Finley  et  al  is  limited  to  flows  over 
smooth  surfaces  with  minimal  free  surface  distortions. 

In  an  investigation  of  self-aerated  flows,  Volkart  (1980)  was  the  first  to  note  the 
drop-producing  liquid  jet  structures  at  the  free  surface  of  an  open  channel  flow.  He  did 
not  examine  the  process  by  which  the  structures  and  their  subsequent  drops  were 
produced,  but  instead  derived  an  expression  describing  the  drop  size  and  velocity 
necessary  for  a  drop  to  entrain  air  as  it  impacted  the  free  surface.  More  recently.  Rein 
(1998)  used  an  energy  balance  method  to  discuss  the  same  problem.  Rein  points  out  that 
Volkart's  (1980)  more  empirical  approach  did  not  adequately  explain  the  aeration 
mechanism  that  led  to  the  formation  of  large  bubbles  in  the  flow,  but  he  also  admits  that 
his  analytical  solution  failed  to  predict  the  small  droplet  aeration  witnessed  in  the 
laboratory  and  in  shallow  channels.  Brattberg  et  al  (1998)  have  performed  experimental 
work  in  supercritical  (Fr  =  9)  open  channel  flows.  They,  too,  were  primarily  concerned 
with  examining  the  degree  of  air  entrainment  within  the  Jet,  and  make  no  mention  of  the 
free  surface  characteristics.  Apparently  all  three  works  assumed  that  the  drop  formation 
is  only  a  result  of  turbulence  existing  near  the  free  surface,  but  they  did  not  discuss  the 
influence  of  the  channel  wall  and  its  boundary  layer  on  the  production  of  this  turbulence. 
Rein  recognizes  that  the  boundary  layer  is  somehow  involved,  but  concedes  that  there  is  a 
general  lack  of  knowledge  of  turbulent  free-surface  flows  in  the  high  Froude  number 
regime  (F  »  1). 

Others  have  also  examined  flows  of  this  nature,  but  their  results  do  not  further 
elucidate  the  physics  important  to  this  study.  For  example,  they  were  concerned  with 
flows  having  much  lower  Reynolds  numbers  than  those  cited  above  (Binnie  1957, 
Benjamin  1957),  extremely  thin  films  (Knuth  1955),  or  considered  instances  where  the 
liquid  flow  is  driven  with  a  large  shear  imposed  by  a  moving  gas  above  the  liquid  (van 
Rossum  1959,  Miles  1960,  Craik  1966). 
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1.  Free  Jets 

Previous  examinations  of  the  factors  that  contribute  to  liquids  breaking  up  into 
drops  are  mostly  found  in  free  jet  studies.  Although  the  mechanisms  involved  are 
somewhat  different  than  those  encountered  here,  much  of  the  physics  noted  in  these 
efforts  is  still  relevant  to  the  actions  which  occur  at  the  free  surface  of  a  liquid  wall  jet. 
As  it  is  encountered  in  most  studies,  the  term  “free  jet”  is  understood  to  imply  that  the  jet 
is  issuing  into  an  environment  with  a  substantially  lower  density  than  that  of  itself.  In 
most  cases  it  is  also  axisymmetric  at  the  point  of  origin.  One  of  the  earliest  analytical 
works  examining  the  breakup  of  such  a  liquid  jet  is  the  linear  analysis  of  Rayleigh  (1945). 
He  found  that  axisymmetric  infinitesimal  disturbances  with  wavelengths  greater  than  the 
jet  circumference  are  amplified,  and  that  disturbances  with  wavelengths  approximately 
equal  to  4.5  Ddr  have  the  maximum  rate  of  amplification.  Weber  (1931)  and  Ohnesorge 
(1936)  expanded  his  results  to  include  the  effects  of  viscosity,  which  Sterling  and 
Sleicher  (1975)  then  modified  to  achieve  better  agreement  with  the  reported  experimental 
data  in  the  nonlinear  region  of  jet  breakup.  Yuen  (1968)  and  Neyfeh  (1970)  both 
developed  a  complex  nonlinear  analytical  analysis  to  account  for  finite-sized 
disturbances.  Each  of  the  above  efforts  met  with  a  certain  degree  of  success,  producing 
some  marginal  improvement  over  the  simple  analysis  of  Lord  Rayleigh. 

Others  have  used  experimental  results  to  describe  the  breakup  process  in  terms  of 
the  Weber,  Froude,  Ohnesorge  and/or  Reynolds  number  of  the  jet.  For  example,  the 
following  is  one  of  the  more  popular  descriptions  by  Ranz  (1958)  and  Reitz  (1978),  and 
separates  the  modes  of  jet  breakup  into  four  regimes  that  are  functions  of  the  Weber 
number: 

•  Drip  mode:  Small  drops  form  at  the  nozzle  exit,  inertia  is  negligible,  and 
surface  tension  and  gravity  govern  the  breakup  of  the  jet. 

•  Rayleigh  mode:  A  low  speed  coherent  jet  exists,  but  its  velocity  is  low 
enough  that  the  shear  interaction  with  the  surrounding  gas  is  negligible. 
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Inertia  and  surface  tension  are  the  predominate  forces,  causing  the  jet  to 
disintegrate  into  drops  with  diameters  greater  than  that  of  the  undisturbed  jet. 

•  First  wind-induced  instability  mode:  The  higher  inertia  liquid  jet  develops 
large  instabilities  resulting  from  the  relative  motion  between  the  jet  and  the 
gas.  Surface  tension  acts  to  stabilize  the  jet  since  the  disturbance  wave 
numbers  produced  in  this  mode  are  greater  than  the  critical  value  (Ruff  et  al 
1989). 

•  Second  wind-induced  instability  mode:  As  described  by  Wu  et  al  (1991, 
1992),  small  drops  are  released  from  the  surface  of  the  liquid,  and  any  small- 
scale  disturbances  resulting  from  pre-existing  turbulence  serve  to  enhance  the 
drop  formation  process.  As  shown  in  Fig.  1  (Hoyt  and  Taylor,  1977a),  small- 
scale  surface  irregularities  appear  as  a  result  of  the  pre-existing  turbulence  in 
the  jet.  These  disturbances  can  produce  fluid  filaments  which  extend  radially 
outward,  and  frequently  break  up  into  drops  (where  Ddr «  Djet).  Both  Hoyt  and 
Taylor  and  Wu  et  al  (1992)  have  concluded  that  aerod3mamic  shear  and  drag 
do  not  appear  to  contribute  to  the  formation  of  these  drops.  Although  the 
main  body  of  the  jet  may  eventually  display  large-scale  irregularities,  the 
majority  of  the  disintegration  process  occurs  through  surface  disintegration 
(Ruff  et  al  1989).  For  both  wind-induced  regimes,  the  jet  length  becomes 
shorter,  and  the  disintegration  more  chaotic,  as  the  jet  velocity  increases. 

•  Atomization  mode:  The  jet  breaks  up  at  the  nozzle  exit,  forming  a  fine  spray 
of  drops.  The  mechanisms  leading  to  atomization  are  still  not  fully 
understood. 

The  velocity  profile  at  the  nozzle  exit  is  also  an  important  jet  breakup  parameter. 
As  in  the  case  of  pipe  flows,  the  length-to-diameter  ratio  of  the  nozzle  determines  the 
state  of  the  boundary  layer  and,  therefore,  the  velocity  profile  at  the  exit.  For  laminar  jets, 
the  greater  the  length-to-diameter  ratio,  the  more  developed  the  velocity  profile  is  at  the 
nozzle  exit.  In  studies  on  the  stability  of  Newtonian  jets,  Rupe  (1962)  and  Grant  and 
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Middleman  (1966)  found  that  the  velocity  profile  does  not  influence  the  breakup  length  of 
the  jet  at  low  Reynolds  numbers.  However,  beyond  a  certain  point  within  the  Rayleigh 
regime,  a  flat  velocity  profile  (i.e.,  slug  flow),  or  a  fully  developed  turbulent  profile, 
produce  more  stable  jets  than  does  a  fully  developed  laminar  profile.  The  results  of 
Sterling  and  Sleicher  (1975),  together  with  those  of  Tseng  et  al  (1992),  support  this 
conclusion  for  higher  velocity  jets  as  well.  The  physical  reason  for  this  is  that  the  greater 
the  variation  in  the  velocity  profile,  the  greater  the  relaxation  the  jet  undergoes  once  it 
exits  the  nozzle  and  is  no  longer  constrained  by  the  no-slip  condition  at  the  nozzle  wall. 
The  redistribution  of  the  momentum  tends  to  promote  the  growth  of  disturbances  that 
ultimately  break  up  the  jet. 

In  summary,  the  surface  disintegration  of  liquid  jets  (which  may  be  analogous  to 
the  surface  behavior  of  liquid  wall  jets)  is  primarily  a  function  of  the  turbulence  within 
the  jet.  At  low  enough  flow  rates,  the  radial  turbulent  component  is  insufficient  to 
produce  the  drop  forming  filaments.  As  the  velocity  increases  into  the  Wind  Induced 
regimes,  the  radial  turbulent  fluctuations  become  significant  enough  to  lead  to  drop 
production  along  the  jet  surface.  The  source  of  turbulence  can  be  provided  by  pre¬ 
existing  turbulence  upstream  of  the  nozzle,  velocity  profile  relaxation  at  the  discharge  of 
the  nozzle,  and/or  turbulent  boundary  development  along  the  inner  wall  of  the  nozzle. 
Further  discussion  of  the  present  state  of  knowledge  regarding  spray  formation  from  a 
liquid  jet  may  be  found  in  the  review  of  Lin  and  Reitz  (1998). 

Although  not  germane  to  this  research,  it  must  be  mentioned  for  the  sake  of 
completeness  that  there  is  a  large  body  of  knowledge  related  to  the  disintegration  of  thin 
liquid  sheets.  Although  the  mechanics  that  lead  the  sheets  to  deteriorate  into  filaments  of 
fluid  is  much  different  from  that  presented  above,  the  breakup  of  the  filaments  into  drops 
is  similar  to  that  of  jets  in  the  Rayleigh  regime.  For  a  fairly  complete  discussion  of  the 
factors  influencing  the  breakup  of  thin  liquid  sheets,  see  the  works  of  Squire  (1953), 
Dombrowski  and  Fraser  (1954),  Brown  (1961),  Dombrowski  and  Hooper  (1962),  Briffa 
and  Dombrowski  (1966),  Ford  and  Furmidge  (1967),  Crapper  et  al  (1973,  1975), 
Mansour  and  Chigier  (1990, 1991)  and  Chubb  et  al  (1993). 
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2.  Turbulent  Boundary  Layer  and  Open  Channel  Flows 

As  mentioned  earlier,  there  has  been  little  investigation  into  the  nature  of  the 
subsurface  characteristics  in  high-speed  thin  liquid  wall  jets.  Experimentally,  it  is  a 
considerable  challenge  to  measure  the  flow  field  in  the  vicinity  of  a  highly  disturbed  free 
surface  such  as  encountered  here.  Using  a  two-dimensional  laser  doppler  anemometer  in 
an  open  channel  flow  arrangement,  Nezu  and  Rodi  (1986)  could  get  no  closer  than  about 
5  mm  to  a  free  surface  possessing  only  small  variations  in  elevation.  Moreover,  the  work 
which  has  been  performed  in  this  area  is  limited  to  less  challenging  conditions  such  as 
low  Reynolds  or  Froude  number  flows  over  smooth  surfaces.  Hence,  it  is  unclear  how 
the  large  free  surface  structures  arise  from  the  apparently  violent  turbulent  interactions 
beneath  the  surface. 

To  glean  some  insight  into  the  subsurface  characteristics  of  the  jet,  it  is  beneficial 
to  revisit  some  of  the  previous  works  concerning  turbulent  boundary  layer  flows  over 
smooth  and  rough  surfaces.  Because  we  desire  to  contrast  the  smooth-  and  rough-wall 
behaviors,  the  nature  of  these  flows  will  be  discussed  with  their  velocity  components 
normalized  by  either  the  free  stream  or  mean  velocity.  The  reason  for  this  is  that  the  wall 
variables  are  strongly  dependent  on  surface  roughness  as  well  as  the  mean  velocity  profile 
in  the  sublayer.  If  comparisons  are  made  based  on  normalizing  the  velocities  (U,  u',  v’ 
and  w')  by  u*,  then  it  is  more  difficult  to  draw  conclusions  about  the  influence  of  the 
roughness  on  these  parameters,  since  the  normalization  variable  is  also  affected.  Also, 
the  dramatic  nature  of  the  flow  witnessed  at  the  free  surface  of  the  high-speed  wall  jets 
examined  in  this  study  suggests  that  the  mean  behavior  in  the  sublayer  does  not 
adequately  describe  the  local  behavior,  and  the  instantaneous  local  shear  stress  at  the  wall 
will  not  be  measured  directly.  For  all  of  these  reasons,  it  is  all  the  more  logical  to  avoid 
using  the  wall  variables  whenever  possible. 

Figure  2  combines  the  experimentally  determined  mean  velocity  profile  data  of 
Klebanoff  (1954),  Pimenta  et  al  (1975)  and  Ligrani  et  al  (1979)  for  turbulent  boundary 
layer  flow  over  a  smooth-wall,  fully  developed  rough-wall,  and  artificially  thickened  fully 
developed  rough-wall,  respectively.  All  three  curves  are  plotted  for  similar  values  of  Rex. 
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The  "fully  developed  rough-wall"  flow  implies  that  the  data  are  taken  at  a  downstream 
location  where  ks"^  >  60.  The  reduced  values  near  the  wall  for  the  rough-surfaces  are  a 
consequence  of  the  added  drag  of  the  roughness  elements.  One  of  the  more  interesting, 
and  perhaps  most  relevant,  aspects  of  turbulent  boundary  layer  flow  ascertained  from  the 
above  analyses  is  the  change  that  occurs  in  the  streamwise  turbulence  intensity,  Unns,  in 
the  near-wall  region.  As  shown  in  Fig.  3,  the  maximum  value  of  the  rough-wall  Unns  is 
about  10  percent  lower  than  that  of  the  smooth-wall,  and  the  maximum  value  occurs  at  a 
point  further  away  from  the  wall.  Grass  (1971)  was  one  of  the  first  to  notice  this  behavior 
(further  pointing  out  that  as  the  size  of  the  roughness  increases,  Unm  decreases  near  the 
wall).  Based  on  continuity  requirements  alone,  one  would  then  expect  that  the  remaining 
turbulence  intensity  values  must  increase.  As  shown  by  Grass,  and  demonstrated  here  in 
Fig.  4,  this  is  indeed  the  case.  Grass’  results  were  for  low  Reynolds  number  open  channel 
turbulent  flows,  but  the  results  of  Pimenta  et  al  (1975)  and  Ligrani  et  al  (1979)  show  that 
this  behavior  is  also  trae  for  higher  Reynolds  number  boundary  layer  flows.  It  has  been 
understood  for  some  time  that  the  turbulence  produced  by  flow  over  a  smooth-wall  is 
largely  a  result  of  the  bursting  phenomenon.  That  is,  the  lift-off  and  breakdown  of  the 
streamwise  low-speed  streaks  followed  by  an  in-rush  of  fluid  from  the  high-speed  to  the 
low-speed  region.  This  sequence  is  often  defined  in  terms  of  an  ejection  (u'  <  0,  v'  >  0), 
and  sweep  (u’  >  0,  v'  <  0)  process.  The  existence  of  these  streaks  was  first  noted  by 
Klebanoff  and  coworkers  in  the  1950s,  and  their  role  in  turbulence  production  later 
confirmed  by  others  (i.e.,  Kline  et  al  1967).  Kim  et  al  (1971)  claim  that  all  of  the 
turbulence  produced  within  y’^  =  100  is  via  the  ejection  and  sweep  process.  Although  the 
streaks  are  also  present  on  rough  surfaces,  the  higher  level  of  turbulent  kinetic  energy 
(TKE)  (Ligrani  et  al  1979);  the  more  isotropic  nature  of  the  turbulence;  and,  the  structural 
difference  in  the  turbulence  components  (Krogstad  and  Antonia  1994),  suggest  that  the 
predominant  turbulence  producing  mechanism  is  altered  by  the  existence  of  roughness.  It 
is  generally  thought  that,  unlike  the  “gentle”  lift-off  of  the  slow  moving  streaks  on  a 
smooth-wall,  the  ejection  and  sweep  processes  on  a  rough  surface  are  much  more  violent. 
Nonetheless,  for  both  surface  types  Grass  (1971)  notes — and  Raupach  (1981) 


9 


demonstrates  statistically — ^that  the  ejected  fluid  behaves  in  a  quasi-coherent  fashion, 
moving  almost  vertically  through  the  boundary  layer.  Raupach  further  pointed  out  that: 

•  Within  y/5  »  0.2,  nearly  all  of  the  stress  results  from  the  sweep  and  ejection 
process,  with  each  making  equal  contributions.  Furthermore,  they  are  the  only 
source  of  high  level  stress  (i.e.,  from  5  to  10  times  uv ). 

•  This  trend  continues  until  very  near  the  surface  of  a  smooth-wall  (y"^  <  20)  when 
the  sweeps  become  dominant  in  the  viscous  sublayer. 

•  For  rough-walls,  sweeps  are  the  predominant  contributor  to  stress  from  just  above 
the  tops  of  the  roughness  elements  to  down  in  the  roughness  canopy.  The  sweeps 
in  this  region  are  capable  of  producing  stresses  as  high  as  20  times  uv . 

•  As  the  number  of  roughness  elements  is  increased  (with  sufficient  interstitial 
spacing),  the  frequency  of  the  ejections  and  sweeps  also  increases. 

The  ejection  frequency  plot  of  Raupach,  Fig.  5,  supports  the  notion  that  ejections  are 
coherent  motions;  and,  it  shows  that  the  ejection  frequency  measured  at  different  points  in 
the  boundary  layer  remains  essentially  constant.  This  suggests  that  the  ejected  fluid 
proceeds  vertically  with  its  motion  only  minimally  influenced  by  the  local  boundary  layer 
conditions.  In  contrast  to  Raupach's  results,  Krogstad  et  al  (1992)  found  that  the 
frequency  of  the  ejection  and  sweep  events  increases  with  roughness;  and,  that  the 
influence  of  roughness  is  not  confined  to  the  near  wall  region  as  it  is  often  assumed. 
The  reader  is  referred  to  the  review  articles  by  Robinson  (1991)  and  Smith  and  Walker 
(1995)  for  more  detailed  discussions  on  the  nature  of  coherent  structures  in  boundary 
layer  flows. 

Another  relevant  aspect  of  the  boundary  layer  research  which  may  be  applicable  to 
this  study  are  investigations  of  how  the  boundary  layer  responds  to  changes  in  surface 
roughness.  Antonia  and  Luxton  (1971)  showed  that  the  inner  boundary  layer  adjusts  from 
a  fully  developed  smooth  profile  to  a  rough  profile  within  205  of  the  start  of  the 
roughened  region.  In  the  case  of  a  rough-to-smooth  transition,  Andreopoulos  and  Wood 
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(1982)  found  that  the  transition  took  over  555  after  the  flow  was  exposed  to  a  150-mm 
length  of  roughness.  A  similar  result  is  also  noted  by  Pimenta  et  al  (1975)  who  found  that 
a  rough-wall  boundary  layer  transitions  from  laminar  to  turbulent  behavior  fairly  rapidly 
near  the  wall,  but  much  slower  in  the  outer  region  (and  that  the  turbulent  variables 
continue  to  evolve  long  after  the  mean  properties  have  reached  their  fully  turbulent  form). 
As  they  suggest,  the  outer  flow  appears  to  have  a  fairly  long  “memory”  of  its  upstream 
conditions.  Gibbings  and  Al-Shukri  (1997)  conclude  that  roughness  tends  to  shorten  the 
length  of  the  transition  from  laminar  to  turbulent  characteristics,  in  addition  to  moving 
the  beginning  of  this  transition  further  upstream. 

With  the  exception  of  Grass  (1971),  all  of  the  previous  studies  were  limited  to 
boundary  layer  flows.  Therefore,  the  question  of  how  the  free  surface  might  affect  these 
results  needs  to  be  discussed.  As  was  done  by  Finley  et  al  (1966)  for  a  shallow  channel,  it 
is  fairly  common  to  use  the  law  of  the  wall  adjusted  by  a  wake  parameter  to  describe  the 
mean  smooth-wall  open  channel  flow  velocity  profile.  However,  unlike  Finley  et  al,  most 
open  channel  flow  investigations  conclude  that  the  law  of  the  wall  can  be  used  to  describe 
the  entire  region.  For  example,  Nezu  and  Nakagawa  (1993)  have  shown  that 


u^  =  — ln(y^)-HA  +  -^sin^f-^ 
0.41  ^  ^  0.41  2ho 


where  H  is  Coles'  wake  strength  parameter  and  varies  from  zero  at  Re  »  SxlO'^  to  0.2  at 
Re  >  2x10^.  For  rough  surfaces  they  claim  the  same  equation  can  be  used  if  the  middle 


term  is  replaced  with 


-^ln(k^)+A 

0.41  ^  *  ’ 


(5) 


and  y/ks  is  substituted  for  y*  in  the  first  term.  As  the  fully  rough  region  is  reached,  the 
actual  value  of  the  middle  term  departs  from  Eq.  (5),  and  approaches  a  constant  value  of 
8.5.  It  is  thought  that  the  constant  A  can  be  universally  defined  as  5.0  to  5.3,  however 
the  results  of  Tominaga  and  Nezu  (1992)  suggest  otherwise,  especially  at  supercritical 
Froude  numbers.  Likewise,  the  value  of  Coles'  wake  parameter  aqipears  to  vary  rather 


unpredictably  for  changes  in  ks ,  Re  and  Fr. 
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As  was  the  case  in  boundary  layer  flows,  wall  roughness  influences  the  Unns  and 
Vnns  distributions  across  the  channel  depth,  but  not  as  dramatically.  This  is  clearly  seen 
by  comparing  data  from  Nezu  and  Nakagawa  (1993),  shown  in  Fig.  6,  with  that  of  Figs.  3 
and  4.  Figure  6  also  shows  very  little  difference  between  the  smooth  and  rough  flows  in 
the  outer  region  of  the  flow,  but  this  result  is  apparently  limited  to  low  Froude  number  (Fr 
<  1)  flows.  As  Fr  is  increased,  the  turbulent  intensities  in  the  outer  region  grow,  and 
Nezu  and  Nakagawa  attribute  this  to  the  formation  of  waves  on  the  free  surface  that  begin 
contributing  to  the  turbulence  in  the  region. 

Near  the  free  surface,  surface  tension  begins  to  diminish  the  vertical  velocity 
fluctuations  as  described  by  Davies  (1972).  Comparing  open  channel  and  closed  channel 
flows,  Nezu  and  Rodi  (1986)  found  that  the  variations  in  Vnns  were  limited  to  the  region 
above  y/ho  =  0.9,  suggesting  that  the  influence  of  a  calm  free  surface  on  open  channel 
turbulence  is  limited  to  the  upper  ten  percent  of  the  channel.  In  this  region  Komori  et  al 
(1987)  found  that  Unns  increases  at  the  expense  of  Vnns. 

3.  Numerical  Investigations 

As  well  as  the  analytical  efforts  mentioned  earlier,  a  number  of  researchers  have 
used  linear  and  nonlinear  techniques  to  successfully  describe  the  characteristics  of 
laminar  jets  with  finite  disturbances  (for  example,  Lee  1974).  To  effectively  predict 
results  in  liquid  wall  jets  using  a  numerical  scheme,  the  code  must  be  capable  of 
modeling  the  turbulence,  free  surface  deformation  and  surface  tension  over  a  three- 
dimensional  domain.  Moreover,  if  aerodynamic  effects  are  to  be  considered,  it  must  also 
be  able  to  do  the  above  while  simultaneously  solving  the  Navier-Stokes  equations  for  two 
fluids  of  large  density  ratio  (»1000).  The  ability  to  incorporate  a  turbulent  model  with 
that  of  a  reliable  surface  tension  model  is  at  the  forefront  of  CFD  research.  The 
Theoretical  Fluids  Group  at  Los  Alamos  National  Laboratory  has  successfully 
demonstrated  a  two-dimensional  laminar  free  surface  code  (Kothe  and  Mjolsness  1992), 
and  work  continues  on  a  three-dimensional  version  incorporating  turbulence  modeling. 
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At  the  present  time,  however,  it  has  been  ascertained  that  the  existing  commercial  codes 
are  simply  not  capable  of  handing  this  complex  problem. 

It  is  important  to  recognize  that  the  uniqueness  and  challenges  present  in  the 
liquid  wall  jet  make  it  unlikely  that  a  numerical  code  will  successfully  consider  all  of  the 
necessary  parameters  in  a  realistic  manner  a  priori.  The  ability  to  predict  the  behavior 
witnessed  in  this  research  will  likely  depend  on  the  programmer  being  able  to  extract  the 
physics  of  the  process  based  on  insight  provided  by  experiments  such  as  this. 

C.  LIQUID  WALL  JETS 

As  suggested  by  Fig.  7,  a  liquid  wall  jet  is  most  simply  described  as  an  extremely 
shallow  high-velocity  open-channel  flow.  The  liquid  jet  discharges  from  a  nozzle  and  is 
bounded  on  its  lower  side  by  a  wall,  which  can  be  smooth  or  rough,  and  above  by  a 
liquid-air  interface. 


If  the  Reynolds,  Weber  and  Froude  numbers  are  sufficiently  large,  the  boundary  layer 
formed  at  the  wall  quickly  reaches  the  free  surface,  causing  it  to  form  large  three- 
dimensional  disturbances.  As  the  agitation  increases,  some  of  the  disturbances  take  the 
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form  of  long  narrow  protrusions  (also  referred  to  as  fluid  filaments,  ligaments,  columns 
or  jets)  that  can  extend  far  above  the  level  of  the  local  mean  free  surface.  In  some 
instances  these  filaments  break  up  into  one  or  more  drops  as  surface  tension  acts  on 
disturbances  imbodied  within  them  to  reduce  the  overall  energy  of  the  filament.  In  other 
cases,  they  simply  retract  back  into  the  body  of  the  wall  jet  if  the  form  of  the  disturbances 
or  the  length  of  time  that  the  filament  exists  above  the  free  surface  is  insufficient  to  allow 
surface  tension  to  act. 

As  the  jet  flows  downstream,  the  momentum  flux  to  the  wall  causes  the  jet  to 
slow  and  thicken.  The  falling  drops  formed  by  the  breakup  of  the  liquid  filaments,  as 
well  as  the  downward  motion  of  the  larger  three-dimensional  surface  disturbances,  draw 
air  into  the  body  of  the  jet.  This  self-aeration  process  adds  to  the  turbulent  nature  of  the 
flow,  and  reduces  the  mass  density  of  the  jet  as  it  progresses  downstream.  All  of  these 
factors — ^which  appear  to  be  functions  of  the  relative  wall  roughness,  and  the  Reynolds, 
Froude  and  Weber  number  of  the  flow — act  to  prevent  the  flow  from  reaching  any  type  of 
stable  state  or  dynamic  similarity  in  the  regions  investigated  in  this  research. 

D.  CONCLUSIONS 

Based  on  an  extensive  review  of  the  relevant  published  literature,  it  is  clear  that 
little  work  has  been  performed  in  the  area  of  interest  to  this  study:  spray  and  drop 
formation  from  high-speed  liquid  wall  jets.  Nonetheless,  some  important  facts  from  the 
related  research  can  be  transferred  to  the  present  investigation  for  guidance  and 
interpretation. 

Since  the  immediate  concern  is  a  fundamental  understanding  of  the  behavior  of 
bow  sheets,  it  is  important  that  the  liquid  wall  jets  retain  a  character  similar  to  that  of  the 
ocean  surface  encountered  by  a  ship  as  it  glides  through  the  water.  Therefore,  prior  to 
encountering  the  wall,  the  liquid  jet  should  have  a  low  turbulence  intensity  and  no 
boundary  layer  formation.  This  makes  the  nozzle  shape  a  critical  element  of  the 
experimental  design. 
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Researchers  in  the  field  of  free  jets  have  concluded  that  aerodynamics  is  not  a 
significant  factor  vis-a-vis  drop  formation  at  the  surface  of  the  jet.  However,  the 
recognized  self-aeration  process  that  occurs  in  many  channel  flows  suggests  this  may  not 
be  entirely  true.  It  will  be  important  to  determine  whether  or  not  the  relative  motion 
between  the  air  and  the  liquid  wall  jet,  as  well  as  that  of  the  individual  fluid  filaments,  is 
influential  to  the  process. 

The  character  of  the  wall  surface  is  important  because  the  wall  is  the  primary 
source  of  turbulence  in  the  liquid  wall  jet.  The  results  of  both  boundary  layer  and  open 
channel  investigations  show  that  wall  roughness  alters  the  mean  flow  properties  and  their 
moments,  in  addition  to  the  formation  process  of  any  quasi-coherent  structures  existing  in 
the  flow.  Some  studies  have  concluded  that  the  coherent  structures  themselves  are  the 
main  source  of  turbulent  production.  While  it  is  obvious  that  the  turbulence  generated  by 
the  free  surface  disturbances  is  important,  it  is  also  clear  that  these  disturbances  are  a 
product  of  the  boundary  layer  growth;  and,  that  they  alone  cannot  account  for  the 
existence  of  the  liquid  filaments  and  drops.  Therefore,  the  question  of  how  surface 
roughness  and  coherent  turbulent  stractures  influence  the  formation  of  filaments  and 
drops  above  the  jet  remains  unanswered. 

Guided  by  the  purpose  of  the  research,  as  well  of  what  is  currently  understood 
about  liquid  wall  jet  phenomena,  the  goals  of  this  study  are  to: 

•  Better  elucidate  the  characteristics  of  the  free  surface  in  high-speed  liquid  wall 
jets. 

•  To  explain  the  physics  behind  the  behavior  observed  by  qualifying  and 
quantifying  the  influence  of  the  dominant  factors  such  as  relative  roughness, 
Reynolds  number,  Froude  number  and  Weber  number. 

•  To  relate  the  subsurface  nature  of  the  jet  to  that  of  its  free  surface. 

All  of  the  effort  that  has  gone  into  investigating  the  problem  points  out  that  there 
are  a  number  of  issues  affecting  the  drop-formation  process  as  well  as  the  techniques  used 
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to  investigate  it.  These  issues  will  be  addressed  in  the  body  of  the  thesis,  but  it  is 
important  to  point  out  some  of  the  more  significant  ones  now,  enabling  the  reader  to  be 
better  prepared  and  fully  aware  of  some  of  the  limits  of  the  research.  For  example,  it  is 
important  to  understand  that  one  is  dealing  with  a  highly-disturbed  free  surface  that 
makes  the  existing  tools  insufficient  for  measuring  many  of  the  properties  of  interest. 
The  thickness  of  the  liquid  jet  is  so  small  that  it  is  impossible  to  assess  many  of  the 
properties  normally  described  in  boundary  layer  or  channel  flows.  The  large  turbulent 
eddies  that  are  formed  change  their  shape  and  breakup  in  a  rather  random  process,  as  they 
move  rapidly  through  the  jet.  The  phenomenon  of  a  liquid  wall  jet  is  never  stationary. 
From  inception,  the  jet  gives  up  energy  to  the  wall  and  through  the  free  surface,  thereby 
continuously  changing  its  own  character.  This  energy  transfer  is  so  rapid  that  the  jet  loses 
the  ability  to  sustain  these  same  energy  transfer  mechanisms  as  it  approaches  the  end  of 
the  wall.  Reaching  this  latter  stage,  the  changes  in  jet  character  bring  about  a  flow  that  is 
of  no  particular  interest  to  this  investigation.  All  of  the  above  described  events  transpire 
in  less  than  two-tenths  of  a  second,  over  a  distance  of  about  one  meter. 

Again,  the  aim  is  to  understand  the  physics  of  the  spray  formation  process  so  that 
action  may  be  taken  to  reduce  the  danger  this  spray  presents,  primarily  resulting  from  its 
ability  to  scatter  electromagnetic  radiation.  Using  all  means  available,  we  have  arrived  at 
the  results  contained  within.  They  include  the  size,  trends  and  behavior  of  the  drops  and 
filaments  formed  at  the  free  surface  of  liquid  wall  jets.  These  jets  were  developed  over 
smooth  and  rough  surfaces  that  were  either  flat  or  curved. 
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II.  EXPERIMENTAL  EQUIPMENT 


A.  RECTANGULAR  LIQUID  WALL  JET 

The  nozzle  which  forms  the  rectangular  jet  was  fabricated  out  of  two  solid  blocks 
of  aluminum  using  the  wire  electron  discharge  machining  (EDM)  process,  and  its 
surfaces  were  covered  with  a  smooth  polymer  coating  to  prevent  corrosion.  Because  it  is 
recognized  that  the  nozzle  characteristics  have  a  significant  impact  on  the  resulting  jet — 
and  that  these  characteristics  are  not  easily  quantified — every  effort  was  made  to 
incorporate  the  experience  of  previous  researchers  in  this  design  (Dombrowski  and  Fraser 
1954,  Hoyt  and  Taylor  1977a  and  1977b,  McCarthy  and  Malloy  1974,  Rouse  and  Abol- 
Fetouh  1950  and  Rouse  et  al  1951,  Wu  and  Faeth  1995).  The  design  goal  was  a  nozzle 
which  produced  a  rectangular  jet  containing  a  minimum  amount  of  turbulence  and  an 
extremely  thin  boundary  layer,  so  that  the  flow  at  the  exit  would  have  a  nearly  flat 
velocity  profile  (i.e.,  "slug  flow").  By  doing  so  we  ensure  that  the  jet  characteristics  were 
a  result  of  the  fluid  dynamic  influences  occurring  outside  the  nozzle  body. 

The  nozzle  contracts  on  its  upper  and  lower  surfaces,  as  shown  in  Figs.  8  and  9, 
but  the  sidewalls  remain  parallel.  Its  area-contraction  ratio  is  13.9  and  the  shape  of  the 
contraction  is  taken  from  that  which  Rouse  et  al  (1951)  determined  to  be  an  optimum 
shape  for  producing  an  axisymmetric  jet  with  a  maximum  breakup  length.  As  previously 
explained,  a  maximum  breakup  length  is  indicative  of  a  jet  possessing  a  minimum 
amount  of  free-stream  turbulence  and  a  flat  velocity  profile  at  the  nozzle  exit.  The 
entrance  to  the  rectangular  nozzle  was  streamlined  to  ensure  that  the  flow  entering  it  was 
smooth  and  undisturbed,  and  the  discharge  coefficient  for  this  nozzle  system  was  0.91 
over  the  range  of  flow  rates  considered.  As  shown  in  Fig.  10,  the  nozzle  is  mounted  to 
one  side  of  a  23,000  liter  U-tunnel  that  was  previously  used  in  the  oscillating  flow 
research  of  Sarpkaya  (1977,  1986).  The  water  discharged  from  the  nozzle  and  formed  an 
8-mm  thick  by  300-mm  wide  rectangular  free  jet.  After  a  "free  flight"  distance  of  7.6  cm, 
the  main  body  of  the  jet  flowed  over  a  sharp-edged  plate,  thereby  forming  the  liquid  wall 
jet.  The  water  eventually  runs  off  the  end  of  the  plate  (except  for  a  small  amount  which 
flows  over  the  sides),  and  falls  down  into  a  catch  tank  which  is  tilted  away  from  the  U- 
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tunnel  at  a  5°  angle,  allowing  the  lower  end  to  act  as  a  weir.  A  five  horsepower  pump 
takes  suction  on  the  bottom  of  the  tank  and  pumps  the  water  back  into  the  opposite  leg  of 
the  U-tunnel. 

The  hydrostatic  head  in  the  tunnel  determines  the  mean  velocity  of  the  water  at 
the  exit  of  the  nozzle.  During  test  runs,  this  head  is  maintained  constant  by  throttling  the 
ball  valve  at  the  discharge  side  of  the  pump.  The  net  positive  suction  head  (NPSH)  to  the 
pump  is  maintained  constant  by  keeping  the  pump  flow  rate  slightly  less  than  that  of  the 
nozzle.  This  caused  the  catch  tank  to  overflow  a  small  amount  at  its  lower  end,  allowing 
the  weir  to  keep  the  tank  level — and  hence,  the  pump  NPSH — constant  at  the  height  of 
the  lower  end.  Providing  a  small  steady  addition  of  city  water  during  operation 
compensated  for  the  spillage. 

With  this  arrangement,  it  was  possible  for  the  system  to  be  operated  continuously 
for  an  indefinite  period  of  time.  Moreover,  the  large  water  volume  of  the  tunnel  kept  the 
water  velocity  through  the  horizontal  section  of  the  tunnel  below  1.5  cm/sec  even  at  the 
highest  jet  velocity  of  nearly  8  m/s.  This  low  internal  flow  rate  resulted  in  long  fluid 
recirculation  times  which,  along  with  a  six-inch  thick  foam  divider  plate  placed  across  the 
horizontal  mid-section,  ensured  that  the  water  entering  the  nozzle  was  quiescent.  This 
fact  was  verified  by  performing  flow  visualization  with  food  dye  inside  the  U-tunnel. 

The  flow  through  the  nozzle  is  controlled  via  a  pneumatically  operated  cylindrical 
gate.  When  the  gate  is  lifted  as  shown  in  Fig.  1 1,  the  water  discharges  from  the  nozzle, 
and  after  the  short  free  flight,  flows  over  a  plate  with  a  sharp  leading  edge.  In  Fig.  12  the 
gate  is  closed  over  the  nozzle  opening.  The  plate  mount  is  designed  to  accommodate 
plates  of  various  geometries  (i.e.,  flat,  curved  or  varied  geometry)  and  is  capable  of 
inclining  a  given  plate  up  to  15°  from  the  horizontal.  Figure  7  provides  a  side  view  of  the 
jet  discharging  over  a  simple  flat  plate  and  demonstrates  how  the  sharp  leading  edge  is 
positioned  to  cut  off  a  thin  portion  (about  2  mm)  of  the  bottom  side  of  the  liquid  jet.  This 
further  ensures  that  the  flow  engaging  the  upper  surface  of  the  plate  is  free  from  any  free 
shear  layer  effects  that  may  have  developed  upon  exiting  the  nozzle. 

It  was  discovered  by  trial-and-error  that  using  a  plate  that  is  slightly  wider  than 
the  nozzle  minimizes  the  side  spray  produced  when  the  surface-edges  of  the  jet  (distorted 
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by  surface  tension)  encounter  the  comers  of  the  plate.  To  permit  flow  visualization  near 
the  wall,  five  0.8  mm  diameter  dye  ports  were  drilled  about  5  cm  from  the  leading  edges 
of  the  smooth  and  roughest  flat  plates  used  in  this  experiment.  All  of  the  flat  plates  were 
fabricated  from  2.0  or  2.5  cm  thick  acrylic  sheets.  The  concave  curved  plate,  Fig.  13, 
was  fabricated  from  a  3  mm  thick  acrylic  sheet  bent  to  a  3  m  radius  of  curvature.  The 
first  13  cm  of  the  curved  plate  is  flat  and  bonded  to  the  top  of  a  2.5  cm  thick  acrylic  piece 
to  prevent  vibration.  A  0.8  mm  wide  dye  injection  slot  (angled  back  at  45°)  was 
machined  into  this  thicker  region  about  5  cm  from  the  leading  edge.  The  curve  of  the 
plate  begins  13  cm  back  from  the  leading  edge,  and  is  held  rigid  by  vertical  supports  that 
fasten  to  the  plate  mount  on  their  lower  edges.  All  of  the  plates  were  cut  and  the  edges 
machined,  resulting  in  a  final  width  of  33  cm  and  a  length  of  either  120  or  180  cm.  The 
leading  edge  of  each  plate  was  milled  to  produce  a  30°  leading  edge  with  the  sharp-edge 
having  a  radius  of  curvature  of  not  more  than  0.4  mm.  This  ensured  that  the  edge  was 
sharp  enough  to  prevent  the  flow  from  separating,  yet  thick  enough  to  prevent  it  from 
being  easily  damaged. 

As  shown  in  Table  1,  three  different  levels  of  plate  (or  wall)  roughness  were  used 
in  this  investigation.  The  roughness  elements  consist  of  sand  which  was  sifted  numerous 
times  to  achieve  the  desired  size  range.  The  roughness  is  applied  to  the  plate  by  coating 
the  top  with  an  extremely  thin  layer  of  epoxy  (applied  with  a  brush),  and  then  gently 
shaking  a  sand-filled  sifter  over  the  surface.  This  resulted  in  walls  with  evenly  coated 
rough  surfaces. 


Roughness 

E>escription 

Minimum 
element  size  (mm) 

Maximum 
element  size  (mm) 

Mean 

element  size  (nun) 

Fine 

0.104 

0.147 

0.13 

Medium 

0.300 

0.355 

0.33 

Coarse 

0.590 

0.840 

0.72 

Table  1.  Plate  roughness  levels. 


The  walls  were  mounted  in  the  experimental  facility  by  bolting  them  from 
underneath  to  a  mount  fabricated  out  of  6.4  mm  aluminum.  As  shown  in  Fig.  14,  the 
front  and  rear  of  the  mount  rests  on  two  horizontal  rods,  and  the  ends  of  these  rods  are 
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individually  pinned  to  fine-control  traversing  mechanisms.  This  allows  the  attached  wall 
to  be  placed  at  angles  of  attack  from  -5°  to  15°  ±0.25°,  and  provides  a  small  amount  of 
lateral  angle  adjustment.  The  slots  in  the  top  of  the  mount  permit  the  wall  to  be  placed  at 
various  axial  positions,  however,  for  all  of  the  experiments  presented  here  the  plate  was 
positioned  to  give  the  jet  a  free-flight  distance  of  7.6  cm. 

The  thickness  of  the  jet  was  measured  along  the  centerline  of  the  plate  at  various 
streamwise  positions  using  a  scale  inserted  from  above.  The  scale  was  fabricated  from  a 
0.23-mm  thick,  4-mm  long,  stainless  steel  razor  blade.  The  device  was  placed  in  the  jet 
with  its  length  parallel  to  the  direction  of  flow,  and  adjusted  so  that  it  produced  no 
noticeable  wake.  The  jet  thickness  was  then  determined  by  reading  the  division  markings 
at  the  front  of  the  scale  (from  video  recordings).  An  attempt  was  made  to  use  a 
capacitance  probe,  but  it  proved  to  be  highly  non-linear  and  the  results  were  not 
repeatable.  The  reasons  for  this  will  be  explained  later. 

B.  AXISYMMETRIC  LIQUID  WALL  JET 

An  axisymmetric  jet  facility  was  constructed  to  overcome  some  of  the 
experimental  and  visualization  limitations  experienced  with  the  rectangular  jet  and  to 
allow  for  different  geometry  wall  jets.  The  nozzle.  Fig.  15,  is  mounted  in  the  horizontal 
section  of  the  U-tunnel,  downstream  of  the  foam  divider  plate.  The  entrance  region  is 
streamlined  using  a  half-toroidal  (or  "donut")  shape  that  mounts  against  the  wall  of  the 
tank  and  possesses  an  inner  diameter  identical  to  that  of  the  nozzle  exit.  The  overall 
contraction  ratio  is  17.6  and  the  discharge  coefficient  of  the  system  was  experimentally 
determined  to  be  0.77.  The  flow  rate  through  this  nozzle  is  substantially  smaller  than  the 
rectangular  nozzle,  so  it  is  not  necessary  to  provide  it  with  a  recirculation  system. 
Instead,  the  nozzle  discharges  into  a  catch  basin  that  empties  to  the  laboratory  drainage 
system.  The  head  in  the  U-tunnel  is  maintained  constant  by  throttling  on  the  city  supply 
line. 

It  is  envisioned  that  a  variety  of  shapes  such  as  cones,  rods  or  wedges  could  be 
placed  at  the  outlet  of  the  nozzle,  but  for  this  work,  only  a  rod  is  used.  Figure  16 
provides  a  side  view  of  the  apparatus.  A  thick-walled  tube  is  mounted  on  the  end  of  the 
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acrylic  test  rod,  and  a  plate  with  a  foam  face  is  mounted  near  the  point  where  the  tube  and 
the  rod  are  joined.  When  the  assembly  is  fully  out,  the  foam-covered  plate  acts  as  a 
stopper — with  the  assistance  of  hydrostatic  pressure — ^to  block  flow  through  the  nozzle. 
When  the  assembly  is  fully  inserted  against  the  back  wall  of  the  U-tunnel,  flow  passes  out 
of  the  nozzle  and  along  the  test  rod.  The  pole  over  which  the  assembly  moves,  and  the 
linear  bearing  that  supports  the  far  end  of  the  test  rod,  are  carefully  positioned  to  ensure 
the  device  is  centered  in  the  nozzle.  This  guarantees  that  the  jet  thickness  is  uniform 
about  the  rod  when  it  exits  the  nozzle.  The  significance  of  this  arrangement  is  that  there 
is  no  stagnation  point  where  the  flow  encounters  the  rod.  Thus,  no  undesirable  spray  or 
other  disturbances  are  produced;  and,  it  provides  a  clearer  visual  presentation  of  the 
events  that  occur  since  one  is  not  viewing  the  region  of  interest  through  many  centimeters 
of  other  flow  disturbances.  The  limitation  of  the  system  is  that  there  is  only  a  short 
region  over  which  the  flow  may  be  examined  before  gravity  begins  to  noticeably  alter  the 
jet  thickness  about  the  rod. 

A  single  test  rod  was  roughened  (using  the  medium  roughness  sand)  as  described 
for  the  flat  plates  except  that  the  rod  was  slowly  rotated  on  a  lathe  during  the  application 
process.  Additionally,  two  1  mm  wide  by  10  cm  long  diametrically  opposed  strips  were 
clear  of  roughness  to  permit  a  laser  light  sheet  to  pass  through  the  rod  for  flow 
visualization.  As  shown  in  Fig  15,  the  nozzle  exit  diameter  is  3.18  cm,  while  the  rod  in 
Fig.  16  has  a  diameter  of  2.54  mm,  giving  a  jet  thickness  of  3.18  mm  (about  half  that  of 
the  rectangular  nozzle). 

Unlike  the  rectangular  nozzle  that  produces  a  jet  with  a  fixed  initial  thickness,  ho, 
the  axisymmetric  arrangement  allows  ho  to  be  varied  by  simply  installing  a  different 
diameter  rod  (a  quick  and  easy  process).  This  provides  a  greater  variety  of  jet  Reynolds, 
Weber  and  Froude  number  combinations  which  can  be  examined.  However,  before 
substituting  the  axisymmetric  results  for  those  of  the  rectangular  device,  it  must  first  be 
verified  that  ho/Drod,  where  Drod  is  the  diameter  of  the  test  rod,  is  small  enough  to  prevent 
the  region  of  interest  from  being  affected  by  the  curvature  of  the  rod. 
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c. 


RECRICULATING  FREE-SURFACE  WATER  TUNNEL 


Both  of  the  above  nozzle  systems  provide  the  capability  to  visualize  and  quantify 
the  surface  features  of  liquid  wall  jets.  However,  the  narrow  jet  thickness  coupled  with 
the  large  free  surface  distortions  made  it  impossible  to  measure  the  subsurface 
characteristics  with  any  degree  of  certainty.  Clearly,  one  could  do  so  with  a  somewhat 
thicker  jet,  provided  that  the  sides  of  the  jet  are  optically  transparent  walls  that  will 
permit  the  use  of  a  digital  particle  image  velocimeter  laser  (DPIV).  A  low-turbulence 
water  tunnel  used  in  previous  works  (Sarpkaya  and  Neubert  1994)  provided  the  most 
direct  means  of  achieving  this  objective.  The  tunnel  has  a  test  section  with  thick 
Plexiglas  bottom  and  sidewalls.  A  35  horsepower  pump  provides  flow  through  a  38  x  70 
X  125  cm  test  section  via  a  flow  straightening  honeycomb  mesh  and  a  settling  chamber. 
For  this  experiment,  a  2.5  cm  thick  clear  acrylic  plate  is  mounted  to  the  bottom.  The 
plate  extended  30  cm  into  the  settling  chamber  and  had  a  semi-circular  leading  edge.  The 
sides  of  the  plate  are  machined  to  exactly  match  the  width  of  the  test  section,  and  the 
desired  jet  thickness  is  obtained  by  clamping  a  vertical  aluminum  gate  to  the  wall  at  the 
entrance  of  the  test  section.  The  gate  is  fabricated  from  a  6.4  mm  thick  aluminum  sheet, 
with  the  upstream  bottom  edge  machined  to  a  3.2  mm  radius  of  curvature,  and  the  trailing 
edge  into  a  sharp  90°  comer.  The  gate  can  be  raised  or  lowered  to  achieve  the  desired  jet 
thickness. 

The  entrance  region  to  the  water  tunnel  jet  was  significantly  different  from  that  of 
the  rectangular  jet  since  it  was  not  formed  using  an  optimum  nozzle  shape,  nor  did  it  have 
an  interval  of  free  flight  before  encountering  the  wall.  The  wall  itself  acted  as  the  lower 
half  of  the  nozzle  and  the  contraction  of  the  settling  chamber  provided  a  small  amount  of 
fluid  acceleration.  The  majority  of  the  area-contraction  was  achieved  by  the  vertical  plate 
at  the  entrance  to  the  test  section.  Overall,  the  fluid  in  the  water  tunnel  experienced  an 
effective  contraction  ratio  of  about  55  before  reaching  the  test  section.  The  recirculating 
pump  of  the  water  tunnel  had  a  variable  speed  controller  capable  of  providing  flows  with 
mean  velocities  from  0.5  m/s  to  10  m/s  for  the  wall  jet  examined  in  this  investigation. 

The  surface  of  the  flat  plate  is  roughened  by  applying  sand  to  the  surface  of  a  3 
mm  thick  acrylic  sheet,  and  then  attaching  this  sheet  to  the  flat  plate  installed  in  the  test 
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section.  A  1.0  mm  wide  strip  along  the  center  region  of  the  roughened  plate  was  clear  of 
sand  to  allow  a  laser  light  sheet  to  penetrate  into  the  flow  region.  A  rough  plate  made 
from  the  coarse  sand,  and  a  smooth  plate  were  examined  in  this  research. 

All  the  analyses  performed  in  the  water  tunnel  were  constrained  to  the  center  of 
the  test  section  to  eliminate  the  influence  of  the  sidewalls  on  the  observed  flow. 


D.  MEASUREMENT  EQUIPMENT 

Two  primary  means  were  used  to  gather  the  data  in  this  research.  The  first  is  that 
of  high-speed  imagers  and  cameras.  The  second  is  digital  particle  image  velocimetry 
(DPIV)  which  was  used  to  obtain  the  subsurface  data  in  the  water  tunnel.  The  constantly 
varying  optical  path  and  narrowness  of  the  region  made  the  use  of  DPIV  and  laser 
doppler  velocimetry  (LDV)  impractical  for  the  rectangular  and  axisymmetric  jets. 


1.  Imaging 

As  listed  in  Table  2,  various  coupled  capacitive  discharge  (CCD)  imagers  (also 
referred  to  as  cameras)  were  used  to  collect  data  for  flow  visualization  as  well  as  images 
that  could  later  be  quantitatively  analyzed. 


Imager 

Speed 

(fi'ames/sec) 

Recording 

Time 

Shutter  Speed 
(sec) 

Recording 

Media 

NAC  HSV-500 

250  or  500 

14  min 

1/250  - 1/10,000 

Redlake  HR  2000 

50  -  8,000 

4-8  sec 

1/50  - 1/40,000 

Camera  RAM 

NEC  TI-22P 

30 

120  min 

1/30 

Table  2.  Description  of  imagers. 


The  data  used  for  quantitative  analysis  of  the  surface  of  the  rectangular  liquid  wall 
jet  was  gathered  using  the  NAC  camera.  It  was  mounted  on  a  tripod  along  side  of  the 
wall  jet  at  the  desired  location  downstream  of  the  leading  edge  of  the  plate  (x/ho  =  0). 
During  the  image  capturing  process,  a  scale  was  placed  in  the  focal  plane  and  recorded 
for  later  dimensional  scaling.  A  typical  example  of  the  jet  free  surface  image  is  shown  in 
Fig.  17  (in  this  and  all  subsequent  figures,  flow  is  from  right-to-left).  After  an  image  or  a 


23 


sequence  of  images  was  captured  on  videotape,  it  was  downloaded  to  a  PC  using  a  frame 
grabber  (Imagraph  Imascan  Prism/Spectram  Board  Set),  and  subsequently  analyzed  using 
the  Optimas  6.1  and  Optimas  MA  3.1  software.  This  software  permits  the  user  to  define 
points,  lines  and  areas  in  an  image,  and  to  track  these  objects  through  a  sequence  of 
images  to  determine  their  positions,  velocities,  accelerations,  etc.  Once  all  of  the  values 
of  interest  have  been  identified,  the  Optimas  software  then  provides  the  output  in  a  text 
format  or  directly  loads  it  into  a  Microsoft  Excel  spreadsheet.  Virtually  all  of  the  images 
obtained  in  this  manner  were  captured  using  a  medium  intensity  diffused  backlighting, 
and  a  high  quality  camera  lens.  The  lighting  was  consistently  set  so  as  to  maintain  a 
minimum  depth  of  field,  ensuring  that  only  features  appearing  in  the  central  region  of  the 
plate  were  in  clear  focus.  In  all  instances  the  total  depth-of-field  was  less  than  8  cm  with 
a  horizontal  field-of-view  of  about  12  cm. 

A  time  history  of  a  given  set  of  surface  stmctures  was  acquired  by  moving  the 
camera  as  near  as  possible  to  the  velocity  of  the  flow.  This  was  also  used  to  capture  the 
transition  of  the  surface,  as  well  as  the  development  of  filaments  and  the  formation  and 
demise  of  drops  in  the  streamwise  direction. 

Additional  flow  visualization  was  performed  by  injecting  food  coloring  through 
small  dye  ports,  or  by  positioning  a  laser  light  sheet  at  various  angles  relative  to  the 
streamwise  direction.  As  discussed  later,  these  techniques  aided  in  verifying  that  the 
flow  in  the  center  of  the  liquid  wall  jet  formed  over  flat  and  curved  walls  is  unaffected  by 
the  side  conditions.  Therefore,  based  on  observations,  the  subsequent  data  analysis  of  the 
rectangular  jet  was  limited  to  the  center  of  the  plate,  so  the  results  can  confidently  be  said 
to  be  those  exhibited  by  a  "two-dimensional"  jet. 

2.  Digital  Particle  Image  Velocimetry  (DPFV) 

As  previously  discussed,  the  DPIV  system  was  used  exclusively  in  the  low- 
turbulence  water  tunnel  to  determine  the  jet’s  subsurface  characteristics  (U,  v',  u',  Unns 
and  Vnns)  in  a  plane  defined  by  the  laser  light  sheet.  The  TSI  DPIV  system  is  comprised 
of  a  PC  containing  a  high-speed  digital  frame  grabber  and  the  TSI  Insight  1.3  software 
for  system  control.  A  separate  synchronizer  coordinated  the  pulses  of  the  dual  YAG 
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lasers  with  that  of  a  Kodak  MegaPIus  digital  imager  possessing  a  1008x1016  pixel 
resolution.  The  software  was  capable  of  capturing  up  to  1,000  sequential  images.  The 
laser  pulse  separation  was  held  at  0.01  ms  while  operating  the  system  at  its  highest  image 
capture  rate  of  30  Hz. 

The  system  was  always  used  in  the  two-frame  crosscorrelation  frame-straddling 
mode.  The  advantage  is  that  the  frame-straddling  mode  provides  the  greatest  system 
bandwidth  (Keane  and  Adrian,  1992).  The  disadvantage  of  this  method  is  that  two 
separate  images  are  needed  for  each  flow  field  calculation.  However,  even  if  the  system 
were  operated  in  a  different  mode  which  made  use  of  its  highest  capture  rate,  the 
maximum  frequency  that  could  be  determined  based  on  the  Nyquist  criterion  was  15  Hz. 
For  the  anticipated  time  scales  present  in  this  flow,  this  data  rate  is  still  too  low  to  obtain 
meaningful  time  history  information  or  perform  useful  spectral  analysis.  Therefore,  the 
sequentially  obtained  PIV  data  must  be  analyzed  as  an  ensemble  of  samples  regardless  of 
the  operating  mode,  so  using  the  system  in  frame-straddling  presented  no  actual  system 
limitation. 

The  flow  was  seeded  with  8  jxm  glass  beads  having  a  specific  gravity  of  about 
1.08,  and  the  back  wall  of  the  test  section  below  the  mean  free  surface  was  blackened  to 
enhance  particle  identification.  The  area  above  the  surface  was  well  lit  to  prevent  the 
camera  from  also  capturing  dust  particles  moving  in  the  air  (which  made  it  more  difficult 
to  distinguish  between  motions  above  and  below  the  free  surface).  Regardless,  the  free 
surface  disturbances  and  occasional  free  surface  reflections  of  the  laser  light  made  it 
impossible  to  acquire  data  close  to  the  fi-ee  surface. 

During  post-processing  of  the  images,  the  Whittaker  peak  search  algorithm  was 
used  to  analyze  the  images  since  it  is  accurate  and  requires  a  minimum  amount  of 
computational  effort  (Lourenco  and  Krothapalli  1995).  The  seeding  was  sufficient  to 
allow  the  minimum  spot  size  of  16  pixels  to  be  used  with  fewer  than  about  5%  of  the 
velocity  vectors  categorized  as  invalid.  The  small  spot  size  (which  has  a  diameter  of 
about  0.5  (im),  combined  with  a  50%  spot  overlap,  resulted  in  fairly  smooth  vorticity  and 
instantaneous  velocity  contour  plots,  and  provided  the  possibility  that  fluid  behaving  in  a 


25 


quasi-coherent  fashion  might  be  identifiable.  The  total  area  examined  during  a  PIV  mn 
was  about  500  mm^  (15  x  33  mm). 


26 


in.  EXPERIMENTAL  RESULTS 


A.  PRELIMINARY  ANALYSIS 

The  overall  characteristics  of  filaments  and  drops  were  examined  prior  to  the  data 
gathering  process.  The  examination  was  carried  out  by  making  high-speed  recordings  of 
the  liquid  wall  jet  from  the  side  and  directly  above.  A  stroboscope  directed  down  toward 
the  surface  of  the  jet  confirmed  that  the  early  stages  of  the  surface  were  laminar,  but  it 
was  also  noted  that  the  surface  contained  many  small  streamwise  capillary  streaks.  Hoyt 
and  Taylor  (1979)  noticed  a  similar  occurrence  in  their  round  nozzle  experiments,  which 
they  attributed  to  minute  machining  grooves  that  existed  on  the  inner  wall  of  the  nozzle. 
Nonetheless,  because  this  apparatus  is  significantly  different  from  that  of  Hoyt  and 
Taylor,  it  was  important  that  steps  be  taken  to  remove  the  streaks,  or  at  least  to  identify 
their  cause.  To  start  with,  the  water  in  the  U-tunnel  was  allowed  to  set  for  many  days 
before  use  so  that  it  could  self  de-aerate.  A  foam-covered  splash  plate  was  also  installed 
just  below  the  pump  discharge  to  reduce  the  amount  of  air  entrained  as  the  recirculated 
water  dropped  back  into  the  tunnel  (Fig.  10).  A  small  viewing  window  and  a  movable 
dye-injection  rig  were  installed  on  the  tunnel  wall  directly  behind  the  nozzle.  These  were 
used  to  perform  flow  visualization  around  the  nozzle  entrance  region  inside  the  U-tunnel, 
and  to  verify  that  the  inner  surface  of  the  nozzle  was  clean  and  free  of  air  bubbles.  Based 
on  the  results  of  these  actions,  it  was  concluded  that  the  small  (<  0.25  mm)  streamwise 
streaky  structures  must  be  the  consequence  of  minute  irregularities  on  the  inner  nozzle 
surface.  Moreover,  their  existence  was  also  seen  to  be  a  function  of  Uo,  the  jet  velocity 
at  the  nozzle  exit.  As  Uo  was  reduced,  the  streaks  eventually  disappeared,  leaving  a 
glass-like  free  surface.  The  flow  visualization,  also  confirmed  that  large  air  bubbles  or 
debris  on  the  wall  of  the  nozzle  were  capable  of  producing  large  (V2  to  2  mm)  capillary 
waves  that  danced  back  and  forth  across  the  free  surface.  Thus,  for  each  mn,  the  nozzle 
surface  was  cleaned  and  polished  before  filling  the  tunnel. 

Laser-induced  fluorescence  (LIF)  was  used  to  investigate  the  free  surface  of  the 
jet  by  directing  a  laser  light  sheet  downward  onto  the  surface,  and  then  positioning  it  at 
various  angles  to  the  flow.  This  examination  made  it  clear  that  the  free  surface 
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characteristics  were  uninfluenced  by  the  edges  of  the  jet  and  wall  over  a  central 
transverse  span  in  excess  of  20  centimeters.  It  further  indicated  that  most  of  the  well- 
defined  filaments  were  axisymmetric  in  nature  and  possessed  no  significant  lateral 
motion:  they  grew  and  receded  in  the  (x,y)  plane  as  they  traversed  with  the  flow  in  the  x- 
direction.  This  is  not  to  say  that  three-dimensional  surface  disturbances  do  not  exist  (Fig. 
18),  but  only  that  most  drop-forming  filaments  were  similar  to  that  of  Fig.  19. 

During  the  early  stages  of  the  investigation,  dye  lines  were  connected  to  the  dye 
ports  that  spanned  the  center  section  of  the  smooth  and  the  roughest  wall.  The  flow 
through  the  lines  was  sufficiently  restricted  so  that  the  suction  provided  by  the  jet  gently 
drew  the  dye  (food  coloring)  near  the  plate  surface.  Hence,  the  movement  of  dye  into  the 
body  of  the  jet  had  little  impact  on  the  flow  field.  This  examination  revealed  that  the 
flow  at  the  wall  became  turbulent  very  near  the  dye  ports  regardless  of  the  level  of  wall 
surface  roughness;  and,  that  the  near-wall  flow  in  the  center  region  was  unaffected  by  the 
wall  edges.  The  same  results  were  witnessed  in  the  curved  wall  flow  also. 

For  all  of  the  rans  except  one,  the  plate  was  positioned  parallel  to  the  direction  of 
flow,  with  the  height  set  so  that  the  sharp  leading  edge  trimmed  off  a  small  amount  of  the 
bottom  side  of  the  liquid  jet  (Fig.  7).  This  ensured  that  the  flow  encountering  the  wall 
was  absent  of  nozzle  shear  layer  effects,  and  of  the  previously  described  capillary  waves. 
This  resulted  in  an  experimental  apparatus  that  was  a  simplified  version  of  a  bow  sheet 
produced  by  a  ship  hull.  The  trimming-off  of  the  lower  side  of  the  jet,  combined  with  a 
small  amount  of  contraction  at  the  nozzle  exit,  resulted  in  a  jet  with  a  mean  initial 
thickness,  ho,  of  about  5.7  mm.  In  the  following  discussion,  the  results  of  the  same  jet 
over  different  surfaces  are  compared  and  contrasted. 

It  was  not  known  a  priori  which  of  the  parameters  (Re,  We,  Fr  or  Z)  governed  and 
which  influenced  the  creation  of  filaments  and  drops  at  and  above  the  free  surface. 
Perhaps  not  entirely  unexpected,  the  smooth  plate  proved  to  be  very  resilient  to  drop 
formation:  although  the  roughest  plate  displayed  drops  for  velocities  as  low  as  3  m/s,  the 
smooth  plate  did  not  do  so  until  the  velocity  approached  7  m/s.  At  this  high  of  a  jet 
velocity,  the  free  surface  for  the  rough-wall  case  was  too  violent  to  allow  the  evaluation 
of  individual  filaments  and  drops.  It  is  primarily  because  of  this  reason  that  the  filament 
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and  drop  analysis  were  performed  at  Uo  =  6.2  m/s  for  all  rough-walls,  and  at  Uo  =  7.7 
m/s  for  the  smooth-plate.  This  higher  smooth-wall  velocity  was  the  minimum  necessary 
to  obtain  a  statistically  significant  sample  of  filaments  and  drops.  The  roughest  wall  was 
also  analyzed  at  a  jet  velocity  of  Uo  =  4.5  m/s  to  examine  the  effect  of  lower  speed  jets 
over  rough  surfaces.  The  effect  of  the  surface  curvature — an  integral  element  of  the  ship 
bow-sheet  phenomena — was  investigated  through  the  use  of  a  rough-plate  (k/ho  =  0.06) 
having  a  3  meter  concave  radius-of-curvature  (Fig.  13).  This  amount  of  curvature,  at  this 
level  of  roughness,  is  similar  to  that  found  on  the  bows  of  large  ships  (Schlichting,  1978). 
The  Ohnesorge  number,  Z,  remained  constant  at  0.0015  because  all  experiments  were 
performed  using  fresh  water  and  a  fixed  jet  thickness. 

In  addition  to  the  filament  and  drop  quantification  analysis,  the  free  surface 
transition  regions  were  visually  identified.  This  analysis  was  performed  for  all  of  the 
above  scenarios,  as  well  as  the  smooth-wall  at  Uo  =  6.2  m/s  and  a  medium-roughness 
wall  inclined  at  a  15°  angle.  Table  3  provides  a  summary  of  all  the  rectangular  jet 
experiments  and  their  relevant  parameters. 
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Table  3.  Summary  of  the  rectangular  liquid  wall  jet  parameters  and  analyses. 


The  Reynolds,  Weber  and  Froude  numbers  could  not  be  varied  independently,  so 
knowing  any  one  of  them  implied  knowledge  of  the  other  two.  Therefore,  specific  test 
runs  are  frequently  described  by  their  relative  roughness  and  Weber  number.  For 
example,  "0.06/3000"  stands  for  the  wall  jet  run  over  the  plate  with  a  relative  roughness 
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of  k/ho  =  0.06  and  a  Weber  number  of  3000.  It  is  further  understood  that  the  terms  jet 
"speed",  "velocity"  and  "inertia"  will  be  used  interchangeably. 

B.  DATA  ACQUISITION  PROCESS  AND  MEASUREMENT  DEFINITIONS 

hi  the  course  of  the  presentation  and  discussion  of  the  results,  the  most  important 
outward  manifestation  of  the  free  surface  (i.e.,  filament  or  ligament)  will  be  referred  to 
quite  frequently.  It  is  for  this  reason  that  we  will  define  a  "filament"  as  precisely  as 
possible  within  the  statistical  bounds  of  its  finite  life.  To  be  sure,  the  definition  is  a 
subjective  matter.  Qualitatively,  a  filament  should  be  long  and  thin,  and  its  two- 
dimensional  image  must  readily  lend  itself  to  a  mental  translation  of  a  three-dimensional 
axisymmetric  object.  Quantitatively,  it  was  required  that  the  slenderness  ratio  (Lf/Df) 
reach  a  value  of  about  1.5  to  2  prior  to  forming  a  drop.  The  only  exception  being  nascent 
filaments  which  were  considered  acceptable  if  they  behaved  axisymmetrically  and 
reached  a  value  of  Lf/Df «  2  before  drop  formation.  In  most  cases,  the  selected  filaments 
were  already  "mature,"  with  slenderness  ratios  greater  than  three,  by  the  time  they  were 
under  consideration  for  analysis.  Based  on  a  summary  of  general  observations,  one  could 
be  quite  certain  that  if  a  filament  possessed  a  three-dimensional  character  (other  than 
small  axisymmetric  or  longitudinal  irregularities),  it  would  make  itself  known  within  two 
or  three  sequential  images.  Filaments  meeting  the  above  criteria  also  had  to  be 
measurable  over  a  period  of  time.  Hence,  their  bases  and  profiles  had  to  be  clear  of  other 
disturbances  that  might  be  present.  In  the  case  of  runs  0.00/4700,  0.02/3000  and 
0.13/1600  (i.e.,  low  k/ho  or  low  We),  most  of  the  filaments  appeared  as  "trees  in  the 
desert."  That  is,  the  surrounding  disturbances  were  typically  much  smaller  than  the 
filament.  In  run  0.13/3000,  however,  the  free  surface  profile  was  more  like  a  forest,  and 
suitable  filaments  often  had  to  be  bypassed  as  their  images  where  blocked  by  other 
filaments  and  large  disturbances.  Run  0.06/3000  was  a  transition  between  these  two 
extremes.  In  rare  instances,  small  filaments  actually  formed  from  large  three- 
dimensional  disturbances  as  shown  in  Fig.  18.  Filaments  such  as  this  were  avoided  since 
it  was  felt  that  there  was  a  greater  risk  of  their  genesis  not  being  directly  linked  to  the 
subsurface  events. 
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The  video  images  used  in  the  quantitative  analysis  of  the  rectangular  liquid  wall 
jet  were  captured  using  the  NAC  HSV-500  imaging  system  operated  at  500  frames/sec 
and  a  shutter  speed  of  5000  sec’’.  The  camera  was  placed  on  the  same  horizontal  plane  as 
the  wall  over  which  the  jet  formed,  and  perpendicular  to  the  streamwise  direction.  Before 
recording,  the  jet  flow  was  started  and  given  sufficient  time  to  stabilize  (5  to  10  min). 
As  mentioned  previously,  the  focal  plane  of  the  camera  was  at  the  center  of  the  wall  and 
had  a  12-cm  wide  field-of-view.  Starting  at  x  =  0  cm,  the  imaging  system  recorded  the 
jet  profile  long  enough  to  provide  sufficient  data  for  analysis,  before  being  repositioned 
12  cm  further  downstream  and  the  recording  process  repeated.  This  sequence  of  events 
was  continued  along  the  length  of  the  jet.  This  particular  camera  position  was  chosen 
because: 

•  It  guaranteed  that  only  filaments  and  drops  along  the  centerline  of  the  jet  were 
in  focus;  hence,  the  results  were  not  influenced  by  the  edges  of  the  jet. 

•  It  was  close  enough  to  ensure  that  the  filament  and  drop  features  could  be 
accurately  defined. 

•  It  was  far  enough  away  to  capture  from  8  to  15  images  of  identifiable 
filaments  and  drop(s)  as  they  passed  across  the  field-of-view  of  the  camera. 
This  was  sufficient  to  ensure  that  the  desired  silhouette  of  the  filament 
exhibited  an  axisymmetric  character,  and  that  the  drops  remained  nearly 
circular. 

Figure  19  provides  an  example  of  a  sequence  of  images  used  to  analyze  a  filament  and 
drop  (the  filament  of  interest  is  indicated  with  an  arrow).  For  this  investigation,  over  800 
filament-drop  pairs  such  as  this  were  analyzed,  which  amounts  to  approximately  8,000 
individual  image  analyses. 

All  of  the  recorded  images  were  obtained  using  the  same  camera,  lens  and 
lighting  arrangement.  After  carefully  reviewing  the  tape,  the  images  of  filaments  and 
drops  whose  motions  were  identifiable  with  little  or  no  ambiguity  were  digitized  and 
downloaded  to  a  PC,  and  the  sequence  of  images  then  analyzed.  The  cross-sectional  area 
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of  the  filament  and/or  drops  were  outlined  with  the  help  of  a  "mouse,"  and  the  base  and 
tip  of  the  filament  were  marked.  The  Optimas  software  calculated  the  cross-sectional 
area,  length,  position,  angle,  velocity  and  acceleration  of  every  feature  identified  in  the 
digitized  image  sequence,  based  on  the  specified  length-scale  and  the  time  between  the 
images  (2  ms).  The  uniqueness  of  each  filament  required,  a  simplified  means  of 
identifying  and  comparing  their  characteristics.  As  explained  earlier,  the  analysis  was 
limited  to  those  filaments  possessing  an  apparent  axisymmetry.  It  was  therefore, 
reasonable  to  identify  their  geometry  with  that  of  a  cylinder  as  shown  in  Fig.  20.  The 
major  uncertainties  associated  with  this  filament  model  were  that  the  base  could  not  be 
precisely  defined,  and  that  the  shape,  as  per  definition,  was  axisymmetric.  However, 
independent  analysis  by  two  separate  observers  indicated  that  the  overall  filament 
uncertainty  is  less  than  10  percent.  The  mean  diameter  of  a  drop  is  determined  from  the 
average  of  its  cross-sectional  area  measurements.  A  complete  discussion  of  the 
uncertainty  involved  with  the  filament  and  drop  analysis  is  presented  in  Appendix  B. 

C.  LIQUID  WALL  JET  ATTRIBUTES 

For  each  of  the  eight  mns  listed  in  Table  3,  the  regions  of  surface  transition  and 
jet  thickness  were  determined  to  help  elucidate  any  general  trends  or  features  that  might 
be  present  in  the  flow.  The  plate  positions  were  carefully  set,  and  the  flow  was  allowed 
to  stabilize  (as  indicated  by  a  constant  water  level  in  the  U-tunnel)  prior  to  beginning  any 
of  the  measurements. 

1.  Transition  Regions 

fii  examining  the  gross  features  of  liquid  wall  jets,  it  was  found  that  the  surface 
characteristics  could  be  divided  into  four  identifiable  regions,  with  examples  of  three  of 
the  four  depicted  in  Fig  21.  The  gradual  nature  of  surface  transitioning  defies  precise 
position  identification,  as  in  all  "transitions"  in  fluid  motion,  nonetheless,  these  locations 
was  measured  as  accurately  as  possible  (±  3  cm)  using  a  number  of  different  camera 
angles,  fi-ame  speeds,  shutter  speeds  and  lighting  arrangements. 
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The  first  transition,  transition  to  roughness  (Fig.  21. a),  is  the  region  in  which  the 
turbulent  boundary  layer  reaches  the  free  surface,  noticeably  roughening  it.  In  this 
region,  three-dimensional  disturbances  appear,  but  filaments  are  not  yet  evident. 
Although  these  disturbances  are  small,  they  do  manage  to  completely  obliterate  the 
streamwise  capillary  waves  in  very  short  order  (2  to  3  cm).  How  visualization  showed 
that  the  turbulent  boundary  layer  began  almost  immediately  regardless  of  the  wall 
roughness,  but  the  distance  necessary  for  it  to  reach  the  free  surface  (about  30ho  to  55ho), 
did  depend  on  k/ho  and  Uo.  Figure  22  shows  that  the  distance  to  the  surface  transition  is 
only  slightly  affected  by  roughness  for  k/ho  <  0.06,  and  is  virtually  unchanged  for 
roughness  levels  above  this.  It  is  also  essentially  unaltered  by  the  reduced  jet  velocity  for 
0.13/1600  (k/ho  =  0.13/We  =  1600),  but  increases  noticeably  for  the  high  velocity  jet  on 
the  smooth-wall.  In  the  case  of  boundary  layer  flow  over  a  smooth  flat  surface,  the 
thickness  of  the  turbulent  boundary  layer  as  predicted  by  Schlichting  (1978)  is. 
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Substituting  Uo  for  Uoo,  and  local  mean  jet  thickness,  h ,  for  the  boundary  layer  thickness, 
6,  Eq.  (6)  is  seen  to  over-estimate  the  distance  to  this  transition  for  runs  0.00/3000  and 
0.00/4700  by  about  60  percent  and  45  percent,  respectively.  Davies  (1972)  credits  van 
der  Hegge-Zijnen  (1928)  with  deriving  the  following  equation  for  predicting  the  turbulent 
boundary  layer  thickness  over  a  rough-wall: 

(yiT 

5  =  0.26-kj^^J  .  (7) 

As  pointed  out  by  Davies,  Eq.  (7)  suggests  that  8  is  no  longer  a  function  of  Uoo.  Again, 
h  is  substituted  for  5  in  an  attempt  to  predict  the  transition  point  for  the  free  surface  of 
the  rough-wall  jets.  Just  as  was  the  case  for  the  smooth-wall,  Eq.  (7)  over-estimates  these 
transition  locations  by  25  to  100  percent. 

From  the  leading  edge  of  the  plate  to  the  beginning  of  the  surface  transition 
region,  the  temporal  mean  jet  thickness  increases  by  15  to  50  percent,  with  an  average 
increase  of  about  40  percent  (Fig.  23).  By  the  end  of  this  very  short  region,  the  thickness 
is  about  50  to  100  percent  greater  than  ho. 


In  the  second  region,  filament  formation,  those  disturbances  that  are  going  to 
produce  filament-like  structures  become  apparent.  As  suggested  in  Fig.  21.b,  the  length 
of  the  filaments,  as  well  as  their  frequency,  tend  to  increase  as  the  flow  progresses 
downstream.  Figure  22  shows  that  this  region  appears  to  be  more  strongly  influenced  by 
roughness  than  was  the  point  of  transition  to  roughness.  It  is  also  apparent  that  the 
potential  energy  loss  incurred  by  a  jet  flowing  over  an  angled  or  curved  wall  is  not 
significant  enough  at  this  point  to  affect  the  location  of  the  beginning  of  filament  growth. 
However,  changes  in  the  inertia  of  the  jet  do  seem  to  have  a  notable  impact  on  the 
location  of  this  event.  For  the  coarse- wall  (i.e.,  k/ho  =  0.13),  increasing  the  flow  rate 
reduces  the  distance  to  filament  formation  from  47ho  to  31  ho,  while  for  the  smooth-wall 
a  similar  increase  reduces  this  distance  from  105ho  to  83ho.  It  is  apparent  in  Fig.  24  that 
there  is  more  variation  in  the  mean  jet  thickness  at  the  beginning  of  the  filament 
formation  region  than  in  the  transition  to  roughness  region.  This  is  a  result  of  the 
continued  growth  in  size  and  quantity  of  the  free  surface  disturbances. 

The  drop  generation  region  is  the  most  dynamic  and  is  the  one  of  primary  interest 
in  this  research.  As  shown  in  Fig.  19,  the  key  elements  of  this  region  are  the  necking  of 
the  filament,  and  the  drop  separation  (pinch-off)  that  ensues  from  it.  As  indicated  in  Fig. 
20,  the  distance  from  the  tip  of  the  filament  to  the  necked  region  is  called  the  disturbance 
length,  X,  and  more  often  than  not,  a  single  drop  forms  from  a  single  disturbance.  Fig.  17 
shows  this  most  common  drop  separation  event,  but  there  are  others.  For  example,  the 
filament  may  separate  near  its  base  forming  a  large  elongated  drop  that  either  contracts 
into  a  spherical  drop,  or  breaks  up  into  two  or  more  satellite  drops.  The  filament  may 
also  disintegrate  into  multiple  drops  of  various  sizes — ^typically  two  or  three,  and 
sometimes  as  many  as  six.  Occasionally,  a  long-necked  region  will  produce  a  smaller 
satellite  drop  as  seen  in  free  jets  (Goedde  and  Yuen  1970).  The  subsequent  drop  analysis 
presented  here  will  concentrate  on  the  characteristic  behavior  of  these  tip-formed  drops, 
because  the  vast  majority  of  filaments  (=80  percent)  form  a  single  drop  at  the  tip. 

The  production  of  drops  in  the  drop  generation  region  depends  on  k/ho  and  jet 
velocity,  as  was  the  case  in  the  filament  formation  region.  Figure  22  indicates  that  at 
higher  values  of  k/ho,  the  dependence  of  Xj/ho  (location  of  the  beginning  of  drop 
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generation)  on  We  is  weak,  but  below  k/ho  =  0.06  it  is  unmistakable.  For  run  0.00/3000, 
drops  are  nonexistent,  however  at  Uo  =  7.7  m/s  (run  0.00/4700)  drops  begin  to  appear  at 
x/ho  *  130.  Although  the  apparent  influence  of  the  jet  inertia  on  the  beginning  of  drop 
generation  in  the  coarse-wall  is  not  as  pronounced  as  it  was  in  the  filament  formation 
region,  it  is  still  statistically  significant.  As  was  the  case  in  other  regions,  the  wall  angle 
and  curvature  have  no  noteworthy  influence  on  the  location  of  this  transition. 

Figure  25  provides  a  comparison  of  the  temporal  mean  jet  thickness  for  each  of 
the  runs  at  the  point  of  drop  generation.  The  variation  associated  with  h/ho  at  the 
beginning  of  this  region  is  slightly  less  than  that  seen  in  the  previous  two  regions,  and  the 
values  are  all  very  similar  (h/ho  «  2.0)  despite  the  large  difference  in  k/ho  and  Xd/ho. 
Nevertheless,  it  still  suggests  that  run  0.13/1600  gives  rise  to  a  jet  that  is  slightly  thicker 
when  it  begins  to  produce  drops.  A  difference  in  h  between  the  angled-  and  flat-wall 
jet — ^both  of  which  are  at  k/ho  =  0.06  and  We  =  3000 — also  becomes  apparent,  with  the 
angled-plate  mean  jet  thickness  being  slightly  less  when  drops  begin  to  form.  It  is 
noteworthy  that  the  value  of  h/ho  where  the  mean  jet  thickness  for  the  smooth  non-drop- 
producing  jet  levels  off,  the  other  flows  are  just  beginning  to  produce  drops  ( h/ho  « 1 .9 ). 

The  last  region  of  decreasing  activity  (not  shown  in  Fig.  21)  is  marked  by  a 
gradual  reduction  in  the  size  and  number  of  three-dimensional  surface  disturbances  and 
filaments.  The  onset  of  this  region  blends  smoothly  with  that  of  drop  generation,  and  as 
such,  its  exact  beginning  cannot  be  readily  identified.  It  appears  as  a  gradual  relaxation 
of  the  jet  free  surface.  The  jet  has  apparently  lost  the  prerequisites  (kinetic  energy, 
surface  energy,  critical  depth,  etc.)  to  produce  filaments.  An  example  of  the  character  of 
this  region  is  provided  in  Fig.  26  that  was  taken  fi'om  the  axisymmetric  jet  apparatus. 

2.  Jet  Thickness 

The  temporal  mean  jet  thickness  increases  as  it  moves  through  each  of  these 
transition  regions.  Figure  27  shows  how  this  thickness  varies  as  a  function  of  x/ho  for 
each  of  the  runs  considered.  For  runs  other  than  0.00/3000,  0.02/3000  and  0.13/1600,  the 
instantaneous  surface  fluctuations  are  very  large  beyond  x/ho  =  60  (or  x/ho  «  135  for  the 
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high-We  smooth-wall  run).  As  the  standard  deviation  values  for  the  jet  thickness 
suggest.  Fig.  28,  these  large  and  rapid  variations  in  the  free  surface  elevation  make  the 
present  measurement  technique  inaccurate  beyond  the  points  presented  in  Fig.  27. 
Furthermore,  it  is  not  clear  as  to  why  h  at  a  given  location  should  be  more  meaningful 
than  the  spectrum  of  the  fluctuations  at  the  same  point.  It  is  believed  that  the  violent 
nature  of  the  jet  free  surface  was,  in  large  part,  responsible  for  the  ineffectiveness  of  the 
capacitance  probe  in  measuring  h . 

At  first  glance  it  appears  that,  with  the  exception  of  run  0.00/4700,  all  of  the  jet 
thickness  profiles  shown  in  Fig.  27  behave  in  a  very  similar  manner:  a  gradual  increase 
in  h  followed  by  a  steep  change,  and  then  a  return  to  a  gradual  increase  (and  eventually 
one  would  expect  the  flow  to  revert  to  a  subcritical  state,  i.e.,  form  a  hydraulic  jump). 
However,  run  0.00/4700  is  not  truly  different,  it  is  simply  an  expanded  version  of  the 
other  curves.  The  region  from  0  <  x/ho  <  100  is  the  gradually  increasing  profile,  with  the 
steep  change  beginning  around  x/ho  =  110.  Recall  from  Fig.  22  that  the  transition 
regions  for  run  0.00/4700  are  substantially  further  out  than  that  of  the  others,  so  it  is  not 
unexpected  that  these  h  variations  would  occur  over  a  longer  distance. 

ff  one  compares  the  shape  of  h/ho  between  the  runs,  it  is  apparent  that  below 
some  level  of  relative  roughness  between  0.02  and  0.06,  the  aforementioned  thickness 
gradients,  dh/dx,  become  smaller.  Conversely,  for  k/ho  >  0.06,  the  initial  rate  of  jet 
thickening  is  rather  insensitive  to  the  variation  in  wall  roughness. 

D.  FILAMENT  AND  DROP  CHARACTERISTICS 

In  the  drop  generation  region,  a  minimum  of  30  filament/drop  events  were 
measured  in  most  of  the  x/ho  regions  recorded.  Near  the  end  of  the  wall,  the  flow  for 
most  runs  was  well  into  the  decreased  activity  region,  therefore  measurements  were  not 
typically  carried  out  beyond  x  =  130  cm.  Those  runs  that  did  continue  generating  drops 
beyond  x  =  130  cm  did  so  at  a  very  low  rate  (the  only  exception  was  run  0.00/4700, 
which  did  not  begin  generating  drops  until  x  =  75  cm).  As  the  end  of  the  wall  was 
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reached,  the  influence  of  secondary  effects,  such  as  splashing  from  the  end  of  the 
experimental  enclosure,  became  significant. 

Although  all  of  the  drops  recorded  were  measured,  the  present  discussion  will 
focus  on  those  formed  at  the  filament  tips.  As  noted  here  and  in  other  studies  (Volkert 
1980),  the  majority  of  the  filaments  give  birth  to  only  one  drop  at  the  filament  tip.  Even 
for  those  that  form  more  than  one,  the  "tip  drop"  is  almost  always  the  first  and  the  largest. 
The  few  drops  not  produced  by  this  process  were  very  small  (D*  ~  0.25  mm),  and  their 
formation  appeared  to  be  similar  in  nature  to  that  of  free  surface  jetting  (Rein,  1996). 

There  are  a  number  of  additional  measurements  that  could  have  been  taken  from 
the  digitized  images.  However,  here  the  analysis  is  directed  toward  those  measurements 
that  appeared  to  be  the  best  indicators  of  the  drop  generation  mechanism  to  be  discussed 
in  detail  later.  Those  indicators  are  the  drop  size,  velocity  and  relative  position;  and,  the 
filament  length,  diameter,  angle  and  tip  velocity  just  before  drop  pinch-off  (in  most  cases 
this  length  is  synonymous  with  the  maximum  filament  length).  It  is  felt  that  these 
measurements  give  the  clearest  possible  indication  of  the  events  causing  the  drop 
formation,  as  well  as  the  energy  associated  with  their  vertical  motion.  Only  the 
measurements  of  velocities  and  turbulence  within  the  jet  body  could  have  made  the 
mechanisms  of  filament  and  drop  formation  clearer. 

In  the  ensuing  discussion,  histograms  will  be  used  frequently  to  compare  and 
contrast  the  results  of  the  six  runs  listed  in  Table  3.  Although  the  smooth-wall  run 
(0.00/4700)  was  at  a  higher  flow  rate,  it  will  be  grouped  with  the  three  flat  rough-wall 
runs.  A  separate  comparison  between  the  high-flow  (0.13/3000)  and  low-flow 
(0.13/1600)  coarse- wall  case  will  be  used  to  help  elucidate  the  influence  of  jet  velocity. 
The  effects  of  wall  curvature  are  also  discussed  separately. 

1.  Drop  Diameter 

Figure  29  is  a  histogram  showing  the  drop-diameter  distribution  of  the  four  flat- 
wall  runs.  It  indicates  that  as  the  roughness  of  the  plate  is  increased,  the  mean  value  of 
Ddi/ho  becomes  slightly  larger  and  the  distribution  broadens  out.  The  variation  of  the 
mean  value  is  clearer  in  Fig.  30.  The  uncertainty  in  Dd^ho,  as  well  as  that  of  the  other 
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normalized  mean-value  characteristics  of  the  filaments  and  drops  presented,  is  less  than 
or  equal  to  the  height  of  the  corresponding  symbols  shown  in  the  figures  (unless 
otherwise  stated).  The  difference  in  Ddi/ho  is  almost  insignificant  among  the  rough- 
plates,  but  Ddi/ho  is  smaller  for  the  smooth-wall  case.  Reducing  the  coarse-wall  Weber 
number  from  3000  to  1600  has  no  apparent  influence  on  the  mean  drop  diameter  size,  but 
it  does  increase  the  production  of  smaller  drops  (Fig.  31).  Figure  32  presents  a 
comparison  of  the  flat  and  curved  plate  distributions  for  k/ho  =  0.06.  The  shape  of  the 
distributions  are  essentially  the  same,  but  the  curved-plate  values  are  shifted  to  the  left. 
This  is  consistent  with  its  lower  mean  value  in  Fig.  30.  In  short,  the  drop  diameter 
relative  to  the  initial  jet  thickness  increases  slightly  with  roughness  and  it  does  not  appear 
to  depend  on  We  (at  least  for  coarse  roughness).  Ddi/ho  is  smaller,  as  expected,  for  the 
smooth-plate;  and,  for  the  same  medium  roughness  (k/ho  =  0.06)  is  also  somewhat 
smaller  for  a  concave  curved-plate.  The  latter  may  be  due  to  the  inhibiting  effects  of 
centrifigal  acceleration  on  the  growth  of  the  filaments  from  which  the  drops  are  bom. 

2.  Filament  Length  and  Angle 

The  distribution  of  the  filament  length  just  before  drop  pinch-off,  Fig.  33,  exhibits 
a  more  apparent  skew  relative  to  that  of  the  drop  diameter.  As  the  relative  roughness 
increases,  the  distribution  broadens  and  shifts  to  the  right,  while  the  influence  of  concave 
curvature  moves  the  distribution  to  the  left  (Fig.  34).  It  is  evident  from  Fig.  35  that  both 
curvature  (i.e.,  centrifugal  acceleration)  and  jet  inertia  have,  as  expected,  a  substantial 
effect  on  the  maximum  length  of  the  filament  (the  former  decreases  it  while  the  latter 
increases  it).  Figure  36,  clearly  demonstrates  that  a  reduction  in  Weber  number  reduces 
the  likelihood  of  long  filament  formation.  Based  on  this  observation.  Fig.  35  implies  that 
if  the  rough-wall  had  been  sampled  at  the  same  Weber  number  as  the  smooth  (We  = 
4700),  one  might  have  seen  filaments  30  to  50  percent  longer.  Also,  the  medium 
roughness  on  a  flat-plate  would  have  yielded  smaller  filaments  lengths  at  We  =  1600 
relative  to  Lf/ho  =  2  at  We  =  3000.  However,  it  is  strongly  anticipated  that  there  are 
upper  limits  to  the  filament  size  as  We  and  k/ho  individually  or  simultaneously  increase. 
Although  not  measured  in  this  investigation,  it  is  clear  from  numerous  observations  (i.e.. 
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see  Fig.  37)  that  increasing  roughness  or  flow  rate  causes  dramatic  increases  in  the  drop 
and  filament  density  (drops/cm^). 

Changes  in  the  flow  inertia  or  plate-surface  condition  over  the  ranges  considered 
did  not  appear  to  influence  the  angle  of  the  filaments.  For  all  of  the  runs,  the  mean 
filament  angle  was  approximately  81°  (measured  CCW  from  the  negative  x-direction  as 
shown  in  Fig.  20).  In  other  words,  most  of  the  filaments  leaned  slightly  backwards  as 
they  moved  downstream,  as  in  Figs.  2  (Hoyt  and  Taylor,  1977a)  and  19.  Additionally, 
for  a  random  sampling  of  196  filaments,  the  change  in  angle  from  the  moment  a  filament 
entered  the  viewing  area  of  the  camera  until  just  before  drop  release  was  measured  at  +4°. 
This  implies  that  there  is  no  apparent  change  in  the  mean  filament  angle,  because  the 
angle  estimation  is  only  accurate  to  within  about  ±5°. 

3.  Effective  FUament  Diameter  and  Slenderness  Ratio 

Figure  38,  shows  the  distribution  of  the  effective  filament  diameter,  Df/ho.  The 
distribution  widens  and  moves  toward  larger  values  of  Df/ho  as  roughness  increases,  but 
only  marginally  so  for  k/ho  from  0.02  to  0.06.  Although  the  change  is  larger  between 
k/ho  =  0.06  and  0.13,  it  is  still  not  as  significant  as  the  difference  between  the  smooth  and 
fine-roughness  walls.  Comparing  these  changes  with  that  between  the  high-  and  low- 
Weber  number  coarse-wall  flows.  Fig.  39,  it  is  apparent  that  the  lower  flow  rate  shifts  the 
filament  diameter  distribution  slightly  towards  smaller  values  of  Df/ho.  Looking  at  Fig. 
40,  it  is  seen  that  the  form  of  the  variation  of  D^o  is  similar  to  that  of  the  Drop  diameter 
(Fig.  30).  Recall  that  D^o  is  estimated  by  dividing  the  filament  cross-sectional  area  (as 
determined  using  the  Optimas  MA  software)  by  the  filament  length.  Although  the 
uncertainty  analysis  shows  a  great  deal  of  uncertainty  for  the  individual  filament 
diameters,  the  mean  values  remain  fairly  good  predictors  of  the  overall  population.  The 
influence  of  curvature  on  the  distribution  is  shown  in  Fig.  41.  It  indicates  that  the  shape 
is  largely  the  same,  but  the  distribution  is  shifted  to  lower  values  of  Df/ho.  Again,  this  is 
very  similar  to  the  drop  diameter  behavior  noted  in  Fig.  32,  and  the  same  summary 
statement  applies. 
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Having  examined  both  the  filament  length  and  diameter,  it  is  now  worthwhile  to 
consider  the  filament  slenderness  ratio,  Lf/Df,  just  prior  to  drop  formation.  Figure  42 
shows  the  Lf^f  distribution  for  the  four  nominal  runs.  The  variance  is  extremely  small 
for  the  smooth-wall,  but  gradually  increases  with  roughness.  The  mean  value  of  Lf/Df  in 
Fig.  43  also  increases  as  roughness  increases.  The  25  percent  reduction  in  the 
slenderness  ratio  for  the  low-Weber  number  coarse-wall  case  suggests  that  the  smooth- 
wall  value  of  4.1  is  somewhat  enhanced  by  the  higher  Weber  number  of  that  flow.  As 
shown  in  Fig.  44,  increasing  the  velocity  of  the  jet  broadens  the  distribution  in  addition  to 
increasing  the  mean  value.  The  influence  of  curvature,  Fig.  45,  is  similar  to  what  was 
noted  before:  the  variance  of  the  distribution  and  the  mean  value  are  reduced  a  small 
amount.  As  noted  earlier,  the  said  increases  with  k/ho  and/or  We  are  expected  to  have 
reasonable  upper  bounds,  perhaps  accompanied  by  new  transformations  such  the  total 
disintegration  (atomization)  of  the  wall  jet. 

4.  Filament  Volume 

As  previously  explained,  the  filaments  are  idealized  as  cylindrical  columns  (Fig. 
20).  This  idealization  can  be  used  to  consider  a  number  of  filament  geometric  properties; 
e.g.,  a  comparison  of  the  effective  diameter,  angle,  disturbance  wave  number  and  volume 
of  the  filaments.  Such  calculations  should  provide  reasonable  results  because  most  of  the 
drop-forming  filaments  were  axisymmetric  in  nature,  and  only  those  which  clearly 
exhibited  this  behavior  were  considered  in  the  study.  For  the  case  of  the  filament 
volume,  one  has 

Vf 

4 

The  uncertainty  associated  with  an  individual  filament  volume  estimate  is  large,  but  the 
mean  value  is  still  good  because  it  is  based  on  about  100  to  150  samples  (see  Appendix 
B).  The  small  uncertainty  regions  of  Fig.  46,  which  displays  the  mean  nondimensional 
filament  volumes,  Vf/ho^,  makes  this  point  clear.  Although  the  large  measurement 
uncertainty  (about  40  percent)  may  significantly  add  to  the  variance  of  the  distributions 
shown  in  Figs.  47  through  49,  they  are  still  of  some  value.  For  example,  even  with  the 
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uncertainty  tending  to  contribute  a  large  amount  of  scatter,  the  plots  still  suggest  that  the 
distribution  is  rather  narrow  compared  to  the  other  measured  parameters.  From  Fig.  46, 
the  volume  appears  to  decrease  linearly  with  roughness  for  k/ho  >  0.02,  but  falls  off  more 
sharply  for  a  value  of  k/ho  less  than  this.  Comparing  the  mean  drop  volume  based  on  the 
Ddi/ho  values  provided  earlier  (Fig.  30)  with  the  filament  volumes  shown  here,  it  is  seen 
in  Fig.  50  that,  on  average,  a  drop  consumes  about  50  to  60  percent  of  the  filament 
volume  at  pinch-off.  The  figure  suggests  that  there  may  be  a  slight  decrease  in  the 
volumetric  ratio,  Vdi/Vf,  with  increasing  k/ho  and  We.  Nonetheless,  it  seems  apparent 
that  the  small  change  in  drop  diameter  with  k/ho  allows  the  drop  to  maintain  a  fairly 
constant  volumetric  ratio  with  the  filament,  despite  significant  changes  in  filament 
volume. 


5.  Disturbance  Wave  Number 

Consideration  is  now  given  to  the  disturbance  wave  number,  just  prior  to  drop 
pinch-off. 


The  individual  wave  numbers  carry  an  uncertainty  of  about  20  percent  because  the 
determination  of  the  wavelength,  X,  requires  more  subjective  interpretation  than  some  of 
the  other  parameters.  It  is  expected  that  this  uncertainty  is  random  in  nature,  so  that  the 
mean  value  calculations  still  provide  reasonable  estimates.  Figure  51  shows  that  the 
mean  values  of  ki  do  not  display  a  clear  trend,  and  that  the  average  value  for  flat-plate 
flows  is  about  1.8.  There  was  also  no  consistent  behavior  seen  in  the  individual  filament 
wavelengths.  During  the  life  of  a  filament,  some  disturbances  grew,  some  shrank  and 
some  were  unidentifiable  until  just  before  drop  pinch-off.  However,  ki  does  display  a 
fairly  constant  mean  value  and  a  Gaussian-like  distribution  (Fig.  52  through  54)  for  all 
wall  jets  considered. 
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6.  Drop  and  Filament  Characteristics  as  a  Function  of  Streamwise 

Position 

All  of  the  physical  attributes  of  filaments  and  drops  were  measured  in  different 
areas  of  the  drop  generation  region.  The  above  discussion  considered  the  region  as  a 
whole,  but  it  may  now  be  useful  to  consider  how  these  values  vary  from  area  to  area 
within  the  region.  Figures  55  through  59  show  only  those  locations  where  30  or  more 
observations  were  made.  In  addition,  all  abscissa  values  are  plotted  in  terms  of 
(x  -  Xj  )/ho  to  make  comparisons  between  runs  clearer. 

The  local  mean  drop  diameter,  Fig.  55,  shows  a  slight  increase  in  drop  size  in  the 
downstream  direction.  With  the  exception  of  the  fine-roughness  (0.02/3000)  and  low 
velocity  (0.13/1600)  runs  where  the  drop  size  is  decreasing,  or  beginning  to  decrease,  as 
Ax/ho  increases.  Regardless,  the  change  is  very  slight  and  the  drop  size  remains  on  the 
order  of  h(x)  throughout  the  region. 

The  behavior  of  the  mean  filament  length  with  position.  Fig.  56,  is  fairly  constant 
given  the  level  of  uncertainty.  However  unclear,  there  is  perhaps  a  mild  downward  trend 
for  the  fine-roughness  and  the  low-Weber  number  flows.  The  large  value  of  Lf/ho  seen  at 
the  end  of  the  coarse-wall  flow  is  a  result  of  the  existence  of  a  few  rather  large  filaments, 
and  the  uncertainty  of  this  value  is  actually  closer  to  20  percent.  For  the  mean  filament 
diameter.  Fig.  57  indicates  that  these  values  are  also  reasonably  constant  over  the  length 
of  x/ho  considered,  exhibiting  a  character  similar  to  Da/ho  (Fig.  55).  The  trend  of  the 
curved-wall,  for  which  Df/ho  is  consistently  smaller,  agrees  with  the  mean  value  results 
of  Fig.  40. 

In  Fig.  58,  which  depicts  the  mean  local  slenderness  ratio  of  the  filaments  for  each 
run,  the  uncertainty  is  up  to  about  15  percent,  so  one  should  be  cautious  about  reading  too 
much  into  the  mild  trends  which  exist.  However,  the  consistency  of  the  reduction  in 
Lf/Df  in  the  streamwise  direction  for  all  runs  suggests  that  the  trend  is  probably  real. 
Finally,  for  completeness.  Fig.  59  presents  the  mean  local  filament  volume.  The 
uncertainty  is  now  in  the  25  percent  range  (and  as  high  as  30  for  some  of  the  larger 
values).  The  curves  generally  agree  with  the  mean  filament  volume  of  Fig.  46:  an 
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increase  in  inertia  or  roughness  tends  to  increase  the  volume  of  the  filament  just  before 
drop  pinch-off. 

In  summary  of  the  foregoing,  one  may  reasonably  argue  that  the  basic 
characteristics  of  the  drops  and  filaments  and  their  volumes  remain  relatively  constant 
over  the  area  of  measurement,  with  small  and  often  statistically  insignificant  variations, 
along  the  streamwise  direction.  Having  said  this,  there  is  undoubtedly  some  transitional 
variations  in  the  said  characteristics  at  the  beginning  and  the  end  of  the  measurement 
regions.  Additionally,  any  variations  are  obscured  or  masked  by  their  Gaussian 
distribution  at  any  given  section  of  the  area. 

E.  WATER  TUNNEL  WALL  JET 

The  previous  sections  examined  the  different  regions,  and  the  filament  and  drop 
characteristics,  of  a  rectangular  liquid  wall  jet.  It  was  not  possible  to  make  measurements 
within  the  jet  because  of  its  extremely  violent  nature  and  shallow  depth.  The  goal  of  the 
water  tunnel  investigation  was  to  help  gain  some  insight  into  the  nature  of  the  flow  below 
the  free  surface.  The  tools  of  choice  for  this  investigation  were  a  digital  particle 
velocimetry  system  (DPIV)  and  a  Redlake  High-Speed  Digital  Imager  (Table  2). 

It  was  possible  to  produce  a  high-speed  jet  with  velocities  equal  to  that  of  the 
rectangular  wall  jet  by  arranging  the  test  section  of  the  water  tunnel  as  described  in 
section  n.D.  Although  the  tunnel  had  the  ability  to  support  a  wide  variety  of  wall  jet 
parameters,  there  were  explicit  limits  on  the  choices  for  ho  and  Uo.  If  ho  was  selected 
too  large,  then  the  experiment  would  leave  the  realm  of  wall  jets,  becoming  a  subcritical 
or  mildly  supercritical  open  channel  flow  (i.e.,  with  a  relatively  undisturbed  free  surface). 
If  ho  were  too  small,  then  it  would  be  impractical  to  visualize  the  interior  of  the  jet  on  a 
scale  providing  any  sort  of  useful  insight.  For  low  values  of  Uo,  the  flow  again  becomes 
more  of  an  open  channel  flow.  However,  for  large  values  of  Uo  for  a  given  ho  (i.e.,  high 
Froude  number),  the  free  surface  is  excessively  disturbed  and  essentially  destroys  the 
optical  path  between  the  side  of  the  test  section  and  the  center  of  the  wall  jet.  High  flow 
rates  at  large  depths  were  also  of  no  use  because  the  surface  transitions  that  are 
distinctive  of  wall  jets  would  not  occur  within  the  short  length  of  the  test  section.  Thus,  it 
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was  necessary  to  achieve  the  best  possible  flow  circumstances  in  order  to  get  a  glimpse  of 
what  might  be  transpiring  below  the  free  surface.  Therefore,  a  single  Uo-ho  setting  that 
would  meet  all  the  necessary  conditions  was  found  after  considerable  trial  runs.  It  was 
determined  that  a  vertical  plate  height  of  1  l-mm  with  the  pump  providing  a  flow  rate  of 
approximately  3  m/s,  produced  a  wall  jet  with  ho  =  9.5  mm  at  the  vena-contracta.  The 
resulting  jet  was  similar  to  run  0.00/3000  when  the  wall  was  smooth  (i.e.,  no  drops 
formed),  but  was  capable  of  producing  a  few  filaments  and  drops  near  the  end  of  the  test 
section  when  the  wall  was  roughened  with  0.72  mm  sand  particles.  Table  4  lists  the 
parameters  of  the  two  water  tunnel  wall  jets. 


Plate 

k 

(mm) 

Umax 

(m/s) 

k 

ho 

(Rex/x)xl0'^ 

(m-') 

Rex  10-^ 

Fr 

We 

ZxlO^ 

Smooth 

0.00 

3.5 

0.00 

3.6 

3.4 

11 

1600 

1.2 

0.72 

2.9 

0.08 

2.8 

9 

1100 

1.2 

Table  4.  Summary  of  the  water  tunnel  liquid  wall  jet  parameters. 


Figure  60  provides  an  idea  of  how  the  water  tunnel  rough-wall  jet  behaved  in 
comparison  to  the  other  jets.  It  reinforces  the  fact  that  the  free  surface  transition  point  is 
relatively  inert  to  changes  in  k/ho  and  velocity.  Considering  the  large  difference  in 
Weber  number  between  this  and  the  previous  rectangular  wall  jet  (1 100  versus  3000),  it 
is  somewhat  surprising  that  the  transition  to  filament  formation  is  not  further  out.  The 
beginning  of  the  drop  generation  region  changes  substantially,  and  is  similar  to  that  of 
run  0.02/3000.  The  smooth-wall  jet  did  undergo  transition  to  roughness  on  the  free 
surface,  but  the  process  was  too  gradual  for  the  location  of  the  event  to  be  determined. 

The  mean  velocity  profile  of  the  jet  was  measured  over  both  smooth  and  rough 
walls  using  a  DPIV  system.  Figure  61  shows  the  results  of  these  measurements  along 
with  those  of  Klebanoff  (1954)  in  a  fully-developed  smooth-wall  boundary  layer,  and 
Grass  (1971)  in  a  fully-rough  channel  flow.  The  jet  velocities  were  recorded  at  x/ho  = 
50,  which  is  roughly  the  beginning  of  the  filament  formation  region. 

The  visualization  of  the  tunnel  wall  jet  was  conducted  using  the  Redlake  high¬ 
speed  camera  positioned  along  side  the  wall  jet.  The  interior  of  the  jet  was  illuminated 
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with  a  diffused  backlight,  or  from  directly  beneath  the  center  of  the  test  section  using  the 
DPIV  laser  light  sheet  as  an  ultra-high-speed  two-dimensional  flash.  The  visualization 
results  for  the  smooth-wall  flow  were  of  no  interest  and  will  not  be  presented. 
Conversely,  the  subsurface  flow  over  the  rough-wall  was  quite  fascinating,  as  seen  in  Fig. 
62.  Here  one  sees  the  light  sheet  highlighting  a  number  of  bubbles  that  have  been  drawn 
into  the  flow.  The  bottom  of  the  images  show  the  roughness  and  the  dark  areas  at  the  top 
are  a  result  of  the  time-dependent  free  surface  variations  existing  along  the  optical  path. 
The  flow  in  the  images  is  from  left-to-right.  Figure  62.a  presents  an  elongated  bubble 
that  has  not  separated  from  the  top,  yet  it  is  nearly  in  contact  with  the  rough-wall.  In  Fig. 
62.b  there  is  a  similar  type  of  event,  except  that  the  photograph  is  taken  shortly  after  the 
bubble  has  pinched  off  from  the  free  surface.  A  number  of  bubbles  of  different  size  and 
shape  are  shown  in  Fig.  62.c.  The  size  of  the  bubbles  generally  vary  from  <  0.25mm  to 
one-half  of  the  jet  thickness,  but  in  most  cases  they  appear  to  be  about  one-fourth  to  one- 
sixth  of  h(x).  A  few  of  the  smaller  bubbles  resulted  from  air-entrainment  at  the  end  of  the 
test  section  (and  subsequently  recirculated  in  the  system),  but  most  were  seen  to  be 
introduced  by  free  surface  disturbances  as  shown.  Some  of  the  larger  bubbles  in  Fig.  62 
suggest  that  there  may  be  a  preferential  angular  orientation  (about  125°  to  135°).  Figure 
63  is  similar  to  Fig.  62,  except  it  is  taken  with  the  jet  back-lit  and  the  camera  angled 
upward  to  get  a  clearer  view  of  the  free  surface.  These  images  were  captured  just 
downstream  of  the  position  where  the  velocity  profile  was  acquired.  It  is  the  appearance 
of  the  bubbles  in  this  region  (dark  blobs  in  Fig  63  .b)  that  prevented  the  recording  of 
DPIV  data  at  points  further  downstream  than  x/ho  «  50,  as  it  was  not  clear  what  influence 
the  bubbles  would  have  on  the  DPIV  velocity-vector  calculations.  Figure  64  is  a 
sequence  of  back-lit  images  showing  how  bubbles  are  introduced  into  the  flow.  The  free 
surface  fluctuations  in  this  region  were  fairly  moderate  (same  region  as  Fig.  63),  however 
this  sequence  leaves  little  doubt  that  these  disturbances  have  a  profound  impact  on  the 
character  of  the  jet.  As  the  bubbles  move  downstream,  most  never  make  it  back  to  the 
free  surface  prior  to  reaching  the  end  of  the  test  section  (x/ho  «  1 10).  The  bubbles  rotate, 
stretch  and  dart  about  randomly  and  rapidly,  and  some  of  those  that  reach  the  bottom  will 
bounce  along  the  wall  as  they  move  downstream.  Although  the  stretching  (Figs.  62.a  and 
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64)  is  a  clear  indication  that  the  bubbles  are  undergoing  shear,  there  was  absolutely  no 
evidence  that  they  felt  any  type  of  drag  from  the  liquid  flow:  the  bubbles  neither 
collapsed  nor  distorted  in  any  sort  of  preferential  way  along  their  trajectories.  Figure  65 
is  a  unique  photograph  that  simultaneously  captures  a  filament  above  the  free  surface  and 
bubbles  below. 

The  mean  velocity  profiles  presented  earlier  are  the  result  of  over  450  horizontally 
averaged  DPIV  image-pair  vector  plots.  The  root-mean-squared  velocity  values  were 
also  calculated  from  this  data  set,  Fig.  66,  but  over  20,000  vector  plots  would  be 
necessary  to  make  them  statistically  significant.  Notwithstanding  the  large  uncertainties 
that  are  recognized  to  exist  (and  which  may  explain  the  anomalies  present  in  the  v^s 
figure),  these  results  are  presented  for  what  they  are  worth. 

Instantaneous  v-velocity  contour  plots  of  the  DPIV  data  taken  at  x/ho  »  50  are 
shown  in  Fig.  67.  This  is  also  part  of  the  same  data  set  that  was  used  to  calculate  the 
mean  velocity  profiles  in  Fig.  61.  The  whited-out  regions  near  the  top  left  comer  are 
areas  where  insufficient  valid  vectors  were  calculated  (usually  due  to  surface  reflections). 
The  circled  region  on  each  plot  annotates  a  fluid  "clump"  with  a  vertical  velocity  greater 
than  0.5  m/s  (v'/U  «  0.18).  It  is  interesting  that  such  large,  seemingly  coherent,  fluid 
stmctures  with  significant  upward  velocities  exist  this  far  upstream.  Considering  their 
present  vertical  position  and  velocity,  as  well  as  the  velocity  of  the  jet,  all  three  of  these 
structures  should  reach  the  free  surface  zone  prior  to  x/ho  =  70. 

F.  DISCUSSION  OF  RESULTS 

The  most  important  contribution  of  this  investigation  is  the  quantification  of  the 
physical  manifestations  of  the  liquid  wall  jets.  The  aforementioned  effort  has  resulted  in 
numerous  individual  measurements  of  filament  length,  diameter,  angle,  disturbance 
wavelength  and  velocity;  drop  diameter(s)  and  velocity;  and,  the  mean  temporal  jet 
thickness.  The  reliability  of  the  information  is  related  to  the  number  of  individual 
realizations,  hence,  the  data  gathering  was  an  enormously  time-consuming  affair  (which 
perhaps  explains  why  it  has  not  been  done  before). 
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In  addition  to  providing  a  basis  for  identifying  the  governing  and  influential 
parameters  in  the  present  study,  these  values  will  be  essential  for  gaining  even  deeper 
insight  from  future  studies  of  this  nature.  There  is  no  doubt  that  the  structures  of  the  free 
surface  affect  the  flow  beneath,  and  the  underlying  characteristics  of  the  flow  affect  the 
structures  of  the  free  surface.  In  fact,  one  may  conclude  that  the  spray  problem  is  the 
determination  of  this  interaction  for  smooth,  rough  and  concave  or  convex  walls  in 
supercritical  wall  jets. 

1.  Free  Surface  Characteristics 

This  work  has  identified  four  separate  surface  transition  regions  that  were 
wimessed  in  all  of  the  drop-forming  flows.  The  first,  where  the  free  surface  disturbances 
begin  to  appear,  leads  to  the  second  region  where  some  of  the  disturbances  continue  to 
grow  into  vertical  columns  of  fluid  (or  filaments).  In  the  third  region,  some  of  the 
filaments  give  rise  to  drops  (usually  just  one).  The  number  and  size  of  the  drops  are 
strongly  influenced  by  the  jet  velocity  and  relative  wall  roughness,  among  other  things. 
The  last  region  is  one  of  decaying  free  surface  activity,  where  the  flow  is  still  highly 

supercritical  (i.e.,  FT  =  u/-yJ^»l),  but  is  below  some  minimum  threshold  necessary 

for  filaments  and  drops  to  be  produced.  As  this  region  progresses,  the  free  surface 
continues  to  calm,  becoming  reminiscent  of  shallow  channel  flows  with  their  boils  and 
kolks.  The  influence  of  roughness  on  the  location  of  the  transitions  was  most  obvious  at 
lower  roughness  levels.  Beyond  a  certain  value  of  k/ho,  the  influence  of  roughness  was 
more  apparent  in  the  quantity  of  the  free  surface  structures,  rather  than  the  location  of 
their  introduction. 

2.  Spray  Generation  Features 

It  is  critically  important  to  recognize  that  the  entire  drop  production  phenomena  is 
of  a  statistical  nature,  and  therefore  necessitates  a  statistical  analysis  for  its  interpretation. 
By  using  such  an  approach,  it  was  possible  to  ascertain  many  of  the  fundamental 
attributes  associated  with  the  formation  of  the  filaments  and  drops.  Figure  68  provides  a 


47 


6 

Lf/Df 

3 


2 

Lf/ho 

1 

0.7 

Ddi/ho 

0.4 
0.4 

Vf/ho^ 

0.1 
0.5 

Df/ho 

0.2 
1 

Vd/Vf 

0 

2.5 

ki 

1.5 

_ 0 _ 0.05  k/ho  0.1 _ 0.15  [ 

Figure  68.  Summary  of  the  mean  value  characteristics  as  a  function  of  the  relative 
roughness. 
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summary  of  the  mean  value  of  those  characteristics  as  a  reference  for  the  following 
discussion. 

In  all  of  the  liquid  wall  jets  examined,  the  thickness  of  the  jet  increases 
substantially  as  the  jet  relinquishes  part  of  its  momentum  to  the  wall.  Roughness 
enhances  this  process,  and  as  such  leads  to  more  rapid  jet  thickening.  Perhaps  more 
importantly,  higher  levels  of  roughness  result  in  greater  temporal  mean  variations  in  the 
local  values  of  h,  even  to  the  point  where  the  measurement  technique  employed  here  is  no 
longer  viable. 

In  the  case  of  the  drop  diameters,  it  was  apparent  that  the  mean  value  remained 
essentially  constant  over  all  flat  rough-walls  regardless  of  roughness  level  (D*  «  0.58ho), 
but  that  the  smooth-plate  values  were  smaller  (Ddr  =  0.51).  The  distribution  of  the  drop 
size  about  the  mean  is  well-described  by  a  normal  distribution  curve,  where  the  standard 
deviation  of  the  curve  increases  with  k/ho.  This  behavior  conforms  with  that  previously 
noted  in  typical  boundary  layer  flows  by  Blake  (1970)  and  Krogstad,  et  al  (1992).  Blake 
found  that  the  larger  interstitial  region  present  around  larger  roughness  elements 
produced  a  greater  range  of  eddy  sizes  and,  similarly,  Krogstad  et  al  showed  that  the 
spectral  distribution  of  energy  containing  eddies  is  wider  in  rough-wall  flows  as 
compared  to  smooth.  Others  have  identified  the  Gaussian  nature  present  in  the  features 
which  describe  the  motion  of  the  large  turbulent  eddies  developed  in  boundary  layer 
flows  (e.g.,  Raupach,  1981). 

The  filament  length  distributions  presented  a  clear  trend  toward  the  lower 
filament  lengths,  which  is  different  than  that  described  for  the  aforementioned  drop 
diameter.  From  a  physical  viewpoint,  this  makes  sense  because  a  filament  must  have  a 
finite  length  to  form  a  drop  (Lf  >  Ddr),  and  very  large  filaments  would  necessitate  a 
turbulent  fluctuation  with  an  exceedingly  high  vertical  velocity.  These  physical  boundary 
conditions  act  to  limit  how  small  (or  large)  the  filaments  can  be,  but  the  lower  energy 
requirement  of  short  filaments  make  their  existence  more  probable.  The  lengths  of  the 
filaments  strongly  increase  with  roughness,  suggesting  that  roughness  enhances  the 
vertical  velocity  fluctuations  in  the  flow,  enabling  the  eddies  to  "stretch"  the  fluid 
continuum  to  a  greater  extent.  Again,  the  influence  of  roughness  on  this  parameter  is  not 
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unlike  that  seen  in  typical  turbulent  boundary  layer  flows:  increasing  k/ho  results  in  more 
frequent  and  higher  energy  quasi-coherent  ejections  of  fluid  from  the  wall  region  (Grass, 
1971;  Raupach,  1981). 

The  constancy  of  the  filament  disturbance  wavelength,  ki,  and  its  Gaussian 
distribution  suggest  that  it  is  what  it  appears  to  be:  an  indication  of  the  diameter  of  the 
drop  which  is  about  to  form.  Nonetheless,  its  mean  value  of  about  1.8  does  provide  some 
additional  evidence  into  the  nature  of  the  filament  formation  process.  Using  a  two- 
dimensional  linear  analysis,  Rayleigh  (1945)  demonstrated  that  ki  must  be  greater  than 
unity  for  the  disturbance  along  a  liquid  column  to  grow.  The  disturbance  of  maximum 
growth  rate  was  shown  to  be  ki  «  0.7,  while  values  greater  than  one  were  stable.  Even 
the  more  contemporary  solutions,  which  include  the  effects  of  viscosity  and  finite 
amplitude  disturbances,  do  not  deviate  from  these  values  greatly.  One  of  the  reasons  for 
this  is  that  the  non-linear  effects  of  disturbance  growth  do  not  become  important  until  the 
last  10  to  20  percent  of  the  length  of  the  jet  (McCarthy  and  Molloy,  1974).  However, 
unlike  jets  in  the  Rayleigh  regime,  the  filaments  examined  here  display  bumps,  wrinkles 
and  kinks,  and  ki  consumes  about  40  percent  of  the  entire  filament  length.  Thus,  it  is  not 
likely  that  their  character  can  be  described  by  a  two-dimensional  analysis,  regardless  of 
whether  it  is  linear  or  non-linear.  These  filaments  also  do  not  conform  to  any  of  the  other 
characteristic  descriptions  provided  by  the  jet  breakup  modes  mentioned  in  the 
introduction. 

Based  on  what  is  known  about  drop  behavior  in  shear  fields  (e.g.,  Clift,  et  al, 
1978),  it  is  evident  that  there  is  little  or  no  aerodynamic  influence  affecting  the  motion  of 
the  drops.  For  the  filaments,  the  best  quantitative  measure  of  aerodynamics  in  the  change 
in  ligament  angle  which  was  shown  to  be  negligible.  The  two  DPIV  experiments 
demonstrated  that  the  flow  of  the  air  above  the  free  surface  could  not  be  differentiated 
from  that  of  the  liquid  near  the  free  surface  zone.  Hence,  the  observational  data 
combined  with  the  statistical  results  suggest  that  there  is  little  or  no  aerodynamic 
influence  on  the  above-free-surface  character  of  the  wall  jet. 

Trends  in  the  geometric  data  (i.e.,  Lf,  Df,  L^f,  Ddr,  ki  and  Vf)  as  a  function  of  the 
streamwise  direction  were  also  considered,  but  the  degree  of  random  scatter  coupled  with 
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the  restricted  number  of  observations  at  each  location  limited  the  amount  of  information 
they  provided.  For  all  runs  except  two,  the  drop  and  filament  diameters  showed  a  slight 
increase  with  x.  The  two  exceptions  were  for  the  lowest  roughness  plate  and  the  lowest 
velocity  (coarse)  plate,  so  it  is  possible  that  their  results  portray  a  transitioning  below 
some  critical  nature. 

It  was  not  feasible  to  measure  the  properties  inside  of  the  rectangular  wall  jet  for 
the  reasons  mentioned  previously.  However,  by  using  the  known  filament  parameters,  it 
is  possible  to  work  backwards  to  get  some  idea  of  the  jet  characteristics  that  gave  rise  to 
their  existence.  For  example,  using  a  simple  energy  balance  that  accounts  for  kinetic, 
potential  and  surface  energy,  one  may  derive  the  vertical  velocity  needed  below  the  free 
surface  to  produce  a  given  filament: 

(8) 

Dividing  both  sides  of  Eq.  (8)  by  Uo  and  performing  a  small  amount  of  additional 
manipulation  leads  to 


where  Fr  and  We  are  the  Froude  and  Weber  numbers  of  the  jet,  respectively.  Equation 
(9)  may  be  used  to  calculate  the  eddy  fluctuation  necessary  to  give  birth  to  each  of  the 
individual  filaments.  This  results  in  mean  turbulent  fluctuations  that  vary  from  about  9 
percent  for  the  smooth-plate  to  nearly  17  percent  for  the  0.13/3000  case,  with  an  average 
value  of  v'/Uo  =  0.12. 

Perhaps  even  more  useful,  Eq.  (9)  can  be  used  to  derive  a  rough  approximation  of 
the  turbulent  fluctuations  necessary  for  the  inception  of  drop  formation.  By  applying 
average  values  for  local  mean  jet  thickness,  filament  length,  and  effective  filament 
diameter,  one  arrives  at  an  expression  for  the  estimated  turbulence  intensity  in  terms  of  Fr 
and  We  alone 
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which  is  similar  to  the  technique  used  by  Merrill  and  Sarpkaya  (1998).  As  plotted  in  Fig. 
69,  Eq.  (10)  more  clearly  demonstrates  the  relationship  between  the  sub-surface  velocity 
fluctuations  and  the  Froude  and  Weber  numbers.  The  importance  of  the  level  of 
supercriticality  of  the  jet  becomes  very  apparent  in  Fig.  69,  which  suggests  that  a  wall  jet 
with  a  "low"  level  of  supercriticality,  requires  very  large  values  of  v'/Uo  (especially  in  the 
outer  regions  of  the  flow)  to  generate  drop-forming  filaments.  It  was  previously  pointed 
out  that  the  relative  roughness  is  a  key  factor  in  drop  formation,  yet  it  does  not  appear 
explicitly  in  Eq.  (10).  The  reason  for  this  is  that  Df/ho  and  h/hoare  the  only  terms 
affecting  the  numerators  in  Eq.  (10)  (see  Eq.  (9)).  Df/ho  can  be  treated  as  essentially 
constant  (Fig.  68)  and  h/ho  does  not  depend  on  k/ho  for  drop  inception  (see  Fig.  25). 
Hence,  a  jet  that  is  sufficiently  supercritical  and  reaches  a  value  of  h/ho  =  2 ,  at  the  levels 

of  Fr  and  We  shown  in  Fig.  69,  embodies  all  of  the  features  necessary  to  produce  nascent 
drop-forming  filaments.  Roughness,  therefore,  does  not  specify  what  level  of  v'/Uo  is 
necessary  for  drop  inception;  rather,  it  increases  the  likelihood  that  a  given  flow  is 
capable  of  producing  that  level  of  v'/Uo.  Regardless  of  roughness,  however.  Fig.  69 
strongly  suggests  that  for  drops  to  be  realizable  events,  the  flow  must  have  reached  a 
minimum  supercritical  state. 

Although  the  above  expression  was  based  on  the  events  observed  at  the  free 
surface,  the  results  make  it  clear  that  the  nature  of  the  filaments  and  drops  is  inextricably 
linked  to  the  subsurface  nature  of  the  jet. 

3.  Internal  Flow  Structure  of  Bow  Sheets 

A  complete  understanding  of  the  drop  formation  phenomena  requires  not  only  a 
knowledge  of  what  exists  above  and  near  the  free  surface,  but  also  (unfortunately)  what  is 
beneath  it,  all  the  way  to  the  wall  and  the  internal  structures  that  form  on  it.  Such  a  total 
awareness  would  necessitate  a  comprehensive  edification  of  the  creation,  convection  and 
diffusion  of  turbulence  within  the  jet.  Recognizing  that  the  free  surface  disturbances 
serve  to  couple  what  is  above  with  what  is  below,  and  that  they  provide  a  path  for  each 
"world"  to  communicate  with  the  other,  knowledge  of  their  behavior  is  a  necessary — ^but 
not  a  sufficient — condition  for  describing  the  subsurface  behavior. 
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From  what  has  been  seen  and  measured,  it  is  evident  that  the  internal  flow 
structure  of  the  wall  jet  (or  bow  sheet)  is  not  the  same  as  either  a  typical  boundary  layer 
or  a  submerged  wall  jet.  The  velocity  profile  in  the  latter  is  largely  influenced  by  the 
strong  shear  at  the  periphery  of  the  jet  (in  addition  to  the  distance  from  the  nozzle, 
velocity  profile  and  turbulence  contained  within  the  jet  at  birth),  a  factor  not  present  in 
bow  sheets.  The  difference  between  this  jet  and  a  typical  turbulent  boundary  layer 
becomes  more  evident  when  one  attempts  to  predict  when  the  boundary  layer  in  the  wall 
jet  will  reach  the  free  surface  using  analytically  determined  equations  such  as  Eqs.  (6) 
and  (7).  The  free  surface  fluctuations  that  develop  in  the  drop  generation  region — ^which 
temporally  change  the  local  jet  thickness  by  as  much  as  50  percent — ^must  alter  the 
boundary  layer  character  in  the  outer  region,  since  the  latter  extends  to  the  free  surface 
long  before  this  point.  The  free  surface  disturbances  in  the  mildly  supercritical  flow  of 
the  water  tunnel  readily  introduced  a  noticeable  amount  of  air  into  the  body  of  the  jet.  It 
is,  therefore,  not  unreasonable  to  suspect  that  this  effect  would  be  even  more  dramatic  in 
the  case  of  the  rectangular  wall  jet  which  was  thinner,  flowed  at  a  higher  velocity  and 
produced  larger  disturbances.  Open  channel  flows  are  acceptably  well-described  by 
classic  boundary  layer  analyses  near  the  wall,  and  the  relatively  docile  free  surface 
variations  caused  by  slowly  rising  boils  and  kolks  do  not  compare  to  the  violent 
interactions  witnessed  here. 

Besides  the  noted  free  surface  eruptions,  the  DPIV  measurements  in  the  tunnel  jet 
point  out  some  additional  differences  between  wall  jets  and  other  boundary  layer  flows. 
Figure  60  shows  that  the  evolving  boundary  layer  of  the  finite  thickness  jet  (for  both  the 
smooth  and  rough  surface)  has  a  great  deal  of  influence  on  the  behavior  of  the  flow  in  the 
free  surface  zone.  In  ordinary  wall  jets  or  boundary  layers,  the  presence  of  the  wall  and 
the  influence  of  the  large  near-wall  velocity  gradient  diminishes  quickly  with  distance 
from  the  wall,  contradictory  to  those  of  the  finite  liquid  wall  jets  over  plane  and  curved 
surfaces  (as  in  the  present  study).  The  turbulent  profiles  of  Fig.  66  suggest  that  higher 
levels  of  turbulence  are  also  manifest  further  from  the  wall  in  agreement  with  the  value  of 
v'/Uo  calculated  using  Eq.  (9)  (compare  with  Figs.  4,  5  and  7).  It  is  important  to  point 
out  that  the  open  channel  data  of  Grass  (1971)  in  Fig.  61  is  at  nearly  the  same  level  of 
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relative  roughness  as  the  wall  jet,  but  it  is  also  a  subcritical  flow,  and  this  is  the  likely 
crux  of  the  difference.  As  the  flow  becomes  less  supercritical  (vis-a-vis  U  and  h),  the 
velocity  and  turbulence  of  the  inner  flow  begins  to  look  like  a  typical  boundary  layer  or 
submerged  wall  jet. 

What  distinguishes  the  bow  sheet  from  other  boundary  layer  flows  is  that  it  is 
highly  supercritical,  and  develops  eddies  with  sizes  on  the  order  of  h.  Additionally,  the 
inherent  fluidity  of  the  free  surface  allows  it  to  continuously  interact  with  the  body  of  the 
jet,  changing  the  entire  dynamics  of  the  internal  flow. 

4.  WaU  Curvature  and  Inclination 

Curvature  and  inclination  of  the  bow  of  a  ship  are  inescapable  realities  that  add 
another  dimension  to  the  bow  sheet  problem.  The  results  of  the  single  concave  curved- 
wall  jet  examined  here  show  that  curvature  reduces  the  density  of  the  spray,  as  well  as  the 
energy  of  the  drop-forming  eddies.  Part  of  the  cause  may  be  attributable  to  the  potential 
energy  sacrifice  required  for  the  flow  to  progress  up  the  wall,  but  it  also  seems  certain 
that  changes  in  internal  flow  characteristics  are  occurring  due  to  centrifugal  acceleration, 
as  suggested  by  the  smaller  mean  drop  size.  A  reduction  in  turbulence  intensity  (which  in 
the  wall  jet  appears  as  fewer  and  shorter  filaments)  has  been  seen  in  other  low  radius-of- 
curvature  wall  jet  flows  (Kobayashi  and  Fujisawa,  1983);  however,  based  on  the  previous 
discussion,  it  is  not  clear  if  the  same  results  can  be  applied  here.  Neither  inclination  nor 
curvature  had  a  noticeable  impact  on  the  location  of  the  free  surface  transitions,  so  it  is 
not  possible  to  ascertain  the  complete  nature  of  any  subsurface  differences  based  on  what 
was  observed.  However,  based  on  the  filament  and  drop  characteristics  alone,  it  is 
certain  that  curvature  has  attenuated  the  drop-forming  events  occurring  in  and  above  the 
free  surface  zone. 

5.  Numerical  Considerations 

The  nonisotropy  of  the  turbulence  and  the  near  impossibility  of  finding  a  suitable 
turbulent  model  that  incorporates  the  near  surface  hydro/aerodynamics  is  an 
insurmountable  task — particularly  when  the  free  surface  possesses  such  random  motion. 
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These  characteristics  also  make  it  unlikely  that  any  realistic  results  will  be  obtained  using 
anything  less  than  a  full  three-dimensional  model.  As  mentioned  in  the  introduction, 
there  are  specific  efforts  being  directed  in  this  area.  To  date,  some  CPD  models  have 
demonstrated  the  ability  to  model  flows  pouring  over  a  weir  and  the  "simple" 
axisymmetric  drop  pinch-off  problems.  It  is  still  undetermined  if  the  ability  to  model  all 
of  the  fluid  dynamic  components  exhibited  in  a  wall  jet  presently  exists.  Regardless,  it  is 
hoped  that  the  quantitative  data  obtained  here  will  provide  a  source  for  numerical 
verification,  and  perhaps  provide  some  direction  for  future  numerical  model-design 
consideration. 


6.  Application  to  Bow  Sheets 

This  project  never  lost  sight  of  its  ultimate  purpose.  By  examining  the 
characteristics  of  high  Froude,  Weber  and  Reynolds  number  jets,  we  wanted  to  simulate 
the  bow  sheet  formed  by  ships,  thereby  providing  a  means  of  predicting  the  factors  that 
contribute  to  spray  formation.  In  so  doing  we  have  identified  two  of  the  primary 
elements  which  govern  their  nature.  Specifically,  increasing  the  supercriticality  of  the 
flow  and  the  roughness  of  the  wall  leads  to  a  greater  manifestation  of  drops  and  spray. 
Hence,  the  very  first  suggestion  is  to  make  ships  smoother,  giving  the  greatest  attention  to 
the  leading  edge  of  the  bow. 

Some  consideration  should  also  be  given  to  the  shape  of  the  bow  (recognizing  that 
there  are  other  factors  which  also  influence  its  selection)  since  fixed  curvature  proved  to 
be  a  spray  inhibitor.  It  is  possible  that  more  complex  shape  considerations  such  as 
varying  the  rate  of  curvature,  or  using  subsurface  appendages  to  change  the  pressure 
distribution  near  the  bow,  could  result  in  thicker  bow  sheets  (or  no  bow  sheets,  as  in  the 
case  of  many  supertankers)  which  are  less  susceptible  to  spray  formation.  This  does  not 
preclude  the  possibility  that  smoothing  coupled  with  bow  shaping  might  result  in  even 
greater  reductions  in  spray  generation  and  all  of  its  attendant  consequences. 
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IV.  CONCLUSIONS 


The  following  conclusions  can  be  made  based  on  the  characteristics  that  have 
been  observed  and  measured  within  the  scope  of  this  investigation: 

•  Spray  generation  from  liquid  wall  jets  (bow  sheets)  requires  a  highly 
supercritical  (in  terms  of  U  and  h)  flow.  There  exists  some  minimum  value  of 
supercriticality,  below  which  drops  will  not  form,  the  local  plate  conditions 
are  below  this  level  (either  because  Uo  is  too  low  or  x/ho  is  sufficiently  large), 
drop-forming  filaments  will  not  develop. 

•  Wall  roughness  is  not  required  for  spray  generation,  but  it  significantly 
reduces  the  minimum  jet  velocity  necessary  for  this  event  to  transpire.  The 
greater  the  level  of  relative  roughness,  the  sooner  and  more  abundant  is  the 
resulting  spray.  Moreover,  one  need  only  reduce  the  roughness  below  some 
reasonable  finite  level  (perfect  smoothness  is  not  required),  to  substantially 
reduce  the  sensitivity  of  a  wall  jet  to  drop  production,  or  to  reduce  the 
frequency  of  the  drop  forming  events. 

•  The  filament  and  drop  production  process  is  of  a  statistical  nature. 

•  The  free  surface  provides  a  path  for  interaction  between  the  above  surface 
events  and  the  internal  flow  dynamics.  This  interaction  guarantees  that  the 
wall  jet  boundary  layer  is  fundamentally  different  from  that  of  typical 
turbulent  flows. 

•  Wall  curvature  appears  to  inhibit  the  formation  of  filaments  and  drops. 

•  Aerodynamics  do  not  influence  the  characteristic  behavior  of  the  filaments 
and  drops;  however,  free  surface  disturbances  may  act  to  introduce  air  bubbles 
into  the  body  of  the  jet. 
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V.  RECOMMENDATIONS 


The  following  research  is  recommended  as  food  for  thought: 

•  Make  quantitative  measurements  of  the  subsurface  flow  region  of  the  jet. 

•  Explore  the  mean  of  the  direct  measurement  of  the  local  wall  shear  stress. 

•  The  wall  jet  should  be  examined  in  other  facilities  that  will  allow  the 
Reynolds,  Weber  and  Froude  numbers  of  the  flow  to  be  uncoupled  so  that  the 
influence  of  inertia,  surface  tension  and  gravity  can  be  analyzed 
independently. 

•  An  effort  should  be  made  to  expand  the  "margins"  of  the  current  information 
database.  That  is,  determine  if  higher  levels  of  supercriticality  lead  to 
additional  phenomenological  transitions;  and,  ascertain  what  the  minimum 
levels  of  supercriticality  are  for  various  degrees  of  relative  roughness. 

•  Thoroughly  examine  the  effects  of  curvature  by  investigating  how  its 
influence  on  spray  formation  is  shaped  by  jet  velocity  and  wall  roughness. 
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APPENDIX  A.  FIGURES 


(1977a). 
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Figure  2.  Mean  velocity  distribution  across  the  turbulent  boundary  layer  formed  by 
smooth-  and  rough-walled  surfaces. 
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Figure  3.  Streamwise  turbulence  intensity  across  the  boundary  layer  formed  by  smooth 
and  rough-walled  surfaces. 
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Figure  4.  Vertical  turbulence  intensity  across  the  boundary  layer  formed  by  smooth- 
and  rough-walled  surfaces. 
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Figure  5.  Nondimensional  sweep  and  ejection  frequency  versus  boundary  layer 
thickness,  ti,  adjusted  for  roughness  element  height.  From  Raupach  (1981). 
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Figure  8.  CAD  drawing  of  the  rectangular  liquid  wall  jet  nozzle. 
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fluid  flow  path  is  as  indicated  by  the  arrows. 
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Figure  11.  Photograph  of  the  outside  of  the  rectangular  nozzle  with  the  cylindrical  gate 
open. 


70 


Figure  12.  Photograph  of  the  outside  of  the  rectangular  nozzle  with  the  cylindrical  gate 
closed. 
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Figure  13.  Photograph  of  the  rough  (k/ho  =  0.06)  curved  plate  used  in  the  rectangular  jet 
experiment. 
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Figure  14.  Photograph  of  a  rough  plate  installed  in  the  rectangular  liquid  wall  jet 
apparatus.  The  plate-mount  without  a  plate  attached  is  shown  in  the  upper  left-hand 
comer. 
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Figure  15.  CAD  drawing  of  the  axisymmetric  nozzle. 


74 


Figure  16.  Side  view  of  the  axisymmetric  nozzle  with  the  test  rod  assembly  installed. 
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Figure  17.  An  image  of  the  free  surface  profile  of  the  rectangular  liquid  wall  jet. 
Rex  =  2.3x10^  Re  =  3.6x10^  Fr  =  26,  We  =  3000,  k/ho  =  0.06. 


Figure  18.  An  image  of  a  three-dimensional  free  surface  distortion. 
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_ t  =  8 _ t=  18 _ 

Figure  19.  A  typical  sequence  of  images  captured  for  the  filament  and  drop  analysis. 
Rex  =  2.3x10^,  Re  =  3.6x1 0^  Fr  =  26,  We  =  3000,  k/ho  =  0.06. 


Figure  20.  Diagram  showing  how  the  filament  characteristics  are  simplified  for  statistical 
analysis. 
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c.  Drop  generation. 


Figure  21.  Photographs  of  the  different  regions  in  a  liquid  wall  jet. 
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Figure  24.  Normalized  mean  jet  thickness  of  the  liquid  wall  jet  at  the  beginning  of  the 
filament  formation  region.  Uncertainty  regions  are  bracketed  and  the  average  value  for 
the  drop-producing  runs  is  1.7. 


Figure  25.  Normalized  mean  jet  thickness  at  the  beginning  of  the  drop  generation 
region.  Uncertainty  regions  are  bracketed  and  the  average  value  is  2.0. 
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Figure  26.  A  laser  light  sheet  illuminated  photograph  of  the  free  surface  flow  over  the 
axisymmetric  liquid  wall  jet.  Row  is  from  right-to-left  and  the  actual  width  of  the  image 
is  about  four  centimeters. 
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Figure  29.  Distribution  of  the  normalized  drop  diameter  as  a  function  of  relative 
roughness. 
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Figure  30.  Mean  normalized  drop  diameter  as  a  function  of  the  relative  roughness. 
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Figure  31.  Distribution  of  the  normalized  drop  diameter  for  the  coarse-wall 
(k/ho  =  0. 13)  at  high-  and  low- Weber  number. 
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Figure  32.  Distribution  of  the  normalized  drop  diameter  for  the  flat-  and  curved-wall 
(k/ho  =  0.06,  We  =  3000). 
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Figure  33.  Distribution  of  the  normalized  filament  length  as  a  function  of  the  relative 
roughness. 
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Figure  34.  Distribution  of  the  normalized  filament  length  for  the  flat-  and  curved-wall 
(k/ho  =  0.06,  We  =  3000). 


93 


50 

40 

^  30 

D 
O 

O 

20 

10 

0 

Figure  36.  Distribution  of  the  normalized  filament  length  for  the  coarse-wall  (k/ho  = 
0.13)  at  high-  and  low-Weber  number. 
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k/ho  =  0.13,  We  =  3000 


k/ho  =  0.02,  We  =  3000 


Figure  37.  Photographs  of  three  different  runs  at  the  same  streamwise  location, 
x/ho  «  90. 
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the  relative  roughness. 
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Figure  39.  Distribution  of  the  normalized  effective  filament  diameter  for  the  coarse- 
wall  (k/ho  =  0. 13)  at  high-  and  low-Weber  number. 
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Figure  41.  Distribution  of  the  normalized  effective  filament  diameter  for  the  flat-  and 
curved-wall  (k/ho  =  0.06,  We  =  3000). 
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Figure  42.  Distribution  of  the  normalized  filament  slenderness  ratio  as  a  function  of 
the  relative  roughness. 
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Figure  44.  Distribution  of  the  normalized  filament  slenderness  ratio  for  the  coarse-wall 
(k/ho  =  0.13)  at  high-  and  low-Weber  number. 
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Figure  45.  Distribution  of  the  normalized  filament  slenderness  ratio  for  the  flat-  and 
curved-wall  (k/ho  =  0.06,  We  =  3000). 
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Figure  47.  Distribution  of  the  normalized  filament  volume  as  a  function  of  the  relative 
roughness. 


Figure  48.  Distribution  of  the  normalized  filament  volume  for  the  flat-  and  curved- 
wall  (k/ho  =  0.06,  We  =  3000). 
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Figure  49.  Distribution  of  the  normalized  filament  volume  for  the  coarse-wall  (k/ho 
0.13)  at  high-  and  low-Weber  number. 
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Figure  50.  Mean  value  of  the  drop  to  filament  volume  ratio  as  a  function  of  the  relative 
roughness. 
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Figure  51.  Mean  disturbance  wave  number  at  the  tip  of  the  filament  just  prior  to  drop 
separation. 
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Figure  52.  Distribution  of  the  disturbance  wave  number  as  a  function  of  relative 
roughness. 
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Figure  53.  Distribution  of  the  disturbance  wave  number  for  the  coarse-wall  (k/ho  = 
0.13)  at  high-  and  low-Weber  number. 
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(X  -  Xd)/ho 

Figure  55.  Normalized  mean  drop  diameter  as  a  function  of  the  nondimensional 
streamwise  coordinate.  Here  x  and  Xj  denote,  respectively,  the  distance  from  the 
leading  edge  of  the  plate  and  the  point  at  which  drop  generation  begins.  Only  those 
locations  with  >30  observations  are  shown.  The  uncertainty  is  approximately  10 
percent. 
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(X  -  Xd)/ho 

Figure  56.  Normalized  mean  filament  length  as  a  function  of  the  nondimensional 
streamwise  coordinate.  Here  x  and  Xj  denote,  respectively,  the  distance  from  the 
leading  edge  of  the  plate  and  the  point  at  which  drop  generation  begins.  Only  those 
locations  with  >30  observations  are  shown.  The  uncertainty  is  approximately  15 
percent. 
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(X  -  Xd)/ho 

Figure  57.  Normalized  mean  effective  filament  diameter  as  a  function  of  the 
nondimensional  streamwise  coordinate.  Here  x  and  Xj  denote,  respectively,  the 
distance  from  the  leading  edge  of  the  plate  and  the  point  at  which  drop  generation 
begins.  Only  those  locations  with  >30  observations  are  shown.  The  uncertainty  is 
approximately  10  percent. 
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Figure  58.  Normalized  mean  filament  slenderness  ratio  as  a  function  of  the 
nondimensional  streamwise  coordinate.  Here  x  and  X<j  denote,  respectively,  the 
distance  from  the  leading  edge  of  the  plate  and  the  point  at  which  drop  generation 
begins.  Only  those  locations  with  >30  observations  are  shown.  The  uncertainty  is 
approximately  15  percent. 
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(X  -  Xd)/ho 

Figure  59.  Normalized  mean  filament  volume  as  a  function  of  the  nondimensional 
streamwise  coordinate.  Here  x  and  denote,  respectively,  the  distance  from  the 
leading  edge  of  the  plate  and  the  point  at  which  drop  generation  begins.  Only  those 
locations  with  >30  observations  are  shown.  The  uncertainty  is  approximately  25 
percent. 
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Figure  60.  Distances  to  surface  transition  (ST),  filament  formation  (FF)  and  drop 
generation  (DG)  as  a  function  of  plate  roughness  for  the  rectangular  and  the  water 
tunnel  wall  jet. 
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Figure  61.  Mean  velocity  distribution  across  the  liquid  wall  jet  formed  in  the 
recirculating  free-surface  water  tunnel  at  x/ho  =  53. 
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Figure  62.  Photographs  of  the  interior  of  a  rough-wall  liquid  jet  formed  in  the 
recirculating  water  tunnel.  Rex  =  2.3x10®,  Re  =  2.8x10"^,  Fr  =  9,  We  =  1 100,  k/ho  =  0.08. 
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Figure  63.  Two  images  depicting  the  free  surface  distortions  and  bubbles  that  are  present 
in  the  rough-wall  jet  formed  in  the  recirculating  water  tunnel. 

Rex  =  1.6x10^  Re  =  2.8x10^  Fr  =  9,  We  =  1 100,  k/ho  =  0.08. 
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Figure  64.  A  sequence  of  images  illustrating  bubble  formation  below  the  free  surface  of 
the  rough-wall  jet  formed  in  the  recirculating  water  tunnel. 

Rex  =  1.6x10^  Re  =  2.8x10^  Fr  =  9,  We  =  1 100,  k/ho  =  0.08. 
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Figure  65.  A  phoiograpn  ox  a  iiiamem  aoove  ana  me  Duooies  present  oelow  the  free 
surface  of  the  rough-wall  jet  formed  in  the  recirculating  water  tunnel. 

Rex  =  2.7x10^  Re  =  2.8x10^  Fr  =  9,  We  =  1 100,  kdio  =  0.08. 
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Figure  66.  Streamwise  and  vertical  turbulence  intensity  across  the  liquid  wall  jet 
formed  in  the  recirculating  water  tunnel  at  x/ho  =  53. 


Figure  67.  Contour  plots  of  v-velocity  (at  three  different  times)  in  the  rough-wall  jet 
formed  in  the  recirculating  water  tunnel  at  x/ho  =  53,  as  acquired  with  DPIV.  The  white 
shade  indicates  v  <  -0.5  m/s,  while  the  black  is  for  v  >  0.5  m/s.  The  circled  regions  show 
"blobs"  with  large  upward  velocities. 
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Figure  69.  The  jet  Weber  number  versus  the  jet  Froude  number  for  different  levels  of 
v'/Uo  as  estimated  by  Eq.  (10).  The  circles  mark  the  three  combinations  of  We-Fr 
examined  in  this  investigation. 
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APPENDIX  B.  UNCERTAINTY  ANALYSIS 


The  uncertainty  of  the  data  presented  herein  was  determined  at  a  95  percent 
confidence  interval.  The  biases  associated  with  the  measurements  appeared 
inconsequential  relative  to  the  random  variance  of  the  data  and  the  precision  level  of  the 
individual  measurements.  Therefore,  the  data  variance  and  the  number  of  observations 
dictate  the  mean-value  uncertainty,  while  the  individual  measurement  uncertainties  are 
precision  controlled.  To  minimize  the  uncertainty,  it  is  important  to  record  as  many 
observations  as  reasonably  possible,  and  to  arrange  the  measuring  equipment  to 
maximize  its  accuracy.  Since  this  is  the  first  investigation  into  the  physics  of  this  type  of 
drop  formation  process,  there  is  no  a  priori  knowledge  of  the  various  parameter 
distributions.  Hence,  the  number  of  observations  of  a  given  feature  is  >  30  whenever 
practical,  so  that  the  implicit  analytical  assumption  of  a  Gaussian  distribution  does  not 
introduce  significant  error  into  the  mean-value  confidence  intervals  (Kline  &  McClintock 
1953). 

The  surface  characteristics  (i.e.,  liquid  filaments  and  drops)  of  the  rectangular 
liquid  wall  jet  are  recorded  on  videotape  using  a  high-speed  camera  (imager).  The 
videotape  is  then  reviewed  to  find  filaments  and  drops  that  can  be  analyzed  with  a 
reasonable  degree  of  accuracy.  The  criteria  used  to  make  this  determination  are: 

•  The  base  of  the  filament  is  clearly  discemable. 

•  The  filament  behaves  in  an  axisymmetric  manner  (as  witnessed  by  its 
behavior  in  numerous  images)  and  has  a  slenderness  ratio  of  about  1.5  to  2. 

•  The  entire  filament,  and  the  drops  it  produces,  are  unobstructed  by  other 
surface  stmctures. 

•  The  filament  and  drops  are  in  focus. 

•  The  drops  are  nearly  circular  in  images  used  to  measure  their  areas. 

•  The  drops  and  filaments  can  be  measured  in  at  least  three  sequential  images. 
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The  images  of  filaments  and  drops  meeting  these  criteria  are  digitized,  and  then 
analyzed  on  a  PC  using  the  Optimas  MA  software.  Optimas  allows  the  user  to  mark 
points  and  lines,  and  to  trace  out  the  areas  of  the  objects  seen  in  the  digitized  image.  It  is 
able  to  track  these  features  throughout  a  sequence  of  images,  and  it  allows  the  user  to 
redefine  the  objects  in  each  image.  Based  on  the  length  scale  prescribed  and  the  time 
between  the  images,  Optimas  calculates  the  two-dimensional  areas,  angles,  lengths, 
centroids,  and  relative  positions,  as  well  as  velocities  and  accelerations  of  the  various 
items.  After  the  image  analysis  is  complete,  Optimas  writes  the  selected  values  to  a  text 
file  or  a  Microsoft  Excel  spreadsheet. 

Using  a  high-quality  Fuji  105mm/fl.8  lens,  and  placing  the  camera  at  an 
appropriate  position,  the  center  of  the  wall  jet  can  be  brought  near  the  minimum  focal 
distance  with  the  lens  f-stop  at  f  1.8.  This  maintains  a  consistently  small  depth  of  field  for 
the  given  lighting,  thereby  ensuring  that  only  those  filaments  and  drops  formed  in  the 
center  of  the  liquid  jet  are  subsequently  analyzed.  This  focal  distance  is  short  enough  to 
make  these  filaments  and  drops  appear  sufficiently  large  for  accurate  evaluation,  while 
ensuring  that  the  camera  far  enough  back  to  capture  a  sequence  of  images  of  a  given 
feature  (typically  8  to  15)  without  having  to  move  it. 

The  uncertainty  analysis  performed  for  this  study  uses  the  constant  odds  (root- 
mean-square)  formula,  and  95  percent  confidence  level  uncertainty  intervals  as 
recommended  by  Moffat  (1982).  Adding  a  bit  of  complexity  to  the  calculations  is  the 
fact  that  some  of  the  parameters  requiring  uncertainty  analysis  are  not  measured  directly. 
For  example,  the  calculation  of  the  drop  diameter  is  based  on  its  measured  cross-sectional 
area.  Nonetheless,  the  primary  uncertainties  associated  with  this  analysis  are  a  function 
of  the  overall  imaging  system  characteristics,  and  the  user’s  ability  to  accurately  identify 
the  free  surface  features  using  the  Optimas  software.  For  the  items  measured  directly 
with  Optimas,  the  identifiable  uncertainties  and  their  estimated  values  are  listed  below: 

•  Object  definition  (OBJ).  This  describes  the  uncertainty  associated  with  the 
user's  ability  to  outline  the  exact  shape  of  an  object.  For  all  objects  measured, 
there  is  some  blurring  around  the  edges  because  of  object  motion  and  limited 
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digital  resolution.  Hence,  this  uncertainty  is  not  a  function  of  the  size  of  the 
object,  but  is  instead  a  fixed  value  applied  to  all  area  objects.  An  analysis  of  a 
clear  glass  sphere  and  a  small  acrylic  rod  moving  through  the  range  of 
velocities  considered  in  this  study  show  0.08  mm  <  6obj  ^  0.20  mm.  For  the 
following  uncertainty  analysis  Sqbj  =  0.14  mm  will  be  used. 

•  Object  not  at  the  focal  plane  (FOCUS).  This  describes  the  uncertainty 
introduced  when  objects  measured  are  not  at  the  precise  focal  plane  of  the  lens 
(not  at  the  center  of  the  wall  jet),  but  still  appear  in  focus  (within  the  optical 
depth  of  field).  By  considering  objects  placed  both  closer  and  further  from  the 
lens,  this  uncertainty  is  estimated  as  5focus  =  4%. 

•  Lens  distortion  (LENS).  Because  imperfections  will  exist  in  the  optics  of  any 
lens,  there  will  be  some  amount  of  distortion — most  likely  near  the  edges  of 
the  image.  However,  for  the  lens  used  in  this  experiment,  the  distortion  is 
negligible:  5lens  ~0%. 

•  Length  scale  calibration  error  (LENGTH).  Optimas  bases  its  dimensional 
calculations  on  the  user's  length-scale  definition.  By  using  Optimas  to  draw  a 
line  over  a  known  distance  on  the  scale,  say  10  mm,  and  comparing  the 
Optimas-calculated  length  with  the  actual  length,  this  uncertainty  is  estimated 
as  8length  =  0.9  %. 

•  Optimas  area  calculation  approximation  (AREA).  Optimas  determines  the 
area  of  an  outlined  shape  by  adding  up  the  pixels  contained  within  the  defined 
region.  If  an  outline  passes  through  a  pixel,  then  the  entire  pixel  may 
contribute  to  the  estimated  object  area  if  a  large  enough  fraction  of  it  is  within 
the  outline.  Therefore,  the  same  shape  at  a  different  location  in  the  image  may 
have  a  slightly  different  area.  By  outlining  three  different  size  shapes,  placing 
them  in  various  locations  on  different  images,  and  having  Optimas  calculate 
the  area  at  each  location,  this  uncertainty  is  estimated  to  be  6area  =  0.3  %. 

•  Observer  interpretation  (OBS).  Different  individuals  will  define  an  object 
somewhat  differently.  The  uncertainty  of  this  factor  is  determined  for  each 
measurement  by  having  two  individuals  independently  analyze  a  number  of 
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identical  image  sequences.  As  shown  in  the  table  below,  the  result  is  a 
function  of  what  is  measured,  and  is  generally  the  largest  contributor  to  the 
precision  uncertainty. 


Filament  Drop  Filament  Filament  Filament  Point 

Area  Area  Length  Angle  Disturbance  (mm) 

0.15  I  0.07  I  0.06  I  0.05  |  0.18  | 

Table  5.  Uncertainty  of  observer  interpretation,  Sobs. 

•  Scale  error  (SCALE).  This  is  the  uncertainty  associated  with  the  inaccuracy 
of  the  rule  used  to  define  the  length  scale  in  Optimas.  By  comparing  the  one 
used  with  that  of  a  six-inch  machinist  rule  of  known  accuracy  (+  0.102  mm,  - 
0.051  mm),  the  scale  error  is  estimated  at  Sscale  <  0.2  mm.  However,  this  is 
only  applied  when  determining  point  position  uncertainty.  In  the  case  of  areas 
and  filament  lengths,  this  error  is  part  of  Sobj- 


As  mentioned  previously,  the  drop  diameter  is  calculated  based  on  the  cross- 
sectional  area  of  the  drop  in  the  image.  For  each  drop,  its  "final"  diameter  is  obtained  by 
averaging  the  values  obtained  for  the  drop  in  each  image  for  which  it  appears.  Because 
the  size  of  the  drop  is  fairly  constant  fi-om  frame-to-fi:ame,  it  is  assumed  that  the 
uncertainty  of  the  drop  diameter  in  each  frame  is  the  same  giving 


SD^_  1  2r5A^^ 

D^  Tn’k^D^  f 


(A.1) 


where  N  is  the  number  of  images  in  which  the  drop  area  is  measured.  This  can  be  as  few 
as  three  or  as  many  as  ten,  but  a  typical  number  of  images  is  five.  Using  this  value  for  N, 
the  relative  uncertainty  for  the  range  of  drop  sizes  considered  in  this  study  is 


0.03  <^^<0.10. 

D., 


It  is  found  that  bDdi/Ddr «  0.04  for  the  majority  of  the  drops. 


Similarly,  the  relative  uncertainty  of  the  effective  filament  diameter  is  estimated 
from  its  cross-sectional  area  assuming,  for  the  purpose  of  the  uncertainty  calculation,  that 
Af  =  Df  Lf.  This  results  in  the  following  relationship: 


SDf  _ 
Df  “ 


A, 


V 


Uf 


1/2 


(A.2) 


This  equation  requires  the  relative  uncertainty  of  the  filament  length  in  addition  to  that  of 
the  filament  area.  Using  four  of  the  measurement  values  cited  earlier  (8obj,  5focus, 
Slength  and  5obs)  the  length  uncertainty  is 

6Lf 


Lf 


«  0.075  for  all  filaments. 


This  then  leads  to  a  filament  diameter  uncertainty  of 

6D, 

0.17< — ^-<0.21. 
Df 


As  expected,  this  particular  diameter  has  a  higher  level  of  uncertainty  than  that  of  the 
drop.  This  outcome  is  attributed  to  the  larger  number  of  uncertainties  involved;  the 
approximation  of  defining  the  filament  area  as  the  sum  of  its  length  and  mean  diameter; 
and,  the  fact  that  the  filament  diameter  is  not  the  average  of  multiple  measurements,  as  is 
the  drop. 

The  relative  uncertainty  of  the  filament  slenderness  ratio 


and  volume 

0.35  <^<0.43 
Vf 

are  functions  of  the  filament  length  and  diameter  uncertainties  given  above. 
Additionally,  the  estimate  of  the  filament  volume  uncertainty  assumes  that  Vf  =  jrD/l^4. 

Neither  the  determination  of  filament  angle  nor  filament  disturbance  length  lends 
itself  to  a  full  uncertainty  analysis.  Since  both  are  simple  user  interpretations — and  are 
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not  used  in  any  analytical  analysis — it  is  appropriate  to  only  consider  the  observer 
interpretation  uncertainties  of  these  values.  As  such,  6tx/a  «  0.06  and  6ki/ki «  0. 18. 

The  wall  jet  analysis  also  measures  the  U-  and  V-velocity  of  the  drop  centroid, 
and  the  points  marking  the  filament  base  and  tip.  Optimas  calculates  velocity  by 
comparing  the  change  in  spatial  position  of  an  object  in  a  pair  of  sequential  images  (U  = 
AX/AT,  V  =  AY/AT),  resulting  in  a  basic  uncertainty  equation  of 

I  1/2 

(A.3) 


U 


^6(AX)Y/5(AT)^ 


AX 


J 


AT 
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Therefore,  the  uncertainty  of  this  calculation  is  a  function  of  the  spatial  precision 
associated  with  an  object,  and  the  timing  precision  of  the  NAC  HSV-500  imaging 
system.  According  to  the  manufacturer,  the  imaging  system  precision  is  better  than  lO"^ 
ms.  In  the  case  of  the  drop  velocity  estimation,  AX  and  AY  measure  the  change  in 
position  of  the  drop  centroid  from  one  image  to  the  next.  Therefore,  the  spatial  resolution 
is  a  function  of  how  well  the  area  of  the  drop  is  defined  in  each  image: 
(AX,AY)=/[Adr(t  =  l),Aj,(t  =  2)].  Using  the  measurement  uncertainties  associated 


with  the  drop  area  calculations,  and  the  given  timing  accuracy  for  a  fixed  AT  of  0.002  s, 
the  drop  U-velocity  uncertainty  is 

0.06  <^^<0.25. 

Ua, 

By  virtue  of  their  cross-sectional  area,  over  90  percent  of  the  drops  identified  have 
relative  U-velocity  uncertainties  of  <  0.10.  The  drop  V-velocity  component  has  an 
uncertainty  range  of 

SV 

0.13  < — ^<0.55. 

V* 

For  the  velocity  of  the  filament  base  and  tip  points  we  must  now  consider  scale  error, 
SscALE,  and  the  observer  interpretation  uncertainty,  6obs,  associated  with  a  point  in  both 
images.  For  the  range  of  U-velocities  witnessed  in  this  study  (4  to  8  m/s),  the  point 
uncertainty  is 
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0.026  S 


SU^i, 

^  point 


<0.04. 


On  the  other  hand,  the  uncertainty  in  the  V-velocity  component  has  a  much  wider  range 
since  the  velocity  magnitude  varies  from  zero  to  four  meters  per  second.  Although  many 
near-zero  vertical  point  velocities  are  recorded,  the  uncertainty  is  not  reduced  to  a 
reasonable  level  (say  SVpomtA^poim  *  0.35)  unless  the  velocity  exceeds  one  meter  per 
second.  Additionally,  the  U-velocity  uncertainties  are  lower  for  both  drops  and  points 
because  they  are  averaged  over  several  images.  Doing  so  reduces  the  uncertainty  by 


i/Vn  ,  where  N  is  the  number  of  image  pairs. 

Now  let  us  consider  the  uncertainty  of  the  jet  velocity  at  the  exit  of  the  rectangular 
nozzle,  Uo.  Defining  this  velocity  by  the  discharge  coefficient  of  the  nozzle  and  the 
hydrostatic  head  in  the  U-tunnel,  results  in  a  relative  uncertainty  equation  of 


8Uo 

Uo 


8C. 


1/2 


(A.4) 


Substituting  in  the  maximum  variation  of  the  discharge  coefficient  (8Cd  =  0.03),  the  mean 
value  of  the  discharge  coefficient  (Ca  =  0.91),  the  U-tunnel  levels  during  experimentation 
1.2  m),  and  the  maximum  error  of  the  manometer  (8^  =  25  mm)  into  Eqn.  (A.4) 
gives: 


8Uo 

Uo 


0.035. 


Combining  this  with  the  result  of  Eqn.  (A.3)  allows  one  to  determine  the  uncertainty  of 
any  velocity  nondimensionalized  by  Uo. 

The  next  uncertainty  to  consider  is  that  associated  with  the  jet  thickness 
measurement.  The  uncertainties  of  this  result  include  Optimas  length  scale  calibration, 
Slength,  scale  error,  8scale,  scale  interpolation,  8intp,  and  the  jet  thickness  variation 
resulting  from  measuring  the  thickness  over  a  rough  surface,  8rough-  The  estimate  of  the 
first  three  uncertainties  is  0.9  percent,  <  0.25  mm  and  <  0.05  mm,  respectively.  The 
surface  roughness  uncertainty  is  estimated  at  two-thirds  of  the  roughness  element  size 
variation  shown  in  Table  1  (Chapter  2).  Assuming  a  normally  distributed  element-size 
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distribution,  this  roughly  equates  to  two  standard  deviations.  Table  A2  shows  the  values 
of  roughness  uncertainty  for  the  four  surfaces  examined. 


Smooth 

Fine 

Roughness 

Medium 

Roughness 

Coarse 

Roughness 

5rough  (mm) 

0.00 

0.014 

0.018 

0.083 

5h/h 

0.046 

0.046 

0.046 

0.048 

Table  6.  Jet  thickness  uncertainty  for  each  level  of  plate  roughness. 


The  uncertainty  equation  for  jet  thickness  is 


h 


0.07  +  6' 


ROUGH 


+  7.92x10"^ 


(A.5) 


with  the  solutions  to  Eqn.  (A.5)  are  also  provided  in  Table  A2  (second  row). 

Finally,  one  must  examine  the  uncertainty  of  the  mean  jet  thickness,  h ,  at  values 
of  x/ho  for  which  the  thickness  was  not  measured.  For  example,  the  jet  thickness  at  the 
points  marking  the  beginning  of  the  three  surface  transitions,  Xtr/ho,  Xp'ho  and  Xj/ho. 
The  uncertainty  of  h  at  these  locations  is  a  function  of  the  mean  jet  thickness  at  the 
locations  where  h  is  measured,  and  the  estimated  streamwise  location  of  the  surface 
transition  point.  From  videotape  analysis,  the  x/ho  values  have  been  identified  to  within 
±  30  mm.  h^ ,  the  mean  jet  thickness  at  the  transition  locations,  is  found  by  linearly 


interpolating  between  the  known  values  of  h.  Differentiating  this  line  with  respect  to  x 
gives  the  slope  of  the  line,  m,  which  is  then  used  in  the  uncertainty  estimate: 

l>/2 


Shaai. 

ho 


5h 

ho 


+ 


m- 


ho 


(A.6) 


Based  on  a  mean  jet  thickness  of  5.7  mm,  5X/ho  is  estimated  at  4.5  for  all  three  transition 
locations.  The  first  term  in  Eqn.  (A.6)  is  found  by  linearly  interpolating  between  the 
confidence  intervals  for  the  known  values  of  h . 

The  confidence  intervals  for  all  of  the  mean-value  parameters  presented  in  this 
study  are  evaluated  by  employing  commonly  used  single-sample  analysis  techniques 
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(Coleman  &  Steele  1989).  That  is,  an  estimate  of  the  population  standard  deviation  is 
determined  using  the  unbiased  estimator, 


Sx  = 


n‘/2 


.N-IS 

The  sample  mean  standard  deviation  is  then  approximated  by 

Ss=sJVn. 

Finally,  this  value  is  multiplied  by  2.0  to  establish  the  95  percent  confidence  level  of  the 
sample  mean  value.  Unless  otherwise  stated,  the  number  of  observations  in  the  sample  is 


>30. 
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APPENDIX  C.  TABULATION  OF  DATA 


This  appendix  contains  a  listing  of  the  data  discussed  in  the  thesis.  The  first  three 
pages  are  a  summary  of  the  statistical  parameters  associated  with  the  measurements. 
They  provide  the  group  parameters  describing  the  statistics  of  the  data  when  arranged  in  a 
histogram  format.  The  sample  statistics  shown  are  based  on  the  individual  observations 
for  each  mn.  The  remaining  pages  list  the  individual  filament  and  drop  measurements 
(arranged  by  run)  as  described  in  the  writing. 
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k/ho 

lile 

Lf/ho 

Q.OO 

O.OE 

O.Ob 

0.13 

0.13 

0.0k 

4700 

3000 

3000 

3000 

IkOO 

3000 

Flat 

Flat 

Flat 

Flat 

Flat 

Curved 

Group  [lean 

1.37 

1.45 

1.^4 

a. 33 

l.kT 

1.5S 

Group  Variance 

0.14 

0.34 

0.47 

O.'IE 

0.3T 

0.3k 

Group  Std  Dev 

0.37 

0.5fl 

D.bT 

O.Tk 

0.k3 

O.kO 

Group  Skewness 

0.30 

0.7‘1 

0.3k 

1.17 

1.45 

0.77 

Group  Kurtosis 

3.53 

3. IS 

E.55 

4.k0 

k.5k 

3.13 

Sample  Mean 

1.37 

1.45 

1.^3 

S.33 

1.70 

1.55 

SampleVariance 

0.15 

0.30 

0.44 

O.TO 

0.3k 

0.34 

Sample  Std  Dev 

0.34 

0.54 

0  .  kk 

0.T5 

O.kO 

0.5a 

Sample  Skewness 

0.fi5 

0.A5 

0.5S 

l.Ek 

1.5a 

0-^4 

Sample  Kurtosis 

l.Ofi 

0.10 

-O.Ek 

5.0A 

3.4a 

a.3a 

Sample  Mode 

l.Et 

0.11 

1.31 

l.flT 

1.5S 

i.Ea 

Sample  Median 

1.31 

1.3E 

i.a3 

S.ll 

1.5a 

1.40 

k/ho 

lile 

Ddr/ho 

0.00 

O.OE 

0.0k 

0.13 

0.13 

0.0k 

4700 

3000 

3000 

3000 

IkOO 

3000 

Flat 

Flat 

Flat 

Flat 

Flat 

Curved 

Group  Mean 

0.510 

0.550 

0.51k 

0.kE5 

O.kOk 

0.41k 

Group  Variance 

o.ooa 

0.015 

0.014 

0.011 

D.0E3 

0.017 

Group  Std  Dev 

O.O'IS 

O.IEE 

O.IED 

0.137 

D.15E 

D.1E1 

Group  Skewness 

-0.3a4 

0.E11 

D.Eia 

Q.7aa 

-o.isa 

a.4k3 

Group  Kurtosis 

3.773 

S.513 

S.17D 

3.357 

E.S51 

3.104 

Sample  Mean 

0.51k 

0.5k0 

O.kOO 

Q.kE4 

O.klO 

0.50E 

SampleVariance 

0.00=i 

0.014 

0.014 

0.011 

D.0E4 

0.015 

Sample  Std  Dev 

0.013 

O-lEO 

0.117 

0.131 

0.154 

D.1E3 

Sample  Skewness 

-0.345 

0.345 

D.EIk 

D.7as 

-0.1k7 

0.55k 

Sample  Kurtosis 

1.1E3 

-0.4S1 

0.011 

0.354 

-D.ai4 

0.751 

Sample  Mode 

0.540 

0.410 

D.5kD 

0.413 

0.k50 

0.540 

Sample  Median 

0.5S0 

o.5ka 

D.51D 

O.kOE 

0.kE5 

0.500 
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k/ho 

Ue 
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Uall  type:  Roughi  Flat 


Lf/ho  Df/ho  angle(deg)  ki  Ddr/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udpop/Uo 
l.MM  0.31=  flfl.T  a.D4  D.4T  □.Qt.Q  Q.D41  D-TS  0.12  Q.'ia 
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Lf/ho  Df/ho  angle(deg)  ki  Dd,./ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
□  .=m  D-BM  t.?.'!  l.'lb  Q.3b  O.D3fl  D-DS?  Q.fl?  D.fl6  D.flS 


Lf/ho  Df/ho  angle(deg)  ki  Ddp/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
1.S3  0.30  me-l  1.73  Q.hQ  -D.017  -Q.D2b  Q.fl^  O.flD  Q.flb 
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Lf/ho  Df/ho  angle(deg)  ki  Ddr/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
B.13  Q.St  titi-t.  B-13  O.flB  -Q.QBfl  0.77  D.flb 

B.D7  0.3S  flT.S  l.SS  O-bl  O.DDt.  D.DBb  D.74  Q.77  □.7b 

□  .Tl  Q.3M  Tb.b  B.3B  D.M3  D.D3D  D-OBb  0.7fl  D.77  0.76 


Iilall  type:  Roughi  Flat 
Uo  =  b  •  2  iti/s 
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k  =0-33  mm 
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DiD-D  2iD 


□  It.  0.003 
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Uall  type:  Roughi  Flat 


Lf/ho  Df/ho  angle(deg)  kj  Ddr/ho  Vtip/Uo  Vo/Uo  Upase/Uo  Utip/Uo  Udrop/Uo 
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Lf/ho  Df/ho  angle(deg)  ki  Ddr/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
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Lf/ho  Df/ho  angle(deg)  ki  Dd^/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
5.3^  Q.m^  51-1  l.fll  D.7b  Q.DOM  Q.QIS  0-77  0.77  D.7M 


Iilall  type:  Rough-i  Flat 
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tr‘tHaab-'bJab-*b-'IHb-'lHtHb-‘b-*njtHb-'b-«tHnJbT‘b-«b-'b-*nJtr'njnJb-'tHlHfcH 

LnLna]O3nJD3«Il^jrtnrULjUrU03br**i}[r’Ln«JJlb-'OaJl-JLnfc-'UJ-C’rLnII'D3nj 

Oaab-’tn-Ufcr^JOaJJ-xlUIlJJCnCnrUaLnrUUJUJDalPOat-'b-'LnjrUJtr'rUUJtnnJ 


aaaaaaaaaoaaaaaaaaaaoaoaaaaaaaaaa 

ujLurum.cr-NjmmijJLU-cujrujrru-crrux'-c'LJLijrLn-irjrLntnnj-c'X’LunjuJnj 
ODT  wDLnbH^oaruujujriJUjjaoacncrtHrFtn-jatjiLnuj-cx'tr'-vJUjnjruLn 


fc-i  tJ  bH  tr*  fc-»  b-* 

b-'aaoaJcroRj-NjjrjjoLncr-j-Njoatn-ojfcr'jr-jcrDaujnrDajrLn-j^j] 
Lua^cr  Ji^tno-Jb-'b-*03bH03^jjr-ox’Lna-oruj3[ro3Lnruirj3oruja 

nj03jrDrjrbH-jruirir*flr-NJUjj[3ruruoj]cnUjruinD3aujDa-j[rcj*-j-JD3D3 


b-»t-'mmb-‘UJrUb-*fc^b-'b-'rUb-‘b-‘bJb-'rUfc-'rururiJb-'fc-*b-‘b-'b-‘fc-'h-*rUIHb-‘b-'b-* 

ir*-jru-crj]b-*En-X]OLn£FLnLnLn[rcrcrcranjacr-jLncr^uj(parLi->j-ncr 

03LnCr*b-*^-0fc-'J]jracn4rarUb-*JlQ3UJLLianJLnLn-JrU-crLnDaO3JD’-JLnO3 


oaoaaaaaaaaaaaaaooaaaaaaaaaaaoaaa 

(JiCPUJUJLn[nLLlUJLnCP-0LiJLn-JLn^ru-sJCrCnLn*«NJ-J-JJlOa-n-J[r'CFLLl[n-C‘ 

rub-'OaiHJiauj-c’tH^uj-j-c’rua-irDaaLnb-'DaDatrnjacr-jfc-'ruir-jtntn 


■  I  I  I  I  I  t  I  i  I  I  I  II 

aaaaoaaaaaaaaaoaaaaaaaaaaaaoaaaaa 


O  O  O  C3  O  O  far*  CD  O  O  CD  O  O  O  O  LxJ  C3  CD  O  CD  OOOCDCDOOO  O  b-*  CD  CD 

b-*ujjicDb-‘CFru-c-fcr'uj-\jaijjcraruao3crtn-vj-c’-c‘tjiLLjruabT*D-cijjir* 

JJLnacDjacr-jnjorruj]tnrujin-jciLnujLLiJi-druijjru-OLnb-*LnJ]-oji 


I  III  I  I  I  I  I  I  II 

□  acDaaatDoatDaaaaaacDaaacDCDaacDCDaaaaacD 


CDaacDaab-*aaacDCDaacDriJCDCDacDCDaaaatDaaab-'CDa 

b-*UJ-Uab-*-C‘b-*jrrUb-*-%]Db-'Lna-JJCDCrJ3UJUJLjLJLJtnCn-Ca[Db^UJb-*-0 

crruirruoaDaCDjosb-'acjaLnirjrjruaajrruorm-cauiruLn-oruirj 


cDDaaoaaoaacDacDaacDaacDCDDCDaaaaaoaaaacD 

J]rUUJCDOLnUJOaxrLLirUb-^LnnJ030lJjaCJ>Cn-vlCnUJaLnaC^LnarUJ30tn 


aacDaacDCDaaaaacDaaaaaaaacDaatDaaaaaaacD 

-Jtr-xJ-O-O-U-UO-ODaOa-O-JCti-JttJ-J-J-J-U-UOa-J-J-JJlCrOaCJaDa-JOaCb* 

UJjD-c:j:,j3-%jnjdLnLiJ-c’ruiriiJJ3Dnjin-o-o-u^[rb-'crtDatn(r"j]Cjacj3cr 


CDoaoaacDCDaaaacDaaaaacDCDCDDCDCDDacDacDCDacDa 


-Nj-J-sJ-U-U-O-U-D-NJ-JOa-O-J-d-J-O-U-J-O-J-U-O-OCr-J-O-vIOaCJaDa-JC^CJa 

•JJ-C'CnJJOs-UUJarLILlJnJLnUJJDOa-JUJb-'jr-aDllCnJCrJlCFtnOa-JQ^EFDaCF 
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Lf/ho  Df/ho  angle(deg)  ki  Ddp/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udpop/Uo 
5.^fl  0-35  7^.3  l.a^  D.7D  0-053  0-017  D-flS  0-fl3  0-7^ 


Lf/ho  Df/ho  angleCdeg)  ki  Ddr/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udpop/Uo 
l-ll  0.36  110. a  2. SO  0.S^  -0.0M7  -0-031.  O-Tb  O-Bl  0.73 


I 

! 


jjmmarur-^mr^crr^r^JiErmcaar^j-muiacoLntna-rucrJ-ni-neOi-Rzrm 


aaoaaaaoaoDaacJCDaaaaaDaaaaaaaaoaaa 


jjiJ3amr^p-aruEramruLnr^Hinjr-uir^ajai-=ii-nmruireacojiJioir 

cajr^car^r^-Dr'-JiJ3r^r-jr^r^r^r“-Dr^r^r^-Dr^r^r^-D-D-D-D-D-D‘=or^ 


aaoaaaaaaadooQoaooDOoaaaaaaaaaaaa 


jLni-qcoruaEar^Jr^amj^muiEra'mar^Lnrnnj^rqjiLnzruiLnn-m 

p-j]e0p“r-p~j3j-j3r-r-p“*-nr“r-p-jr-r-r-jijir-r>-r-j]jij3jaj3jap-r^ 


oooaoooooocDOOoaaooooaaooooooaoooo 


ruacrnjr-ajorRr^r^mLnLjnrrcoJr^majr-ruj3H=ir-coj3Jii-nrutJnrqr- 

orujoHia-mJ-aaDCDrururur-amaLntnrrunjruuntnatncjaa 

oaaoaar^aoaaoDoaocDr^onaDaaaaaoaaaa 


oooooooooaaooocsoooooooooooaooooo 

I  It  I  i  III  t  III  II 

jcraa-r^arur^^Ln«orumi-RircrLn.3“j)j’aar^r-.3-r^mr^r“r-aia 

Dmr-c3ru«orurur=ia.rrrn»^ru-ii*jiai-ni-R-j]ec3crnjruar“jmjinjr^nj 

aoaaoor^aoooaaaaaoiHDaaaoaoaaaaoaa 


aaaaaaaaaaaaaaaaaaaaoaaaaaaaDoaaa 

I  I  I  I  I  I  I  I  III  I  I 


njJii-=?Jir-r“r^r=ir^n-*-nar“r=ii-nj.3"crtrr=iaea[rmErrup-iT'tnunj-rneo 

cai-nj-i-njiaia“r-JiLnLn«ocaLnr^ear^rur^^jjmr-J3J3^j"j3J-j33'iJiLn 


oaOOOOOOOOOOOOCDOOOOOOOOCIIOOOOOOOCDO 


r=ir-r=i3‘j]cazraamr^^r^«t]mnjru.3'uirqcDmr^njnjHitn[rnjH?p-njm 

i-qirmJir^i-nDrur-r-ap-r-oorr-oar-iFr-i-n^jjcFJ-i-njjajia-tFD 


rur^rui^Rjr^rijriJr^rRrur^r^r^r^i^njrnr^r^r^rur-^r-^i-^rurur^rur^njnjnj 


r-aimLrj.3‘P-coiJiir 


[ra-mp-mjir^aj] 
uimnjcaeaeOr^Bra 
i-R  1-^  rR  rR 


rrai-nrrr-aLnir 


CDnjp-mcaaLna 

crDeo.3-eor-crea 


jotorumji-Rj 


a-Eomn-rRruLDj 

njr-rup-crcr*or^ 

rR  I-R 


rumarurRJiru-D 


mcor^mrurujoi 
aoEor-anjjop- 
rR  r-R  i-R  rR 


mcOrRJli-R^ea3“m 

Lnrummi-nmrua^a" 


trrumjDETrRP- 
m  rr  Ln  m  LH  Lxi  tn 


comjcrmBor-njjr-azrcarRi-nji 

mirrumma“mj*mma-mrRiruim 


aaaaaaaaaaaaaaDaaoaanaoaaaaaaaaaa 


p-crrRaeaai-nj3Ji-njjjrR.3"i-Rirjar^rRarumcrjatr)Lnc^j[]LnrReaj] 

Drrar^zrcomaJirurrcumi-RLnLnr-eocarujcreanjmrrr-rRr-rRrur-a 


rU  rR  rR  rR  rR  RJ  rR  rR  rR  rR  rR  fU  HJ  rR  rR  rR  rR  O  RJ  rR  m  rR  rR  rR  rR  rR  rR  rR  rR  rR  rR  rR  rR 
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•lOO-O  hOO 


Iilall  type:  Roughi  3  meter  radius-of -curvature 
Uo  =  b  •  S  m/s 
ho  =  l3»0  mm 
k  =  D  •  33  mm 


p-trr=?^‘Jiaaarijnr-*-nr^r=?cGJ»-n-3'a"iJimmruD"mr-Lnr-jru 

ca^^trDrErtr^trcrir[rir^eo«^3cra"crcr^^^[r^cr£r^[r 


aaaaaaaaaaaaaDDoaoaaaaaaaaaaa 


crircrjir^rucrcrLnr-jjr^Lnojima-jaLnj-^a-rumaLnoiri-n 

^Coeocrcrir^i^cocrtrcrcrirtrcrcrcrcrcrcaeocrartrcroocr 


aoaaooaaaooaaoaaaoaaoDaooDr^aa 


j3P-LnrucroH=irurujLncrmircrnm3“LneOr=!Ln«oam3“eoi-na 
coeccocreoir'iFcoEr  ^[r"coc<][r[rcrcn^^ca^[rcnircp«ojF 


oaaaaaaaooaDDaaoanocDaoaaaaaaa 


mrucoomo^^-rqmmaimmLnnjr^r-^a^mrumamr^a 

rumrq^aaii-Ra-r=imr-qar^rur^iramj}r^ar^Lnomi-nmcj 

aaaaaaaaaaoooaaoaaaaaaaoaaaa 


oooooooooooooaoooooooooooooo 

II  I  I  I  I  I  I  I  I  t  I  II 


mmamiximeai^nji-RaamH?r=imjiiHmrumr^i-Ra[roram 

mrro^mruarrr^marumrnm^r=»m:j*r^njj3mammmm 

aaaaaaoaaaaaaaoaoaoaaoDaaaHia 


aaaaoaaaaoaaacjoaoaaaaaDaaoaaa 

II  III  I  I  II 


ruaa“Er^Lncri-=iamLnr“[rrujiLnrca!T"a[raLjncor-mr-=»i-Rr- 


aaoaaaaaaaaaaaaooaaoaaoaaoaaa 


j3jmr^j:racrrnDirr-=iruanjLncDaeOr^j-mruLnrrmr-uia3 

P*j}rurrcrruP“Ln^EOcrtr[rrur-irj3r-ar-^^r-»-nJr^i>ruKOr^ 


r-=1r=1r^r=lr^mrqnjr^r^arRr^rUr-=»r^r^r^runJr^f-qr^r^rUr=?nJrRrU 


i-^jitr^maaar-cGjr-eGivr-morrmruzrLnmrrzra-rrrua 


jjar-*ojmrut-ncammmLni-njaHi[rnRrRJimji3"ErLnnjirDr 

LncomJicrcor-r^mnjp-aLnBOLnLneoj3Eocrr-r-cOr^Ln*ar^Erj 

HR  HR 


amcotoDa-a-JETtor-i-RruHRtrnjr-aeomcrirr-mrutrmtRDJj 

mrurunjmmrunjrurunjmmmajnjrumrurnrururummrummrn 


aaoaaaaaaoaaaoaaaoaaaaooaaaaa 


p-mm:Ta^arrmuni-RLnmj[icrr^conji-najiaHRr-iraDair 

croeojtnJP-camcaLnirr-crirmruirHRcamr-accjeocrmm 


HRi-Rr^rUHRHRaar-RDr^Di-RarUi-RnRi-Rr-Rai-Ri-RruaT-RanJnRnJ 
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2ED-D  OEO 


b50-0-  ‘lEG 


t-ifcJUjt-'fc-'runjb-'h-'DDaujujab-‘b-*b-'b~*njafc-*b-^b-'rLib-'njtr*fc-'fc-«riJb-*a 

m  m  m  m  m  m  m 

□aCrUJjaUlOrUCDDaOa-D^aaJOsCnCrtH^tJLJ^IHJIb-'LnLUDairLnUJtrOjD 

UJJ3Crh-*0^-OtHt>LnLnJ3-cnj-irD3CPLn-C'J3^UJC3Lnafc-»b-'£na-C’a-O£r’4r 


aaaooaaaaaaaoaaaaaaaoaaaoaaaoaoaa 


ujijJLULUjr^ijummmujfc-‘uiijamruriJLiJtJLjujrunjLLiruLijnjruuj 

iHnjjTjLUb-'moaD^cr-ooa-Jh-'-cr-jji-dtnJ-a-canjcrruLnriJ 


-c  nj  m  uj  ru 

LU  ^  l>  tn 


fc-*  fcH  tH  fcr*  b-* 

-uoLna-OLn-vicr-vJOaajLncnjiocFcroa-ooaao-oji-orcrjinDacr-Ji 

b-*tnc>>a-c‘-oa3ax'Cn(rcro[r-cr-jjij3^ujLncrabJcrj3Jitnnjacn4r(r 


Ln^ir"fc-‘ru^4rnjrub-*D3LLia-oocrLn*jib^cFcr'^^-c-OLiJ-j-o-ob^D>a-c’ 


b-'lr-»bJ|Hfc-‘b-'h-'b-^fcHrUt-‘fc-‘rUb-*nJb-*lHb-*b-‘fc-'LUb-'t-*b^bJfc-*tJlHnJtHb-*fc-‘fc-‘ 

tnLnDatntn-C’-rOaLna-JOaarUb-'b-'O-lI-JOaaJlJJLnOaLJCnoaiF-OErJ] 

JO^oruDscnLurnjfcHnjjiijj-uaoaruLnDarujioaLnLjuatjjb-'-ob-'oa.Do 


aaaaoaaaaaaaaaaacziaaaaaaoaaaaanaaa 

LnLnLntr^-vJ-0-crtnjr[rru-J[rijJLnjrtnDrLnruUJLn-trLnUJUJtnLnLnJrjr-c 

ODjTcrxrcrrutzib-'LUOEJanjruLatnacrb-'jLnir'jrjr-aoair'jraib-'aDanj 


*  I  I  I  I .  I  I  I  I  II  I  I  I  I  I 

oaoDoaaaciaDaoDaaaaaoaaaaaoaoaaaoD 


aaoaaaacDaaaaaaaaanaaoaaaaaaaaaaa 
□  jTOjrcr^cDjroLnrunjru-cuno-c-jruLn-cijjaruujijjru^njcrnjnj 
D3a-Ni&Lij-c‘Cb>x’-c‘rLJD30aaDrutt*DaJ3ajacrjr^inru03fc-'03LU-c-j 


I  I  I  I  I  <  I  I  I  t  I  I  II  I  I  I  I  I  I  I  I  I 

aaaaaoaaaaaaaaoaaaaaaoaaaaaaaoaa 


aaaaoaaaaaaaaaaoaoaaaaaaaaaaoaaa 

a-c*D3irtnLniHUja-oruujnjLU-c:iHErtJiiHir-c’(jiCDbHUjain-obHcrb-»ci 

tt>-nLnJ3UJLn-CLnUJJtnLJfcH03C>*rUIHb-'[r-C'U-iC^ljLlUJ03jajr-JC3-Cirir 


aaoDaaaaaaaaaoaaaaaaaaaaoaaaaaooa 


D>D303J10a-OOiJCJsJJIJ10»DaD303jljaOa030i0303D3D>J3C330303030>03JD 

ruuJOab-*4r-ofc-‘x‘-ir^runj-0J0^J3fc-^aj]fc-'tnjr^[nLjJ4rtnLnttaLLlO3jr-j 


aaaaaaatHaaDaaaDaaaaaaaaaaacijaoaaaa 


Lunjrutr-jnj-rrajjrntHbJ-c-Daujrijnj-c^Daxr^nj-cujkHDaDaLnoaDaaj 


aaaaaaaaaaaaaaaaaaaaaaaaaaaoaaaaa 


lutrtnaLnaujb-'crx'aru-jtHJDODJikHLJin-jajajrujjitnaLn-oLnoa 
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Lf/ho  Df/ho  angle(deg)  ki  Ddp/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
2.7M  Q.B?  bfl.T  l.as  □•ST  D-OIS  a.DOT  Q.ab  Q.flfl  Q.ati 


Lf/ho  Df/ho  angle(deg)  ki  Dd^/ho  Vtip/Uo  Vo/Uo  Ubase/Uo  Utip/Uo  Udrop/Uo 
1.7T  Q.34  fl2.7  1.7'1  0-54  D-Qia  Q-DSt.  0-77  Q.7t.  0-77 


i-qcTLHJimLn^cinjrumjDJr^njmr^J 
MDr-eOeOeacOP-r-cOcOr-eOeOcOctJeOeG 
•  •■■•••••••****** 

aoaoacDoacaDaaaaooa 


rrcr:TLntrr=imi-nr-EanjJijar^trm 

eOP-rN.p-p-i^taor-<3ea^r-r"jir“ 


ooaanoaocDDaaaaaa 


rqmjicOr^meaeaa:3“Jir^mr^jairjamac:ar^irirr^«Do«oa-crrutrrur^ 

etjcaeoeoiosop-r-creop-sacaeacar^cocasap-^p-p-jicaeoeacar-cor^r^r^ 

aoaaoaaaaaaaaaaoaaooaoaooaoaaaaaa 


eajcotnr=irur-mrqHir^3-mcrnRirmmcrr-Lnarrrujuicammn-aieom 

r^r^cocoecieor-r-saeor-tfleor^eor-coeor-p-r-cacoeccor^^caeor-r^jt^ 


aoaaaaaaaaaaaaoooaQaaaaaaaaaaaDaa 


r-njeGi^r-[rr=?a«<3mr-=?mr^r^mirrun-mi-nr^car^zrrui-R^J3mnjmm 

ruLnm<a"eo*-nrnruaDHi^arumcr*-nr^i-=fr-R^r^r^r-rui-RrR^t^arum 

aoaaaaaaaaaaaaaaaaaoaacjaDaaaoaaa 


aaoaaaaoaaaaaaaaaoaacDtDaaaaaaaaaa 

till  I  I  I  II  I  I  II  I  ■  ■  ■  * 

3-aj]Lnr-[roD«tinjeo^j^^ir=orrjBonjaanjao3"-Dru3-cGKO 

jja-Lnr^eor-Druar^ainjamnjr-jar^amruru  r-  ru  ru  a  m  ui  o  a  a 
aaaaaaaaaaaoaaaaoaaaaaaaaaaaaaaa 


aaaaaaaaaaaoaaaaaaaaaaaaaaaaaaaaa 

,111  It  I  II  III  II  I  II  II 


rraamJieo=o[r=oi-njamjianjj]miri-qErrur=!aErcaa-i-qnjLnaicr^=o 

jcojijimru^^^^Tj-LnirjjiruLnLnuizrjja-uiirLn^^meairtnLn 

aaaaaaaaaaoaaaaaaaaaaaoaaaaaaaaaD 


j3runja-jp-£rammi-R*omiJiEOcr-iip-anjcrLnr-ir"*jir^mmcrnjEOir-D 

r^i-ni-=iJ3-3"=ap“jnjr-jj3rnjDJiErr"ap-Drr-r^[ni-n«om3“ruea^mcojr 

r^r^mr^riJr^r=1i-Rr^r^r^i-^rUi-Ri-^i-=lr^riJi-^r^r=irUr^r^r^r^ruriJi-Rr^rUi-^i-^ 


-3“zrLnm<3^r^jj-3“^^«or^^-=>’D'=^‘rLieoLnjacr^mr>--jjmp-p-=otc«o 


r^zrtam=amar^m«o(rr>-cojj3r^rujrurreoirzrr-mr^i-nErnjaairtr 
r-j]Lnj"trr-=?a»-nr^Lncococor--j3a^crp--u^tro"ncjJ3<JJnjjeojir- 
i-q  1-^  r-R  rR  r-R  HI 


anjH»cGi-n:3'^or-eajiDjruLnmr^Lni-Rmr^r“^r^j]-=rEr-iiHi3-Lnrna 

a-^jmmmmmnjrummmmmmrnmmmmrunjmrurnnjrumrrrnm^ 


oaaaaaoaaoaaooaoaaaaooaaooaaaoaoa 


cojir-«^H?r-mLnmajJ3r-ar-r“ruuieanjmj3truip-cri-n«43r^mH?eceGeo 

ajr^m^HiErco^eoacoruru«J3HimruLni-=ii-nr^ai-^crcomLnmmnjEaj 


i-qHiaHiruHiaHiHianjDrur^HiHJHii-RHRrijHiaH^i-^aHiHii-^i-RrijHii-RHi 
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bOQ-D-  9DD 


□  O.QOb  0.7^  a.7M  0.73 


aOruab-'lHb-*tHrUb-»rUb-»kHfc-‘nJb-'b-‘b-^tHrUb-*b-'tr‘fcr-‘t-«b-'brJrUb-'rUb-'hHnj 

jDJb-^j3jrt-*nj-crcrnjLnruaruatjj-c:ruarub-*fcr*oDaLntHruruLnru-irx: 

crrujjjTb-'c^aa-dOaLnLuijJCDfc-'jujLn-c'airuLLJirLnj^ruLntr'UJ-jEra 


aaaaaaaaaaaaaaaaaoaaaaaaaaaaaaaaa 


mixifLjjLuujLuijJLJLuunriJijjLuujLLJUjrijUJ^^tiJUJb--'njriJtjuujuJUjrux^LU 

Lnb-'Oatnainjnja-Ob-'nj^ijjiH-c-o-irmnjULimiHJDaJi^ma-cratJsrunj 


bH  b-«  fcHtHfcH  tHtjfcH 

aijDruj3-xj-Jcrnjcrjrubr'j3-vJcr-Njb-»rufc-'Lntrc>*J3-jaLn-vjDa-dD303tn 

cj>03-jnj03-u(nb-'oai03a-trru-jnjruujao3UJ-jao3rLi[rtHcr-cro3-ct-*nj 

jDmajaD"fc-‘b^jro3njb-*rub-^ruLnLncro3iJLiuj^b-‘^03rucnujujj]mtt>(rLn 


b-*mmb-^fc-»mb-*iiJb-'b-'fc-jmbHb-*rub-'njb-‘b-‘rijfc-'riJb-'runjfc-*rufc-'fc-'fc-jnjfc-'fc-' 

DrLnfc-*-j->Jijjj]CPjiD3(pac^ji-jj3acrjruo3UJ(rLLj&-ofc-*tp[pinLiJ-j-j 

Ln03[ro3Ln-j[n-ctH-josiji-jujLnajj]03aao3jnja-j-jujtniiJcrfLUD3Uj 


oaaDocDaoaaaaaoaaoaaaaaaaaoaaaaaaa 

-C'UJCTLnLnUJtnrULnCroaUJLn-cUJir-crUltntncrjrLnruruUJljJLncr-C-c-JLn 

Luirrucr-c’Ln-cc^runjtr-crErjicraiHUJUJLjj-uujjDDactJjD-jDaatnfcHfcJtr' 


*  *  I  (  I  III  III  I 

dooocDooooaaaaooooocDoaooaoaoacDoooo 


aaoaaofcHaaoaaaaaaaaaaaaaoaaoDDaaa 

ru-trcrab-'aoruaDa-Nj-n:oujujb-'03t-'ruLncjairtt>iHrua-cosruErnj 

oax^jTjrmoaaDijr-ccrxrjrcrrucr-cirarucu-cajrircDOa-cDDs-ca 


■  I  I  I  I  I  I  I  III  II  I 

aaaaaaaaaaoaaaaaanaaaaaaaaciaacjaa 


aaaaaaaaaoaaaaaaoaaaaaaaaanofc-'aaa 

Ln-ctnfc-»ab-»03b-'rucrinujbJnjrunj-jnjfc-'crnjb-'Drcraujnjtrnjmcrh-' 

b-'lJJOLnLnOUJIH-JDrLn-unjnj-C'-JnJ[rcrUJ-Ct-‘tnLn-C'UJ4rCnrU^LiJ03 


aaaaaaaaaaaaaoaaaaaaaacDocDaacDaaaa 


ujaaj]uja-catnaa[rtnj}ijjcD-cira^ijLjmoanjruo3[r^[ra-cuj 


aacDoaaaacDaaooaaaaaaoaaaDaaaoacDaa 


ruruoujLnnj-cLnujoaLn-a-crjajoaujirJUjaorLnir'ruajoaLnb-'irnj-Njj] 
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